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ABSTRACT 



A quantitative analysis was carried out on the performance of 
turboprop aircraft within a microburst windshear. The objective of 
the analysis was to provide specific flight procedures for optimal 
navigation through the windshear. The microburst windshear model 
used in the analysis embodied the severe characteristics of the 
microburst encountered by Delta Flight 191 during an approach to 
landing at Dallas/Ft. Worth, 2 August, 1985. Different escape 
strategies were tested using the flight performance characteristics of 
the U.S. Navy's P-3 “Orion" and T-44 "Pegasus" aircraft. The three 
flight phases investigated were approach to landing, takeoff, and the 
low altitude ASW mission. Results from the analysis were coupled 
with the pilot's view-point from which conclusions were drawn. The 
results of the analysis support a constant-pitch-angle escape 
procedure. The same procedural steps can be used for both aircraft 
in any configuration or situation with the difference being the degree 
of pitch to employ. The conclusions are in a format for integrating 
specific microburst escape procedures within the NATOPS programs 
for the P-3 and T-44. 
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I. INTRODUCTION 



The microburst windshear is documented as a serious hazard to 
aviation. The microburst has been a contributing factor in several 
airline accidents in the past two decades. The seriousness of these 
events has lead to a significant amount of research on the effect of 
microbursts upon airline transport aircraft. However, little has been 
documented on the effect of a microburst on llght-to-medium weight 
turbopropeller aircraft. The U.S. Navy operates a fleet of P-3 "Orion" 
turboprop aircraft in all-weather, long-range patrol missions. The 
P-3 aviation community requested a study of the effects of a 
microburst upon the aircraft's dynamic flight performance and the 
appropriate escape procedure if one were encountered. Training of 
the P-3 pilot starts in the T-44, a Beechcraft King Air twin turboprop. 
Therefore, the ultimate goal of this thesis was to derive viable 
microburst windshear escape procedures which can be Incorporated 
into the Naval Air Training and Operating Procedures Standardization 
(NATOPS) Program for both the P-3 and the T-44 aircraft. 

A. MICROBURSTS AND AVIATION 
1. Microburst Definition 

The microburst is a form of downburst windshear. 
Downbursts are associated with strong convective activity (namely 
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thunderstorms). Downbursts are formed by a cold air mass dropping 

from within a thunderstorm producing strong downdrafts. The 

downdraft turns outward as it approaches the earth. Fujita is a 

prominent pioneer in identifying and classifying this type of 

windshear [Ref. 1]. He defines the microburst as [Ref. 2:p. 8]: 

A small downburst with its outburst, damaging winds 
extending only 4 km (2.5 miles) or less. In spite of its small 
horizontal scale, an intense microburst could induce damaging 
winds as high as 75 m/sec (168 mph). 

Conclusions drawn from the Northern Illinois Meterological Research 

on Downbursts (NIMROD) and Joint Airport Weather Studies (JAWS) 

[Ref. 2] programs show that the microburst structure consists of a 

down flow with a structure of one or more vortex rings as 

illustrated in Figure 1. The microburst may be wet or dry; 

precipitation may or may not be present in the downdraft. The 

design of a microburst structure was not understood, or even 

accepted as a viable phenomenon prior to 1982. Through the 1980s, 

meterological researchers and airline accident investigators came to 

the conclusion that the microburst weather phenomenon is a 

significant hazard to aviation and has been a contributing factor in 

many airline accidents. The accumulation of the NIMROD and JAWS 

project results and Inflight data conclusively shows that ring vortex 

flow is prevalent in microburst windshear [Ref. 3]. 
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Microburst windshear is uniquely different from other types of 
windshear. The other prevalent type of windshear hazard to 

aviation is the thunderstorm gustfront [Ref. 4], A gustfront 
windshear is associated with massive settling of rain-cooled air 
resulting from thunderstorm activity. A sudden wind-shift with 
gusty conditions that precedes a thunderstorm is characteristic of a 
gustfront. The main difference between a gustfront and a 

microburst is the size of the area effected. A gustfront is widespread 
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Figure 1. Wind streamlines characteristic of a 
microburst. 



across the face of one or more thunderstorms. A microburst is small 
and isolated within convective activity. The local magnitude of the 
horizontal and vertical windshear is what matters to an aircraft 
penetrating convective activity. Although smaller, a microburst can 
be much more of a detriment to an aircraft's flight performance. 
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2. The Effect Of Microbursts On Aviation 



The aviation community has always been aware of the 
potential danger of windshear. Classification of different types was 
not addressed until 1976, when The Federal Aviation Administration 
(FAA) published an advisory circular discussing the effects of 
windshear on flight performance. It was updated in 1979 [Ref. 5] 
where the term "downburst" was used to describe a strong 
downdraft associated with the center of a thunderstorm. There was 
considerable debate on the hazard of the downburst. Many argued 
that the downdraft rapidly weakened as it approached the ground. 
The research by Fujita provided overwhelming evidence to the 
contrary, thereby defining the microburst. By the mid-eighties, the 
term "microburst" was readily being used by the aviation 
community. However, the potential effect of a microburst upon 
aircraft performance was not readily understood by aircrew. 

The danger of a microburst became very apparent following 
the crash of Delta Flight 191 (DAL 191) at Dallas/Ft. Worth, August 2, 
1985. The crash ensued following the penetration of a thunderstorm 
during the approach to landing. The aircraft involved was a 
Lockheed L-1011, a three engine heavy airline transport aircraft. The 
National Transportation Safety Board (NTSB) investigating the DAL 191 
accident published in their report the following [Ref. 6:p. 81]: 



4 



The National Transportation Safety Board determines the 
probable causes of the accident were the flight crew's decision 
to initiate and continue the approach into a cumulonimbus 
cloud which they observed to contain visible lightning; the lack 
of specific guidelines, procedures, and training for avoiding and 
escaping low-altitude windshear; and the lack of definitive, real- 
time windshear hazard information. This resulted in the 
aircraft's encounter at low altitude with a microburst-induced, 
severe windshear from a rapidly developing thunderstorm 
located on the final approach course. 

Several pertinent questions arise from this conclusion. Airlines 

routinely fly through "bad" weather every day, How do a flight 

crew and Air Traffic Control (ATC) authorities recognize microburst 

type weather? Can any aircraft navigate successfully through a 

microburst or is it beyond the performance capabilities of the 

average airliner? Is there a significant departure between normal 

piloting techniques and microburst escape procedures? A large 

amount of research has been done to answer these questions. 

However, the focus has been directed predominantly toward airline 

transport aircraft. 

The advent of Digital Flight Data Recorders (DFDR) provided 
conclusive evidence of the structure and the potential effects of a 
microburst windshear on aircraft. The DFDR records the last 30 
minutes of voice transcripts, as well as several aircraft performance 
parameters. Aircraft body axis angles and accelerations, airspeed, 
altitude, and time are a few of the recorded parameters. These 
records allow the reconstruction of the wind field vectors encountered 
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by the aircraft [Ref. 7]. DAL 191 was equipped with a DFDR. 
American Airline Flight 539 (AA 539) followed DAL 191 and penetrated 
the same convective weather 110 seconds later [Ref. 6] as depicted in 
Figure 2. This second aircraft was a McDonnell/Douglas MD-80, also 
equipped with a DFDR. The availability of digital flight data allowed 
a comprehensive database to be built on the wind field generated by 
that microburst. 
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Figure 2. Flight paths of the two airlines through the 
Dallas/Ft. Worth microburst. 

The center of the thunderstorm activity was located 
approximately 10,300 feet from the approach end of Runway 17L at 
Dallas/Ft. Worth Airport. The DFDR data show that DAL 191 initially 
encountered a headwind followed by a very abrupt tailwind [Ref. 7]. 
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A significant downward airflow also accompanied the horizontal 



windshear which resulted in a loss of airspeed accompanied by a 
sharp nose down pitch. A very high sink rate resulted. The aircraft 
exited the bottom of the thunderstorm approximately 100 feet above 
ground level (AGL). Ground Impact ensued with the aircraft in a 
slight nose up attitude and at a high airspeed [Ref 5], 

In retrospect, several previous airline accidents can be 
attributed to windshear of this nature. Fujita had analyzed several 
takeoff and landing airline accidents involving windshear prior to the 
DAL 191 crash [Ref. 2]. However, the Dallas/Ft. Worth accident 
showed the need to address the microburst hazard and gain a more 
thorough understanding of how to alleviate the danger. 

3. Reducing The Hazard To Aviation 

One approach to reducing the hazard of a microburst on 
aviation is recognition and avoidance. This solution is very difficult 
given available aviation weather avoidance equipment. The primary 
windshear equipment available today is the aircraft weather radar 
and ground based Low Level Windshear Alerting System (LLWAS). 
On-board weather radar provides a picture of the amount of 
convective activity ahead, but cannot recognize microburst activity. 
Ground based LLWAS provides warning when wind speed sensors 
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located about the airfield measure significantly different wind 



vectors. However, the microburst must occur at the location of the 
LLWAS system. An LLWAS system was installed at the 
Dallas/Ft. Worth airport at the time of the DAL 191 accident. The 
microburst occurred on the final landing approach corridor two miles 
from the center of the airport. Because of the location of the LLWAS 
sensors around the airport perimeter, the system did not give any 
windshear warning prior to the accident. An LLWAS alert was 
sounded as the thunderstorm passed overhead the airfield several 
seconds after DAL 191 crashed. Neither weather radar nor LLWAS 
provide real-time Information on microburst windshear. 

Very recent technology has produced some viable, and 
expensive, means to recognize and avoid related microburst weather 
activity. Ground based Doppler Radar has been successfully proven 
at Denver's Stapleton airport [Ref. 8]. On-board forward looking 
devices such as Infrared sensors, lasers, and doppler radar are being 
developed to provide several seconds of warning prior to a 
microburst encounter. But until avoidance systems are made readily 
available, survival of a microburst encounter will be dependent on 
the use of a successful escape procedure by the flight crew. 
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The second approach in reducing the microburst hazard is 



developing viable escape procedures. Research has shown that quick 
recognition of a microburst encounter is paramount. The reaction 
time of a flight crew coupled with the type of escape maneuver will 
govern the success of escape [Ref. 9]. The escape procedure the flight 
crew executes depends on the flight phase and configuration of the 
aircraft. Much research has been performed concerning aircraft 
performance and escape procedures with the focus on airline 
transport aircraft operations to be described below. 

B. MICROBURST ESCAPE PROCEDURES 

A large amount of research has been done on defining the optimal 
escape maneuvers for airline transport aircraft. Much of that 
research is based upon control theory. Comparisons of different flight 
strategies for microburst encounters have been addressed [Ref. 10, 
Ref. 11]. The total energy concept is readily used to compare the 
outcomes of different maneuvers [Ref. 12, Ref. 13]. One very 
important aspect in determining the optimal escape procedure is the 
ability of a flight crew to execute the appropriate maneuver utilizing 
the available flight performance information available. 
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The FAA has generated an exhaustive Windshear Training Aid 
[Ref. 14] aimed at modern day transport aircraft. It addresses the 
crew training requirements and suggests viable microburst escape 
strategies. Analytical research, as well as tried and tested 
procedures, have verified the suggestions presented by the FAA. All 
the major airlines have incorporated microburst recognition and 
escape into their recurrent training. The majority of the civil flight 
simulators possess some degree of windshear emulating a microburst. 
As of today, the aviation community recognizes the microburst as a 
aviation hazard and has taken significant steps in reducing the 
vulnerability. 

The different escape strategies fall within three general categories: 

1) Optimum aircraft performance through the airmass. 

2) Optimum aircraft performance related to an inertial reference. 

3) A combination of the two. 

Optimum aircraft performance is defined as the performance 
providing the largest cushion of escape from a microburst encounter. 
The parameters measured for evaluating an escape strategy's 
outcome, are airspeed and altitude gained (or lost). 

The first category includes maneuvers that trade altitude for 
airspeed. The second category uses measured data such as body axis 
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acceleration and windshear magnitude in calculating the optimal 
flight path [Ref. 12, Ref. 10, and Ref. 11], Obviously the second 
category requires significant airborne computing capabilities not 
normally available. The FAA's Windshear Training Aid recommends 
an escape strategy that falls within the third category. The 

recommendation is to maintain a constant pitch attitude with 
maximum engine thrust applied. In general, increasing pitch 
attitude toward a 15° pitch angle shows optimal performance on most 
of the large airline transport aircraft. The constant pitch attitude Is 
maintained with disregard to airspeed. 



C. EFFECTS OF A MICROBURST ON FLIGHT INSTRUMENTS 

How a flight crew perceives and navigates through a microburst 
encounter is strongly influenced by their flight instrument 
indications. The primary flight Instruments available to pilots are: 

1) Attitude Indicator (AI) - a gyroscopic/inertial stabilized aircraft 
attitude indicator. Modern design allows accurate attitude 
information in the highest levels of turbulence. 

2) Airspeed Indicator - a pitot/static system instrument which 
displays airspeed in knots. It is sensitive to airmass pressure 
changes and position error. 

3) Pressure Altimeter - a static system instrument which relates 
altitude, in feet, to static pressure. It is sensitive to airmass 
pressure changes. 
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4) Vertical Speed Indicator (VSI) - a static system Instrument 
which measures the change in static pressure in feet per 
minute. It is unreliable in turbulence due to mechanical and 
pressure lag. 

All aircraft have these primary instruments. Inertial navigation and 
true airspeed computers do provide other flight performance 
references usually through flight directors. However, the majority of 
the flight stations rely primarily on the basic flight instruments. 

Microbursts are always associated with some degree of turbulence. 
Also, pressure differentials can be expected within convective 
activity. Ground and airborne data show that the pressure within 
a microburst varies about ± 3 millibar [Ref. 2, Ref. 3], This pressure 
change equates to about ± 80 feet in altimeter variation. 

The attitude indicator will be the flight instrument least affected 
by microburst atmospherics. This instrument will always give 
reliable aircraft attitude information. Due to the minor static 

pressure changes, the airspeed indicator and altimeter can still be 
relied upon as performance instruments. The airspeed indicator will 
be one of the primary instruments to indicate a windshear 
penetration. The VSI would probably be erratic and unreliable due 
to turbulence and instrument lag. 
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There exits a means of measuring angle-of-attack (AOA) in most 



larger aircraft. Transport aircraft utilize AOA as a stall warning 
device (stick shaker). Navy aircraft incorporate a means to read AOA 
in units. AOA measuring devices may have a high damping factor in 
their measurement, resulting in substantial lag during turbulent 
conditions. Therefore, AOA indication may not be a reliable 
performance indictor during a microburst encounter. 

D. PROBLEM STATEMENT 

The P-3 has no means to recognize and then avoid microburst 
windshear. It must rely solely on escape. The published guidelines 
for microburst escape are based upon airline transport aircraft which 
significantly differ from the P-3. In this study, the results of 

different escape procedures based upon the available flight 
instruments are compared using the parameters of altitude, airspeed, 
and specific energy. Characterizing the effect of a microburst 
windshear upon the aircraft’s performance followed by the 
implementation of the appropriate escape procedure form the 
objective for the research. The final conclusions must provide viable 
escape procedures that are based upon the primary flight 
instruments. 
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The microburst windshear model selected must conform to a 



known, measured phenomenon. The windshear experienced by DAL 
191 provided a suitable database for emulation. 

A series of equations of motion for the aircraft were developed 
and designed to be controlled through pitch angle and thrust inputs. 
These two parameters emulate the pilot’s available flight control 
inputs. 

Three flight phases critical to a microburst encounter were 
considered. Approach to landing, takeoff, and on-statlon (P-3 only) 
phases require operation low to the ground in all-weather conditions. 
Viable escape procedures using attitude, airspeed, altitude, and angle 
of attack were analyzed for all three flight phases. 

Three questions were considered for each flight phase: 

1. What is the optimum microburst escape procedure given 
available flight information to the pilot? 

2. Does the optimum escape procedure change with gross weight 
or available thrust? 

3. What flight instrument indications, flight control feeling, and 
dynamic response would be expected during the optimum escape 
maneuver? 

Viable microburst escape procedures were derived for the P-3 and 
T-44 aircraft. 
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II. MATH MODELS AND EXPERIMENTAL DESIGN 



This analysis is based upon mathematical representations of 
microburst windshear and aircraft performance. Microburst 
windshear, inertial reference, aircraft performance, and critical flight 
phases form the major segments of the total algorithm suite. Most 
of the mathematical simulation performance was carried out on a 
Macintosh SE® computer utilizing MatLab® software. 

A. MICROBURST WINDSHEAR MODEL 

A double ring vortex was chosen to emulate microburst 
windshear. Fujita [Ref. 2:p. 14] demonstrated that the microburst 
was a settling airmass with vortices generated near the earth. This 
phenomenon produces varying degrees of three-dimensional 
windshear when close to the earth's surface. Figure 1 showed the 
general concept of a microburst. 

A strong microburst windshear was experienced at Dallas/Ft. 
Worth on August 2, 1985. A Delta Airlines L-1011 (DAL 191) crashed 
while on final approach during a microburst windshear encounter. 
An American Airlines MD-80 (AA 539) performed a missed approach 
110 seconds following the DAL 191 crash. Both aircraft were equipped 
with DFDRs which register a time history of the aircraft’s 
parameters. Body axis accelerations, velocities, and Euler-angle 
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values were used to calculate the flow-field winds generated from 
the microburst. Wingrove and Bach [Ref 7] provided the analytical 
means to calculate the windshear. From these data, an insight was 
gained of the magnitude and characteristics of the windshear. 

A double-vortex-ring model was chosen for the windshear model. 
Schultz [Ref. 15] devised a double-vortex-ring mathematical model 
which closely approximated the windshear experienced by the AA 539 
flight. The aircraft flew through the microburst windshear at 2500 
feet AGL. The original algorithm was developed and tested against 
recorded flight data. The double-vortex-ring algorithm was modified 
in the current study with the addition of source flows to better fit 
low altitude windshear. This modified windshear model was fit to 
the DAL 191 flight data by the application of a parameter sensitivity 
scheme. Ring location and vortex intensity were found to differ 
from the AA 539 model. However, a multiple microburst structure 
can be expected during strong convective activity [Ref. 2:p. 35]. 
MatLab® programing was utilized in this study and is listed in 
Appendix A. 

1. Vortex Ring Model Development 

The development of this vortex ring model comes directly from 
the work of Schultz [Ref. 15]. Source flow was integrated into this 
algorithm to more closely match the DAL 191 windshear. The wind 
component in the y direction was not considered in the current 
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model. The y component Is not easily modeled by symmetric vortex 
rings [Ref. 15] and has no direct impact on the longitudinal dynamics 
of an aircraft. Figure 3 illustrates the geometric aspects of the 
vortex ring model. Imaging of the vortex rings leads to inviscid 
ground effect upon the wind field. The wind field becomes horizontal 




Figure 3. Vortex ring mathematical representation, 
next to the ground plane as seen In nature. Boundary layer effects 
were neglected in this model. Bowles and Oseguera have performed 
research in boundary layer effects on windshear [Ref. 16], From 
Figure 3, the following relationships can be seen to hold: 
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where G, P, and I are ground, primary, and image reference 
respectively. The velocity components for each vortex were obtained 
by differentiating the vortex ring stream function 
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where r is the radial distance from the z axis: 
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The stream function for a vortex ring filament was obtained 
from the evaluation of an elliptical Integral [Ref. 15]: 



0.788* 



••SW 7 2 

0.25+ 0.75V 1-* 



(6) 



where r\ and T 2 represent the closest and farthest distances to the 
point of interest from the ring's filament; k is the ring filament 
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strength; and X is a scaling term. Their algebraic relationships are as 
follows: 



r^Vz^ir- 
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where R is the radius of the filament. 
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Equation 6 was differentiated to obtain expressions for 3r and 
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( 16 ) 
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Equations 10 and 11 were substituted into Equation 4 to get the 
Individual velocity components produced by each ring at a given 
point in space (with respect to the ground reference). 

A velocity damping factor (£) was calculated for each ring to 
prevent erroneously high values of velocity near the filament cores. 
The value ranged from zero to one; it is zero when the point of 
interest is at the viscous core and approaches one at Increased 
distances. The algebraic relationship for C is as follows [Ref. 15]: 



where ru is the closest distance from the point of interest to the / th 

ring's viscous core. Schultz obtained the weighting factors, S and p, 
by a visual comparison of the vortex ring strength and through a 
parameter estimation scheme [Ref. 15]. A total damping factor (Z) is 
obtained for the point of Interest from the product of the damping 
factors for each of the four vortex rings: 




(19) 
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The low altitude windshear model required the addition of two 
point source flows located 10,000 feet above and below the ground 
plane at the center of the microburst. The radial flow is represented 
by the following expression: 



r 2 * R ! (21) 
where A is the source strength and Rj is the radial distance from the 

point source to the point of interest. The horizontal and vertical 
wind vectors contributions are estimated by: 
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Combining the upper and lower (mirror Image) sources produce: 
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(23) 



In summary, wind velocity for a point of interest in space 

referenced to earth was calculated in the following manner: 

1) X (horizontal) and z (altitude) positions, measured in feet, 
are inputs for Equation 4. Calculations are made for the 
four vortex rings (large primary, small primary, large 
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image, and small image). A horizontal wind velocity (vxj) 
and a vertical wind velocity (vzj) are calculated for each 
ring. 

2) A damping factor (^) was calculated for each ring using 
Equation 19. 



3) The induced velocities for each ring were added and 
multiplied by the total velocity damping factor (Z): 



Wx.flt, 



1-1 u-l 



Xvx, 



(24) 
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where Wx and Wh are windshear values (ft/s). 

4) For the low altitude model, the source flow velocities were 
added to the vortex velocities to give the final horizontal 
and vertical windshear components. 

2. Windshear Model Fit 

The unmodified double vortex windshear model was applied to 
AA 539 flight data by Schultz, using a parameter minimization 
scheme to fit the double- vortex-ring model to the actual DFDR data. 
Table 1 lists the parameters varied and the final results. Schultz 
concluded that the difference between the model results and recorded 
wind data was an rms (root mean square) of 16 ft/s for the entire 
system. Figure 4 shows the comparison of recorded winds to the 
model prediction. 
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The AA 539 windshear model was fit in the current study for 
model verification only. The windshear that DAL 191 encountered 
was the windshear to emulate. This windshear model was used to 
TABLE 1. 2500 FEET ALTITUDE WINDSHEAR MODEL PARAMETERS 



parameter large ring small ring 



ring radius, ft 


8503.3 


1701.7 


core radius, ft 


20004.1 


323.9 


ring circulation, ft'2/s 


431968.8 


57204.9 


x position, ft 


0.0 


50.0 


y position, ft 


3350.4 


830.9 


z position, ft 


3400.6 


2333.6 



analyze aircraft performance during takeoff and landing flight 
segments. Therefore, the double- vortex-ring model was modified by 
th parameters listed in Table 2. The data from DAL 191 showed that 
TABLE 2. LOW ALTITUDE WINDSHEAR PARAMETERS 



parameter large ring small rinR 



ring radius, ft 


7000 


1300 


core radius, ft 


2004.1 


323.9 


ring circulation, fV2/s 


431968.8 


131571.3 


x position, ft 


2500.0 


300.0 


v position, ft 


-300.0 


1.0 


z position, ft 


3400.6 


800.0 


x dir. source strength 


1355396 ft"2/s 


z dir. source strength 


3049200 ft~2/s 



the two rings were not concentric about the center [Ref. 3]. The 
small vortex ring was lower to the ground and displaced in the x 
direction from the large ring. Point source flows located 10,000 feet 
above and below the surface plane were added to the double-vortex- 
ring model to increase the outflow. A simple parameter sensitivity 
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Vertical Vind ft/s Headwind (ft/s) 





Disi from MB Center (ft) 



Figure 4. Vortex ring model comparison to recorded 
flight data [Ref. 3]. 
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scheme was then used to estimate the parameters that best fit the 
modified windshear model to the recorded DAL 191 flight data. 

The sensitivity scheme (listed in Appendix B) used an iterative 
process of varying each parameter individually, then calculating the 
x direction, z direction, and total system rms difference between 
model and recorded winds. The initial parameter matrix was built 
up as an 8x3 matrix. The center column was composed of the "best 
guess" values. Columns one and three were values obtained from an 
interval surrounding the initial guess value. First guess parameters, 
including the corresponding intervals, were obtained by graphical 
comparison of model and recorded data. The parameter matrix was 
refined after each successive run of the sensitivity scheme. Final 
parameter resolution was less than Z% change for the total system 
rms. Table 2 lists the best-fit parameters to the DAL 191 flight 
windshear measurements. Figure 5 is a graphical representation of 
the current model’s vortex rings and source flow point relative 
locations. The wind vectors produced from the model are shown in 
Figure 6. Figure 7 compares the model windshear with the 
measured wind velocities from DAL 191. 

Total system rms was calculated as 16.5 ft/s for the microburst 
windshear model. The final model shows a close approximation of an 
actual microburst, serving as a good windshear model for testing 
aircraft performance during low-level encounters. 
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Figure 5. Geometric relationship of the vortex rings and 
point source within the microburst model. 
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Figure 6. Vector plots of the wind-fields generated by 
the microburst windshear model. 
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distance from MB center (ft) 





Distance from MB Center (ft) 

Figure 7. DAL 191 recorded winds compared to 
mathematical model. 
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B. INERTIAL REFERENCE FRAME 



The concept of aircraft specific energy was used to compare the 
outcomes of different windshear escape strategies. The aircraft's 
total energy is defined as the sum of the airmass-relative kinetic 
energy and the inertial potential energy. A vectorial relationship of 
an aircraft's motion through the airmass is presented in Figure 8. A 

A h 




x 

1 ► 

Figure 8. Inertial reference coordinate system, 
set of coupled, non-linear differential equations can be developed 
from this vectorial relationship. These differential equations can be 
solved numerically where specific energy can be determined at any 
point in inertial space. 
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1. Total Energy Concept 



The aircraft's altitude and airspeed at any given point can be 
used to determined total energy from the following relationship: 

E = H>mV2 + mgh (26) 



where m is mass, V is airspeed, and h is the altitude. The specific 
energy (sometimes called energy height), Es, is defined by: 



F V 2 h 

Es " "2g" + h 



(27) 



The time rate of change of specific energy is equal to aircraft 
acceleration plus rate of climb. Differentiating Equation (27) with 
respect to time gives specific power: 



Es« 




+ h 



(28) 



Values for V, V, h, and h can be obtained at any Instant of time 
from the relevant equations of motion. 

2. Equations of Motion 

The relationship of an aircraft's motion through a moving 
airmass relative to earth leads to the development of six coupled, 
non-linear equations [Ref. 12]. Using the coordinate system in 
Figure 8, the following equations of motion evolve: 

x»VcosY a +Wx ( 29 ) 

hoVsin^+Wh ( 29 ) 
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( 30 ) 
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3x 3h 3t 



V * -J-cosa- - g sin-y - W x cosy -Wh sin-* 
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(32) 



T L g Wx^ W h 

Y " HV Slna + W V C0SV a + — Sln > a - — cos ^ a ( 33 ) 



3Wx 



mV V 
3Wh 



Note that 3t and 3t were set equal to zero due to steady state 
conditions being assumed for the windshear model. This assumption 
is valid due to the minimal time the airplane is exposed to the 
microburst (30 to 60 sec.). Note that the equations are for point 
mass analysis. No dynamic cross-coupling is considered. 

The equations of motion were solved using a predictor-corrector 
numerical scheme. Euler First-Forward, Euler Half-Step and 
Richardson Extrapolation [Ref. 17] were combined in the following 
manner: 



' , !l|- y n* 4ly n 


(34) 


y (1> .y ♦ — y 
y n + l/2 2 y n 


(35) 
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y n+l ^n+1 y n+l 


(37) 
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Second-order accuracy Is expected from this numerical scheme. By 
providing input values for aircraft constants and initial values for 
the inertial reference, the equations of motion were solved to obtain 
the aircraft performance. 



C. AIRCRAFT PERFORMANCE MODEL 

A MatLab® program was written to take aircraft performance 
parameters and use the set of motion equations to measure aircraft 
response in the windshear. Initial values were determined from 
aircraft performance parameters and the initial position. Iteration of 
the set of motion equations was based upon the time interval 
chosen. Values for x, h, V, v, a, 6, yj, y a , E s> and Es were 

tabulated after each iteration. The control of the model was taken 
from a pilot's point of view. The controlling variables were 6 (deck 
angle or pitch angle) and T (thrust). Both variables can be changed 
between iterations. 

Certain aircraft performance parameters were provided for 
calculating lift, drag, and AO A. The parameters were Clq, Cl«, Cd 0 > S, 

K, W, maximum thrust available, and maximum AOA. Cl and Cd 



were calculated as follows: 



C ■ C + C tt 

L Lo Let 



C * C + K(C) 
D Do L 



( 38 ) 

( 39 ) 
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Note that the lift-curve slope was referenced from the longitudinal 
axis of the aircraft (see Figure 8). Lift and drag were calculated 
using the familiar relationships: 



l*qsc l 



(40) 



d»qsc d 



(41) 



where the dynamic pressure, Q, is based upon equivalent airspeed: 



Qb 



1 \r 2 

2 p o V 



( 42 ) 



Initial values for airspeed, thrust, and theta were chosen to 
match the aircraft requirements for the phase of flight under 
scrutiny. Initial values for x and h were chosen referenced to the 
microburst windshear center. The iteration of the model was based 
upon the time step, At. 

Three aircraft were analyzed for their response in a microburst 
windshear. They were the U.S. Navy P-3 (Lockheed L-188), the U.S. 
Navy T-44 (Beechcraft King Air H-90), and a generic 3-engine “heavy" 
airline transport. The P-3 is a four-engine turboprop with gross 
weights in the medium range (75,000 to 135,000 lbs.). The T-44 is 
powered by two turboprop engines and falls within the category of a 
"light- twin" transport. The 3-engine heavy airline transport is 
powered by turbofans and as the name implies, is capable of high 
gross weights. Tables 3, 4, and 5 delineate the performance 
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parameters of the three aircraft. The lift curve slope for the P-3 is 
greater than the theoretical prediction. This result is due to the 
influence of the propeller induced flow-field. (Power-off Cl«= 5.7 rad-i.) 



TABLE 3: P-3 PERFORMANCE PARAMETERS [Ref. 18]. 



Performance 


Landing 

Approach 


Takeoff 


On-station 


Parameters 


Configuration 


conflKuratfon 


Configuration 


W (lbs) 


89,500; 114,000 


89,500; 120,000; 


120,000 


S (iV2) 


1300 


135,000 

1300 


1300 


K 


0.05041 


0.05041 


0.05041 


Cpo 4 engine 


0.0567 


0.0551 


0.0213 


3 engine 


0.0630 


- 


- 


C L a (rad-0 


6.38 


6.38 


6.25 


Clo 


0.800 


0.800 


0.350 


a max (rad) 


0.244 


0.244 


0.209 


TmaxObf) 4 eng. 


33400 


33400 


33400 


3 eng. 


25050 


- 


- 


q (rad/s) 


0.0873 


0.0873 


0.0873 


Vref (ft/s) 


236; 262 


- 


354 


VI (ft/s) 


- 


204; 214; 229 


- 


V2 (ft/s) 


- 


220; 227; 239 


- 



As mentioned, the AOA used in the calculations was based upon the 
longitudinal axis of the aircraft. Maximum thrust available and 
gross weight were adjusted to meet the scenario requirements. 
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TABLE 4. 3 -ENGINE HEAVY AIRLINE TRANSPORT 
PERFORMANCE PARAMETERS [Ref. 10]. 





Landing 




Performance 


Approach 


Takeoff 


Parameters 


Configuration 


Configuration 


W (lbs) 


362,000 


462,000 


S (ft*2) 


4578 


4578 


K 


0.059 


0.059 


Cdo 


0.108 


0.098 


Cl <2 (racH) 


4.96 


4.96 


Clo 


0.532 


0.532 


« max (rad) 


0.317 


0.314 


Tmax(lkf) 


126,000 


126,000 


q (rad/s) 


0.0873 


0.0873 


Vref (ft/s) 


227 


- 


VI (ft/s) 


- 


238 


V2 (ft/s) 


- 


255 



TABLE 5: 


T -44 PERFORMANCE PARAMETERS [Ref. 19]. 
Landing 


Performance 


Approach 


Takeoff 


Parameters 


Configuration 


Configuration 


W (lbs) 


8,280 


8,280 


S (ft A 2) 


210 


210 


K 


0.0503 


0.040 


Cdo 


0.120 


0.100 


Clo (rad-i) 


6.24 


6.24 


Clo 


0.587 


0.0523 


« max (rad) 


0.244 


0.227 


Tmax(lbf) 


3,023 


3,023 


q (rad/s) 


0.0873 


0.0873 


Vref (ft/s) 


203 


- 


VI (ft/s) 


- 


152 


V2 (ft/s) 


- 


202 
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For each Iteration of At, the following aircraft, Inertial, and 



windshear model values were calculated and recorded: 

1) x - Distance from microburst center in feet. 

2) h - Absolute altitude in feet. 

3) V - Equivalent airspeed in knots. 

4) 6 (theta) - Pitch angle in degrees. 

5) a (alpha) - Angle-of-attack in degrees. 

6) y a - Flight path in degrees referenced to the airmass. 

7) Yi " Flight path in degrees referenced to the earth. 

8) V g - Ground speed in knots. 

9) ROC - Rate of climb in ft/sec. referenced to the earth. 

10) Thrust - instantaneous thrust in lbf. 

11) Es - Specific energy in feet. 

12) Es (Esdot) - Time rate change of specific energy in ft/sec. 

13) Wx - Horizontal wind speed in ft/sec. 

14) Wh - Vertical wind speed in ft/sec. 

From the collection of data, plots were made of altitude, theta, alpha, 
airspeed, specific energy, and time with comparison to distance from 
microburst center. 
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Four basic escape maneuvers were examined for critical phase of 



flight. Constant airspeed, constant altitude, constant pitch angle (6), 
and constant angle-of-attack (a) escape maneuvers were evaluated 
quantitatively (engineering standpoint) and qualitatively (piloting 
standpoint). 

D. CRITICAL FLIGHT PHASES MODELED 

There are three phases of flight in which a microburst encounter 
becomes critical for airplane survival. Approach to landing, initial 
takeoff and climb out, and low-altitude maneuvering are the three 
phases of flight considered in this analysis. All three phases occur at 
an altitude of less than 1000 feet AGL. The low-altitude 
maneuvering flight phase is applicable to the P-3 community. The 
P-3 routinely operates 200 feet above the water in all weather 
conditions during Anti-Submarine Warfare (ASW) missions. 

1. Landing Approach 

The landing approach scenario is based upon a 3 degree 
glideslope descent to landing. In this scenario, the microburst center 
is placed 10,300 feet from the end of runway. The simulation starts 
500 feet (x=-500 feet) before the microburst center on glideslope 
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(h=566 feet). This situation closely represents the scenario of DAL 191. 



The aircraft is exposed to the windshear immediately at start, t=0. 
The simulation is run for a programmed length of time or until 
ground impact occurs. Five aircraft/weight combinations were 
analyzed for the approach to landing scenario. They were a P-3 at 
89,500 pounds gross weight, a P-3 at 114,000 pounds gross weight, a 3- 
englne heavy airline transport, a T-44, and a P-3 at 89,500 pounds 
gross weight with one engine shut down. 

Starting aircraft parameters were such that the given pitch 
angle (theta) and thrust would maintain a 3° glideslope at target 
approach airspeed (Vref). Pitch angle and thrust were maintained 
until a loss of airspeed equated to Vref minus 20 knots. At such 
time, the aircraft performed one of the following programmed escape 
procedures: 

1) constant airspeed - Maximum thrust was applied and a pitch 
angle set to 0°. This pitch angle was maintained until the 
airspeed equaled Vref. Pitch angle was then adjusted to 
maintain the airspeed at Vref + 5 knots. 

2) constant altitude - Maximum thrust was applied and a pitch 
angle was set to obtain a positive rate of climb. This pitch 
angle was maintained until the target altitude (altitude at 
which the maneuver began) was established. Pitch angle was 
then adjusted to maintain the target altitude by ± 20 feet. 
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3) constant theta - Maximum thrust was applied and a pitch 
angle (theta) was set and maintained. Specifically, theta 
values of 5°, 8°, 10°, and 15° were used. 

4) constant alpha - Maximum thrust was applied and a pitch 
angle was adjusted to get a given angle-of-attack (alpha). AOA 
values of 12, 15, and 20 units were used for the P-3 model. 

The above escape maneuvers were constrained by certain limits. 

Pitch rate (q) was set at 5°/sec. Thrust application rate for the P-3 

model was set at 0.5 maximum thrust/sec. Thrust application rate 

for the generic 3-engine heavy and T-44 aircraft was set at 0.2 

maximum thrust/sec. This rate accounts for engine "spool-up" time. 

Maximum pitch angle was limited not to exceed maximum AOA. 

Appendix C contains the MatLab® program used for the approach to 

landing scenario. 

2. Takeoff 

The takeoff scenarios primarily explored the effects of 
microburst windshear on the takeoff performance with penetration 
at liftoff. However, some analysis was performed on the microburst 
center distance from liftoff point. Generally, the simulation began 
with the aircraft lifting off 1200 feet from the microburst center 
(x=-1200 feet) at the appropriate liftoff speed (Vj). Initial pitch angle, 



39 



theta, was such to achieve takeoff safety speed (V 2 ) at 50 feet, no 



windshear. Maximum thrust was used for all cases. Execution of 
the particular escape maneuver was begun when the rate of climb 
(ROC) is less than or equal to zero or the airspeed fell below V 2 . Five 
aircraft/ weight combinations were looked at. They were a P-3 at 
90,000 pounds, 114,000 pounds, and 135,000 pounds gross weight, as 
well as a 3-engine heavy airline transport and a T-44. 

Four escape methods were considered for a microburst 
encounter at takeoff. When the ROC was less than or equal to zero, 
or the airspeed fell below V 2 , one of the following programs was 
executed: 

1) constant airspeed - If airspeed was less than V 2 at initiation, 
theta was reduced to 0°. If airspeed was greater than V 2 at 
initiation, theta was increased. In both cases, theta was 
manipulated to maintain airspeed = V 2 + 5 knots once V 2 was 
achieved. 

2) constant altitude - Theta was varied to maintain altitude ± 20 
feet about the target altitude. The target altitude was the 
altitude at which the escape maneuver began. 

3) constant theta - Theta was held at a programmed constant 
value throughout the maneuver. Theta values used were 5°, 
10°, and 15°. 

4) constant alpha - Theta was varied to maintain a constant 
AOA value. 
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For the above maneuvers, a pitch rate of 5°/sec was used and 



maximum thrust was maintained. Also, theta was reduced any 
time critical AOA was exceeded. 

Two variants of the takeoff scenario were analyzed for 
performance sensitivity. Both used the 120,000 pound P-3. For one, 
V[ 0 f was increased 18 knots to 145 knots. This equated to an increase 
of the rotate speed (Vr) to 145 knots. The second variant moved the 
microburst center from 1200 feet to 2000 feet from the point of 
liftoff. This resulted in the aircraft gaining airspeed and altitude 
before encountering the severe horizontal windshear. The second 
scenario closely simulated a microburst encounter after takeoff. 
MatLab® programming for the takeoff scenario is listed in 
Appendix D. 

3. P-3 On-Station 

The on-station microburst encounter scenario used a P-3 at 
120,000 pounds gross weight in a 4-engine loiter configuration. No 
wing flaps or landing gear were extended. Only aircraft reaction and 
performance were analyzed. No specific escape maneuver was used. 
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This scenario tried to emulate an autopilot maintaining altitude 



through pitch authority and the flight crew controlling thrust. 

The microburst encounter began with the P-3 at 200 feet AGL 
and 5000 feet from the microburst center (x=-5000 feet). The initial 
airspeed was 210 knots (prescribed loiter airspeed). Theta was varied 
throughout the encounter to maintain altitude. Thrust was 
maintained at loiter power until 40 knots of airspeed was lost. At 
this point, maximum thrust was applied while still maintaining 
altitude. The scenario was ended at 5000 feet on the opposite side 
of the microburst center. 
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III. RESULTS AND DISCUSSION 



The microburst and inertial reference math models were 
validated by comparing actual flight data to model data. The 
inertial model was further tested for the stability of the coupled 
differential equation scheme. Any weakness noted was considered 
when the results were analyzed. Each critical phase of flight was 
studied for aircraft response when performing a particular escape 
maneuver. Mathematical results were combined with other 
observations to build a foundation for conclusions. 

A. WINDSHEAR AND INERTIAL REFERENCE MODEL VALIDATION 

The windshear math model was developed from fitting a 
vortex/source flowfield to recorded flight data. The wind field 
recorded by the AA 539 DFDR (Digital Flight Data Recorder) showed a 
definite vortex flow. However, the flight data from the DAL 191 DFDR 
showed a different vortex flow arrangement. In addition, a strong 
outflow at the surface, and Increased vertical sink, required the 
addition of source flows to the windshear model. This use of point 
source flows led to a much closer fit of the original vortex model to 
the wind field DAL 191 experienced. The source flows have no range 
damping terms. Therefore, the DAL 191 emulation windshear model 
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becomes invalid at distances greater than 6000 feet (x=6000 feet) 
from the microburst center. All analysis was easily done within this 
distance limit. 

It is important to note that the aircraft performance results are 
insensitive to the exact modeling of a particular windshear. The 
rapid change of a headwind to a tailwind is the governing factor In a 
microburst. Secondary are the vertical down drafts that can be 
experienced above 100 feet AGL (above ground level). The windshear 
math model contains both characteristics to the same degree as the 
windshear experienced by DAL 191. 

Validation of the inertial reference model was scrutinized for 
proper aircraft response to changing conditions and the effect of the 
time step size. Stability of the solutions obtained from the 
differential equation numerical scheme was the greatest concern. 
Results generated from the inertial model were studied for the light 
P-3 under stable and turbulent conditions. Also, a comparison of the 
3-engine transport model was made with the DAL 191 DFDR data 
during the final seconds of the fateful flight. 

The first approach was to check the response of a light P-3 
initially stabilized on a 3° descent path with no windshear. The input 
came from a subroutine that would not vary pitch angle (theta) or 
thrust. No excursions were noted in descent path, airspeed, or 
alpha, as shown in Figure 9. The time step used was 1 second. 
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The second validation was to observe the stability of the inertial 
model solution with different size time steps. A time interval of 0.5 
seconds was the target time step for running all analyses. Therefore, 
time steps of 1, 0.5, and 0.1 second were investigated. The light P-3 
on approach to landing was again used. This time, the windshear 




x<ft> 

Figure 9. Light P-3 on a stabilized descent path, 
model was incorporated to induce changes in the flight environment. 
Aircraft control was through theta and thrust inputs. Control 
feedback was set to vary theta to keep within ± 100 feet of the 
descent path and to vary thrust to maintain V re f within + 5 knots. 

The results are depicted in Figure 10. Note the close correlation 
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between the 0.1 and 0.5 second time step cases, indicating that a 0.5 
second time step provides stability and adequate accuracy for the 
desired results. 




Figure 10. The effect of time step size on the inertial 
reference calculations. 

The final validation was comparing the actual aircraft response to 
model predictions. DAL 191 DFDR data were used to compare an L-1011 
flight path, airspeed, and alpha data to that of the 3-engine heavy 
airline transport model. Inputs to the aircraft control were the 
recorded theta and thrust of DAL 191. The windshear experienced by 
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DAL 191 was also Incorporated for this comparison. Figure 11 
graphically displays the close characterization of the actual flight 
path by the model. Airspeed and alpha do not reflect actual aircraft 




Figure 11. Comparison of aircraft response in windshear 
between the DAL 191 Flight L-1011 and the 3-englne 
transport model. 

response as closely as the P-3 model, but do follow the same general 
response. The lag in airspeed and alpha is probably due to the model 
Cl and Cd equations being only of second order. 
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B. WEAKNESSES AND OTHER INDEPENDENT VARIABLES 

Some weaknesses in the math models and unobserved 
independent variables exist in the total analysis. The most notable 
are: 

1) Cl and Cd value errors at angle of attack near stall. 

2) Unobserved effects of dynamic pressure and AOA upon engine 
thrust. 

3) Unobserved effects of rain and turbulence on aircraft 
aerodynamic performance. 

4) Unobserved effects of coupled longitudinal and lateral dynamic 
modes excited by turbulence and pilot induced oscillations. 

None of the weaknesses or unobserved variables is believed to 

significantly impact the results of the analysis. The overall concept 

was to compare outcomes of different microburst escape strategies 

given the same parameters. 

A weakness in the aircraft model is the ability to predict the 
effects of flow separation at low airspeed and high angle of attack. 
These effects would be an increase in drag and a decrease in 
generated lift. Table look-up or higher order equations could be used 
for the Cp and Cl expressions. The model does limit the maximum 

theta not to exceed maximum angle of attack for that aircraft. 
Keeping this weakness in mind, the results obtained from the math 
models are valid when determining the ’’best" escape maneuver. 
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Engine thrust is somewhat effected by airspeed and AOA The 
turbojet and turbofan-type engines can be significantly effected while 
the turboprop is effected to a much lesser degree. Only nominal 
values for thrust were used in the analysis. 

The effects of turbulence and rain on aircraft performance were 
not programmed Into the models. Studies by Wlngrove and Bach 
[Ref: 2] analytically determined that the rain had negligible effect 
upon DAL 191. The NTSB Report (National Transportation Safety 
Board) came to the same conclusion. What must be considered, is 
the effect of turbulence upon a pilot's ability to control the aircraft. 
The effect of turbulence on the ability of a pilot to execute a 
particular escape maneuver was kept in mind during the writing of 
the final conclusions. 

All the aircraft models are limited to one axis of freedom. The 
effects of coupled lateral and directional modes upon the longitudinal 
response were not modeled. From past analysis and personal 
aviation experience, this limitation should not significantly impact 
aircraft performance. Again, aircraft control may become difficult if 
an inherent mode is excited (such as the phugoid at low airspeed and 
high angle of attack). Escape maneuvers that may produce such 
dynamic modes were noted. 
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The effect of changing dynamic pressure on stick forces was 
evaluated. This is an important consideration from the pilot's 
perspective. A loss of 20-30 knots airspeed can induce a significant 
nose-down pitch force. 

C. DATA REDUCTION AND ANALYSIS 

Aircraft performance for three critical flight phases was evaluated 
when exposed to a DAL 191 type microburst windshear. Approach to 
landing and takeoff flight phases were scrutinized for the best escape 
maneuver that could be applied. The on-station flight phase was 
examined for the performance of a P-3 while encountering the 
microburst windshear and attempting to maintain altitude. 
Calculated data were recorded and converted to applicable units (eg., 
ft/s to knots). The data were then presented in a graphical format 
for analytical and subjective comparison. 

For all cases, tabulated data were converted to graphical form. 
Strip graphs for key dependent variables were produced for each 
aircraft performing a specific escape maneuver. They are listed in 
Appendices F and 6. Combination graphs of flight path, airspeed, and 
specific energy compare escape maneuvers for each aircraft in a 
given microburst encounter. The abscissa axis for all the graphs 
represents the distance from the microburst center in feet. A 
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three-step approach was carried out to deduce the "best" escape 
manuever for a particular encounter using these graphs. 

First, the inertial reference frame of flight path performance was 
analyzed. This involved using the combination flight path graph. A 
list is generated containing the highest to lowest altitude obtained at 
a specific point from the microburst center. For the approach to 
landing scenario, only the maneuvers that resulted in flight paths 
staying above the descent path (landing glideslope) were chosen. 
Altitude values were compared for each valid escape maneuver at 
x=5000 feet. For the takeoff scenario, only the escape maneuvers 
that did not lead to ground impact were listed. Here, altitude values 
were compared at x=4000 feet. 

The second step was to cross reference each selected maneuver's 
airspeed from the combination airspeed graph. This supplied an 
insight of the airmass reference performance. The airspeed values at 
the point of interest, and throughout the microburst encounter, 
allowed a subjective analysis of the validity of the results. Airspeeds 
that fell deep within the power-on stall region were scrutinized with 
the strip graphs comparing calculated AOA values. Any maneuver 
that resulted in a very low airspeed, high AOA condition was rejected 
from the list. 

The third step evolved using the specific energy combination 
graph for a final resolution. Specific energy incorporates both the 
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inertial and the airmass reference. The eligible maneuvers that had 
the highest specific energy value at the point of interest analytically 
are the "best" escape maneuvers for that aircraft configuration, 
given that type of microburst encounter. 

1. Approach to Landing Microburst Encounter Analysis 
Figure 12, Figure 13, and Figure 14 depict the results obtained 
for different escape maneuvers executed by a P-3 at 89,500 pounds 
gross weight. Table 6 compares the results observed at x=5000 feet. 
The constant airspeed and 12 unit AOA escape maneuvers were 
rejected initially due to resulting flight trajectories below the descent 
path (3° landing approach glideslope) as depicted in Figure 12. The 20 
unit AOA escape maneuver was subsequently rejected because of low 
airspeed, high AOA observed in Figure 13. Figure 14 shows that the 
5° theta escape maneuver results in the highest specific energy value 
of the remaining list. 



TABLE 6. RELATIVE VALUE RANKING FOR A P-3 AT 89,5000LBS. 



maneuver 


attitude 


airspeed 


specific energy 


20 unit AOA 


1 


rejected 


— 


const, altitude 


2 


5 


5 


15° theta 


3 


4 


4 


15 unit AOA 


4 


3 


3 


10° theta 


5 


2 


2 


5° theta 


6 


1 


1 



Figure 15, Figure 16, and Figure 17 depict the results obtained 
for different escape maneuvers executed by a P-3 at 114,000 pounds 
gross weight. Table 7 compares the results observed at x=5000 feet. 
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Figure 12. Flight path data comparison of a P3 at 
89,5001bs performing different escape maneuvers during 
an approach to landing microburst encounter. 
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Figure 13. Airspeed data comparison of a P3 at 89,5001bs 
performing different escape maneuvers during an 
approach to landing microburst encounter. 
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Figure 14. Specific energy data comparison of a P3 at 
89,5001bs performing different escape maneuvers during 
an approach to landing microburst encounter. 
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The constant airspeed, 12 unit AOA, and 5° theta escape maneuvers 
were rejected initially due resulting flight trajectories below the 
descent path as depicted in Figure 15. The 20 unit AOA and constant 
altitude escape maneuvers were subsequently rejected because of low 
airspeeds, high AOAs observed in Figure 16. Figure 17 shows that the 
15 unit AOA escape maneuver results In the highest specific energy 
value of the remaining list. 



TABLE 7. RELATIVE VALUE RANKING FOR A P-3 AT 114,000LBS. 



maneuver 


altltyide 


airsDeed 


- ? — — - 

SDecific enerav 


20 unit AOA 


1 


rejected 


— 


const, altitude 


2 


rejected 


— 


15° theta 


3 


3 


3 


10° theta 


4 


2 


2 


15 unit AOA 


5 


1 


1 



Figure 18, Figure 19, and Figure 20 depict the results obtained 
for different escape maneuvers executed by a 3-engine heavy airline 
transport. Table 8 compares the results observed at x=5000 feet. 
The constant airspeed, 5° theta, and 10° theta escape maneuvers 



TABLE 8. RELATIVE VALUE RANKING FOR A 3-ENGINE AIRLINE 



maneuver 


TRANSPORT. 

altitude 


airsDeed 


specific energy 


const, altitude 


1 


2 


2 


15° theta 


2 


1 


1 



were rejected initially due resulting flight trajectories below the 
descent path as depicted In Figure 18. Figure 20 shows that the 15° 
theta escape maneuver results in the highest specific energy value. 
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Figure 15. Flight path data comparison of a P3 at 
114,000113$ performing different escape maneuvers during 
an approach to landing microburst encounter. 
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Figure 16. Airspeed data comparison of a P3 at 114,0001bs 
performing different escape maneuvers during an 
approach to landing microburst encounter. 
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Figure 17. Specific energy data comparison of a P3 at 
114,0001bs performing different escape maneuvers during 
an approach to landing microburst encounter. 
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Figure 18. Altitude data comparison of a 3-engine heavy 
airline transport performing different escape maneuvers 
during an approach to landing microburst encounter. 
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Figure 19. Airspeed data comparison of a 3-engine heavy 
airline transport performing different escape maneuvers 
during an approach to landing microburst encounter. 
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Figure 20. Specific energy data comparison of a 3-engine 
heavy airline transport performing different escape 
maneuvers during an approach to landing microburst 
encounter. 
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Figure 21, Figure 22, and Figure 23 depict the results obtained 
for different escape maneuvers executed by a T-44. Table 9 
compares the results observed at x=5000 feet. The constant 
airspeed, 10 unit AOA, and 16 unit AOA escape maneuvers were 
rejected initially due to resulting flight trajectories below the descent 
path as depicted in Figure 21. The 25 unit AOA and constant altitude 
escape maneuvers was subsequently rejected because of low airspeeds 
and high AOAs observed in Figure 22. Figure 23 shows that the 5° 
theta escape maneuver results in the highest specific energy value of 
the remaining list. 



TABLE 9. RELATIVE VALUE RANKING FOR A T-44. 



maneuver 


altitude 


airspeed 


specific energy 


const, altitude 


-1 


rejected 


— 


15° theta 


2 


3 


3 


25 unit AOA 


3 


rejected 


— 


10° theta 


4 


2 


2 


5° theta 


5 


1 


1 



Figure 24, Figure 25, and Figure 26 depict the results obtained 
for different escape maneuvers performed by a P-3 at 89,500 pounds 
gross weight operating on 3 engines. Table 10 compares the results 
observed at x=5000 feet. The constant airspeed, 12 unit AOA, and 5° 
theta escape maneuvers were rejected initially due to resulting flight 
trajectories below the descent path as depicted in Figure 24. The 20 
unit AOA and constant altitude escape maneuvers were subsequently 
rejected because of low airspeeds and high AOAs observed in 
Figure 25. The 15° theta escape maneuver was also rejected for the 
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Figure 21. Altitude data comparison of a T-44 
performing different escape maneuvers during an 
approach to landing microburst encounter. 



64 




Figure 22. Airspeed data comparison of a 
performing different escape maneuvers during 
approach to landing microburst encounter. 



T-44 

an 
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Figure 23. Specific energy data comparison of a T-44 
performing different escape maneuvers during an 
approach to landing microburst encounter. 
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Figure 24. Altitude data comparison of a P-3 at 89,5001bs 
operating on 3 engines performing different escape 
maneuvers during an approach to landing microburst 
encounter. 
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Figure 25. Airspeed data comparison of a P-3 at 
89,5001bs operating on 3 engines performing different 
escape maneuvers during an approach to landing 
microburst encounter. 
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Figure 26. Specific energy data comparison of a P-3 at 
89,5001bs operating on 3 engines performing different 
escape maneuvers during an approach to landing 
microburst encounter. 
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same reasons. Although the airspeed is only within the flow 
separation region (steady state stall buffet), analysis of the AOA (see 
Figure F34) shows that critical alpha was sustained on the latter 
parts of the maneuver. Figure 26 shows that the 15 unit AOA escape 
maneuver results in the highest specific energy value of the 
remaining list. 



TABLE 10. RELATIVE VALUE RANKING FOR A P-3 AT 89,5000LBS 
OPERATING ON 3 ENGINES. 



maneuver 


altitude 


airsDeed 


SDecific enerRV 


const, altitude 


1 


rejected 


- 


15° theta 


2 


rejected 


— 


20 unit AOA 


3 


rejected 


— 


10° theta 


4 


2 


2 


15 unit AOA 


5 


1 


1 



2. Takeoff Microburst Encounter Analysis 

Figure 27, Figure 28, and Figure 29 depict the results obtained 
for different escape maneuvers performed by a P-3 at 90,0001bs gross 
weight. Altitude, airspeed, and specific energy values are compared 
at x=4000ft. All six escape maneuvers resulted in ground plane 
clearance as depicted in Figure 27. The 15 unit AOA escape maneuver 
was subsequently rejected because of low airspeed, high AOA as 
observed in Figure 28. Figure 29 shows that all escape maneuvers, 
except the 15 unit AOA escape maneuver, resulted in a grouped single 
value. This shows that all relevant escape maneuvers provide the 
same final specific energy. 
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Figure 27. Flight path data comparison of a P-3 at 
90,0001bs performing different escape maneuvers during 
a takeoff microburst encounter. 
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Figure 28. Airspeed comparison of a P-3 at 
9 0,000 lbs performing different escape maneuvers 
during a takeoff microburst encounter. 
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Figure 29. Specific energy comparison of a P-3 at 
90,0001bs performing different escape maneuvers 
during a takeoff microburst encounter. 
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Figure 30, Figure 31, and Figure 32 depict the results obtained 
for different escape maneuvers performed by a P-3 at 120,000 pounds 
gross weight. Value comparison is made at x=4000 feet. The 
constant airspeed, constant altitude, 12 unit AOA, 15 unit AOA, and 5° 
theta escape maneuvers were rejected initially due ground impact as 
depicted in Figure 30. The 15° theta escape maneuver was 
subsequently rejected because of low airspeed and high AOA as 
observed in Figure 31. The only viable, however marginal, 
performance observed was for the 10° theta escape maneuver. 
Figure 31 shows that the 10° theta maneuver resulted in an airspeed 
drop of 30 knots below V 2 at one point followed by a sustained 

airspeed of 20 knots below V 2 . 

Included with the 120,000 pound P-3 data is the resulting flight 
performance when V r was increased. A large gain in altitude and 

specific energy over other performance profiles is observed. 
Increasing the rotate speed by 18 knots coupled with flying a constant 
10° theta provides a 2 fold increase in altitude and a 21 % Increase in 
specific energy. This increase in specific energy is important when 
compared to the spread of specific energy values of different escape 
maneuvers. 



1 
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Figure 30. Flight path comparison of a P-3 at 
120,0001bs performing different escape maneuvers 
during a takeoff microburst encounter. 
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Figure 31. Airspeed comparison of a P-3 at 
120,0001bs performing different escape maneuvers 
during a takeoff microburst encounter. 
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Figure 32. Specific energy comparison of a P-3 at 
120,0001bs performing different escape maneuvers 
during a takeoff microburst encounter. 
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Figure 33, Figure 34, and Figure 35 depict the results obtained 
for different escape maneuvers performed by a P-3 at 135,0001bs 
gross weight. It can be graphically observed that a P-3 loaded to 
maximum takeoff weight is at the mercy of a strong microburst 
encountered at takeoff. Although the 15° theta maneuver misses the 
ground, the low airspeed and high AOA experienced disqualifies it as 
viable. 

Figure 36, Figure 37, and Figure 38 depict the results obtained 
for different escape maneuvers performed by a T-44. Altitude, 
airspeed, and specific energy values are compared at x=4000 feet. All 
five escape maneuvers resulted in ground plane clearance as depicted 
in Figure 36. All five escape maneuvers provided adequate margin of 
airspeed and AOA as depicted in Figure 37. Note that power-on stall 
speed for the T-44 Is less than 90 knots. Figure 38 shows that all 
escape maneuvers resulted in a grouped single value of specific 
energy. As seen above, all escape maneuvers for a given aircraft 
provide the same final specific energy during a takeoff microburst 
encounter. The only exception is when V r is increased. 

The results of one variant to the takeoff encounter are 
depicted in Figure 39, Figure 40, and Figure 41. Here, the microburst 
center is moved an added 800 feet from the liftoff point. This results 
in a delayed encounter with the severest part of the windshear. A 
P-3 at 120,00 pounds gross weight was used to perform the different 
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Figure 33. Flight path comparison of a P-3 at 
135,000 lbs performing different escape maneuvers 
during a takeoff microburst encounter. 



79 



£ 

£ £ 
C7* CT« 

< l > a* 
-o -o 

o in 




Figure 34. Airspeed comparison of a P-3 at 
13 5,000 lbs performing different escape maneuvers 
during a takeoff microburst encounter. 
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Figure 35. Specific energy comparison of a P-3 at 
13 5,0 00 lbs performing different escape maneuvers 
during a takeoff microburst encounter. 
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Figure 37. Airspeed comparison of a T-44 
performing different escape maneuvers during a 
takeoff microburst encounter. 
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Figure 38. Specific energy comparison of a T-44 
performing different escape maneuvers during a 
takeoff microburst encounter. 
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Figure 39. Flight path comparison of a P-3 at 
12 0,0 00 lbs performing different escape maneuvers 
during an after takeoff microburst encounter. 



85 

















< 

o 


T? 

a* 

o. 


< 

o 


a< 


a. 

"O 

Z? 






< 


Vf 


< 


O 




O' 










v— 




JZ 


jZ 


Kfl 




VJ 


k. 




c* 


c* 


c 




c 


> 




a. 


a* 


=1 


w 


=« 




VJ 


*o 


T? 


CM 


<z 

o 


in 


o 

c 


c 

o 


O 


in 


■» — 


o 


■» — 




o 


— 


■* — 




Figure 40. Airspeed comparison of a P-3 at 
120,0001bs performing different escape maneuvers 
during an after takeoff microburst encounter. 
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Figure 41. Specific energy comparison of a P-3 at 
120,0001bs performing different escape maneuvers 
during an after takeoff microburst encounter. 
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escape maneuvers including the increased V r technique. Table 11 



compares the results observed at x=4000 feet. The 12 unit AOA 
escape maneuver was rejected initially due to ground impact as 
depicted in Figure 39. The 15 unit AOA and 15° theta escape 
maneuvers were subsequently rejected because of low airspeeds and 
high AO As observed in Figure 40. Figure 41 shows that the close 
grouping of the specific energy values still exists for the 
non-incr eased V r maneuvers. Again, the increased V r resulted in a 

significant improvement in altitude and specific energy. 



TABLE 11. RELATIVE VALUE RANKING FOR A P-3 AT 120,000LBS 
WITH AN AFTER TAKEOFF ENCOUNTER. 



maneuver 


altitude 


airspeed 


specific energy 


increased V r 


l 


2 


1 


15° theta 


2 


rejected 


— - 


10° theta 


3 


4 


2 


const, altitude 


4 


3 


3 


15 unit AOA 


5 


rejected 


- - - 


const, airspeed 


6 


1 


4 



3. On-Station Encounter Analysis 

The analysis showed that a P-3 at 120,0001bs gross weight 
successfully navigated a strong microburst during loiter operations. 
Figure 42 graphically depicts the results. Altitude deviation was no 
greater than ± 20ft. Theta input never exceeded 10°. AOA remained 
below critical angle. Additional thrust was not applied until 40 knots 
of airspeed was lost in the encounter (a drop from 210 to 170 knots). 
Power was added at x=2700 feet and a time lapse of 23sec from the 
initial point of the encounter. 
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Figure 42. P-3 at 120,0001bs gross weight encountering a 
microburst windshear while on-station loiter. 
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D. OBSERVATIONS 



The results obtained from the analytical analysis for each 
encounter are integrated with other relative information to answer 
the three posed questions: 

1) What is the optimum microburst escape procedure given the 
type of encounter and the available flight performance 
information? 

2) Is the optimum escape procedure effected by change of gross 
weight or available thrust? 

3) What flight instrument indications, flight control feeling, and 
dynamic response would be expected during the optimum 
escape procedure? 

The answers to these questions are directed to the P-3. However, 
comparison is made to the T-44 and 3-engine heavy airline transport 
where appropriate. Other important issues must also be addressed 
for comparison of results to other wlndshear studies. These issues 
include the impact of wing loading (W/S), thrust to weight (T/W), and 
early liftoff speeds for the P-3. 

1. Available Aircraft Flight Instruments 

The available flight information and type of presentation are 
extremely influential on the choice of the optimum escape procedure. 
Both the P-3 and T-44 have "conventional" flight instruments [Ref: 8 
and Ref: 9], Conventional implies gyroscopic-stabilized attitude and 
heading indication in combination with pitot-static airspeed, altitude, 
and vertical speed indication (VSI). Angle of Attack (AOA) indication 
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is installed on the P-3 and T-44. However, the AOA is heavily 
damped and is primarily designed for steady-state flight conditions. 
Standard radio navigation equipment is installed on both aircraft 
encompassing an Integrated Landing System (ILS), VHF Omni Radio 
(VOR), and Tactical Navigation (TACAN). The P-3 in addition has 
Inertial Navigation Systems (INS), but no information is presented at 
the flight station. Both aircraft incorporate autopilots and flight 
directors. The P-3 autopilot is controlled by Control Wheel Steering 
(CWS) input only. The T-44 autopilot is controlled through thumb 
wheels and can couple to the flight director. Flight directors on both 
aircraft can provide navigational steering. However, the primary 
design of the P-3 flight director is tactical while the primary design 
of the T-44 flight director is navigational. In summary, both the P-3 
and T-44 have Attitude Indicators (Al), damped AOA indication, and 
pitot-static instruments for reference while executing an escape 
maneuver. 

2. Optimal Escape Procedure 

Flying a constant theta with reference to the Al seems to be 
the optimal escape procedure during a microburst encounter. This 
conclusion is supported analytically and qualitatively for the P-3 and 
T-44. Analytically, superior performance was obtained in a few 
circumstances from constant AOA maneuvers with an alpha value 
between approach and stall. However, qualitatively the AOA 
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indication in both aircraft docs not lend itself as a viable reference 



for two reasons. First, the AOA indication may be significantly in 
error because of instrument lag in turbulent conditions. Second, P-3 
and T-44 pilots seldom use AOA indication during the approach phase. 
The transition from attitude/airspeed to attitude/AOA reference 
during an intense situation as a microburst encounter, is more than 
can be expected from a pilotage standpoint. The only viable 
alternative to flying constant AOA was a constant theta maneuver. 
Note that quantitatively in all circumstances, a constant theta 
maneuver was superior to or closely matched with any other 
attempted escape maneuver. 

Theta values ranging from 5“ to 15° were identified as optimal 
for an approach to landing encounter. The theta value of 10°, if not 
optimal, provided suitable recovery for the P-3 in the configurations 
analyzed. The light, the heavy, and the 3-engine P-3 successfully 
navigated the microburst encounter utilizing the 10° theta escape. 
This result is important from the standpoint that 10° is easy to 
remember and easy to read on the AI. Also, 10° theta worked well 
for the T-44 in the approach configuration (approach flap setting). A 
theta of 15° was optimum for the T-44 in the takeoff configuration 
(flaps up). This is important because future P-3 pilots receive their 
first training in T44s. 
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3. The Effect Of Specific Energy On Survival Probability 

Escaping from a takeoff microburst encounter has a lower 
probability for success compared to the approach to landing 
encounter. The type of escape maneuver selected in a takeoff 
encounter has a lesser effect on the outcome. This result is caused 
by the excess airspeed and additional altitude available during the 
approach phase. (Note that the P-3 community flies an unusually 
fast approach speed compared to their civilian counterparts. The P-3 
approach V re f airspeed is determined from 1.35V s +5 knots [Ref. 20], 

which equates to 1.4V S for a gross weight of 114,000 pounds. The 
industry standard for V rc f is 1.3V s + an added factor as dictated by 

the type of aircraft.) Specific energy can be used for an interesting 
comparison between an approach to landing and a takeoff microburst 
encounter. A P-3 weighing 89,000 pounds has specific energy values 
of 1230 feet and 1040 feet at x=500 feet and x=4000 feet respectively 
during an approach to landing encounter (reference Figure 14). 
During a takeoff encounter, a P-3 at 90,000 pounds has a specific 
energy values of 810 feet and 980 feet at the same x distances from 
the microburst center (reference Figure 29). Both executed a 10° 
theta escape maneuver and both successfully navigated the 
microburst. Note that the approach to landing phase lost energy 
while the takeoff phase gained energy. 
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An Important aspect of specific energy comparison is seen in 
the increased rotate speed data. Figure 41 graphically depicts a 
120,000-pound P-3 starting with an additional specific energy of 200 
feet compared to the normal rotate speed profiles. The exit specific 
energy was 175 feet greater than the comparable 10° theta 
maneuver. The 175 feet were translated almost completely to 
altitude (reference Figure 30 and Figure 31). Data support that an 
increased rotate speed is far more beneficial than any particular 
escape maneuver during a takeoff microburst encounter. 

4. Stick Force vs. Off-Trim Airspeed For The P-3 

The expected change in flight control "feel" for the P-3 is light 
compared to large aircraft. Flight tests show that stick force versus 
change from trim airspeed results in a shallow gradient for this size 
aircraft [Ref. 18:p. 37]. In the takeoff/approach flap configuration (18°), 
with an aft C.6. of 29 % M.A.C., a stick force of an 11-pound pull is 
needed with a 30-knot decrease from trim airspeed. Stick force 
increases to a 16 pound pull with the most forward C.G. of 16S M.A.C. 
Flight test data [Ref. 18] also show that the stick force to trim 
airspeed gradient is not appreciably effected by landing gear position 
or gross weight. The same flap setting (takeoff/approach) is used for 
landing approach and takeoff. Thus, the elevator "feer 1 is expected to 
be the same for either approach to landing or takeoff microburst 
encounters. 
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The P-3 Flight Simulator (Device 2F87F) was used to confirm 
the flight test data relating to off-speed to stick force and landing 
gear drag. The simulator allows one or more flight parameters to be 
frozen. This feature is used by freezing the altitude to 500 feet and 
trimming the aircraft with maximum power at reference airspeed. 
The airspeed Is then reset and frozen to the low airspeed value 
expected during a microburst encounter. The control force required 
to maintain a pitch attitude was then directly measured and printed 
out. The P-3 flight simulator data (Appendix H) showed that the 
elevator pull needed at Vtrim minus 30kts decreased by 5 pounds 

when full power was applied. Level acceleration maneuvers were 
performed with the landing gear extended, retracted, and during 
retraction. Simulator data showed that the landing gear retraction 
cycle does not increase drag. 

5. Weight, Wing Loading, and Thrust to Weight Effects 

It is obvious from the analysis that a heavy P-3 is effected by 
microburst windshear to a higher degree than a light P-3. The T-44 
was relatively unaffected by the windshear. The question arises 
whether weight, wing loading, or thrust to weight has the largest 
effect on microburst survival. Table 12 compares the weight (Wt), 
wing loading (W/S), thrust to weight (T/W), and specific energy values 
for three different aircraft during an approach to landing microburst 
encounter. The specific energy loss (E s loss) column shows the energy 
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(in feet) lost during the encounter. The specific energy spread (E e 

spread) column shows the spread (in feet) of the different escape 
techniques measured after the encounter. Apparently, the higher 
the thrust to weight ratio, the less energy is lost. Light wing loading 
equates to an individual escape procedure having less effect upon the 
result. Weight, as a separate parameter, is irrelevant. In other 
words, increased thrust to weight decreases the effect of the 
windshear upon the aircraft and light wing loading allows increased 
efficiency of energy transformation. 



TABLE 12. WEIGHT, WING LOADING, AND THRUST TO WEIGHT 
EFFECTS ON P s VALUES. 



aircraft 


wt (lb?) 


W/$ 


IZW 


Ec loss 


Ec spread 


P-3 


89,500 


68.85 


0.37 


380 


145 


P-3 


114,000 


88.08 


0.29 


455 


155 


T-44 


8,280 


39.43 


0.37 


285 


50 


3 -eng. 
hvy. 


362,000 


79.07 


0.35 


345 


100 



6. Early Liftoff Speed For The P-3 

The analyses did not look at the effect of a microburst 
windshear encountered during a takeoff roll. The FAA Windshear 
Training Aid [Ref. 14] points out that if windshear is encountered past 
V t , the aircraft should be committed to flight no later than with 

2000 feet of runway remaining. Although the airspeed may fall and 
stay below V r , most aircraft can go airborne. In this regard, a 

theoretical liftoff speed was calculated for the P-3. Equations 38, 40, 
and 42 were used in conjunction to solve for V. The Cia data were 
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extrapolated from flight test data [Ref. 18]. Table 13 compares the 
recommended V r speed [Ref. 20] to the theoretical liftoff speed for a 

rotate pitch attitude of 10°. Ground effect was not considered in the 
calculations. 



TABLE 13. THEORETICAL LIFTOFF SPEEDS 



weight (lbs) 

80,000 

90.000 

100.000 

110,000 

120,000 

130.000 

140.000 



V. (kts) 


theoretical liftoff (kts) 


115 


90 


115 


95 


115 


101 


117 


105 


123 


110 


129 


115 


134 


119 



8. Summary 

The original three questions are answered for the P-3: 

1) The optimum approach to landing microburst escape 
procedure is to set and maintain maximum power while 
simultaneously setting a 10° pitch attitude. For the takeoff 
microburst encounter, increase rotate speed to 140 knots 
then pitch to and maintain 10°. 

2) The optimum escape procedures remain the same for all 
gross weights. The approach to landing encounter escape 
procedure is the same for four or three operating engines. 

3) The airspeed will rapidly decay and remain abnormally low 
during a microburst penetration. An elevator force of 5 to 
10 pounds can be expected to maintain a 10° pitch attitude 
with full power. Flight tests show that the P-3 has no 
unusual short period, phugoid, nor cross-couple dynamics to 
contend with. 
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The answers for the 1-44 should be similar to the P-3. The 10° pitch 
attitude was optimal for the approach to landing encounter. A 
pitch attitude of 15° is optimal for takeoff, the difference being that 
the flaps are up in the takeoff configuration. Not enough 
performance data were available to predict the expected stick forces. 
However, personal flight experience has shown that the 1-44 exhibits 
longitudinal forces and responses very similar to the P-3. 
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IV. CONCLUSIONS 



The aim of this thesis was to analytically produce microburst 
escape flight procedures. Study was directed at two turboprop 
aircraft the Navy routinely operates. The P-3 and T-44 were closely 
scrutinized for the optimal escape procedure if confronted with a 
microburst penetration. The microburst modeled was patterned 
after that encountered by a Delta Airlines L-1011 at Dallas/Ft. Worth 
on August 2, 1985. The conclusions drawn are presented within the 
format of the FAA Wlndshear Training Aid [Ref. 14]. Viable flight 
procedures and associated precautions are presented for both aircraft. 
The precautions should be adhered to if an encounter with a 
microburst windshear is a possible expectation. NATOPS change 
recommendations reflecting these conclusions will be forwarded to the 
appropriate aircraft model managers. 

A. MICROBURST ESCAPE PROCEDURES FOR THE P-3 

A target pitch attitude of 10° is optimal for all weights and with 
four or three engines operating. This target pitch attitude is used for 
approach to landing and takeoff microburst encounters with the flaps 
in the Takeoff /Approach position. The airspeed at which the P-3 lifts 
off is the critical factor in the takeoff scenario. 
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1. Approach To Landing Microburst Encounter 

Precautions: 

1) Avoid thunderstorm conditions. Delay the approach if 
possible. 

2) Use Takeoff /Approach flap setting. Do not extend flaps 
to Land position until the runway is made. 

3) Use a precision approach procedure. 

4) The approach speed, V re f, should be equal to 1.35V3 + 5 
knots. Additional airspeed is not normally warranted. 
This V re f will provide a comfortable margin above stall 
speed. Note that additional airspeed increases landing 
distance (5 knots faster, 102 farther). 

5) Determine expected descent rate. 

6) Consider using the autopilot (ASW-31) with the altitude 
hold switch in the off position. The autopilot will 
provide wing rock damping and pitch attitude 
augmentation In the face of turbulence. It will also 
provide pitch attitude hold if a microburst escape 
maneuver Is warranted. Insure autopilot disconnect 
passing 200 feet AGL. 

7) Attain a stabilized airspeed approach before passing 
1000 feet AGL. Minimize power lever movement 
beyond this point. Maintain glideslope with pitch. The 
airspeed Indicator will serve as a windshear indicator. 

8) Strong consideration should be given to executing the 
microburst escape procedure if one of the indications is 
observed: 

a) a rapid and sustained airspeed loss of 20 knots below 
Vref I 

b) a descent rate 500 feet per minute greater than the 
predetermined value; 

c) greater than 1 dot low on the ILS glideslope associated 
with airspeed 10 knots below V re f; or 

d) a "well below glideslope" call on a PAR associated 
with airspeed 10 knots below V ref . 
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Escape procedure: 

1) Apply maximum power (power levers to the stops). 

2) Set and maintain a pitch attitude of 10° on the AI. Do 
not attempt to recover airspeed. 

3) Once a positive rate of climb is established, select 
landing gear up. 

4) Do not raise the flaps until the airspeed has increased 
above 140 knots indicating exit from the windshear. 

2. Takeoff Microburst Encounter 

Precautions: 

1) Delay the takeoff if able. 

2) Select the longest suitable runway. 

3) Perform takeoff planning for adverse conditions as 
prescribed by the NATOPS Flight Manual. 

4) Increase rotate airspeed to 140 knots, or the airspeed 
that is attained with 2000 feet of runway remaining. 
This is determined by using the Four Engine 
Acceleration Chart In NATOPS and following this 
procedure: 

a) Subtract two thousand feet from the available 
runway length. 

b) Enter the chart with the adjusted runway length, 
pressure altitude, temperature and gross weight. Exit 
the chart with the corresponding airspeed value. Use 
140 knots if the chart value is higher. 

c) Corrections for runway slope, winds, or standing 
water is not required for THIS prediction. 

5) Thoroughly brief the takeoff procedure, voice calls, and 
windshear indications among the flight crew. 
Penetration of a windshear during takeoff will be 
indicated by a loss of airspeed or no airspeed 
acceleration. 
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Procedure: 



1) Abort the takeoff if windshear is indicated before 
reaching refusal airspeed. Reaching refusal speed 
ground roll distance with airspeed significantly below 
OR above refusal speed indicates windshear. 

2) Continue the takeoff if windshear is experienced after 
refusal speed. Rotate the nose when: 

a) predicted increased rotate speed, usually HO knots, is 
attained, or 

b) 2000 feet of runway remain. 

3) Rotate to a pitch attitude of 10°. DO NOT delay rotating 
the nose because of low airspeed. 

4) Increase pitch attitude toward 15° if ground impact is 
imminent. 

5) Raise the landing gear when above 100 feet AGL with a 
positive rate of climb. 

6) Do not raise the flaps until normal climb airspeed is 
regained, indicating clear of the microburst. 

3. On-Station Loiter Microburst Encounter 

No immediate effect on the aircraft performance was 
observed. However, certain precautions should be taken and 
procedures followed if flying through convective activity at low 
altitude: 

1) Use the autopilot altitude hold in dual axis mode. 
Note that if the autopilot is not used, the pilot can 
expect abnormal elevator control force changes. 

2) Minimize power lever movement. Do not pull off 
power for high airspeed indication. 
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3) A penetration of a microburst will be Indicated by a 
significant loss of airspeed. The first and foremost 
reaction should be the addition of full power. Monitor 
autopilot input and altitude hold. 

4) If a climb is deemed warranted, set 10° pitch attitude 
using the CWS function of the autopilot. 



B. MICROBURST ESCAPE PROCEDURES FOR THE T-44 

Target pitch attitude escape procedures are Just as effective for a 
light turboprop as the T-44. In its case, the target pitch values are 
different between takeoff and landing. This change is owed to 
different flap settings. Many of the same precautions and 
procedures are the same between the T-44 and P-3. 

1. Approach To Landing Microburst Encounter 
Precautions: 

1) Avoid thunderstorm conditions. Delay the approach if 
possible. 

2) Use Approach flap setting. Do not extend flaps to full 
down position unless required and until the runway is 
made. 

3) Use a precision approach procedure. 

4) The approach speed, V re f, should be equal to 120 knots. 
Additional airspeed is not warranted. This V re f will 
provide a comfortable margin above stall speed. Note 
that additional airspeed increases landing distance (5 
knots faster, 10£ farther). 

5) Determine expected descent rate. 
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6) Consider using the autopilot in the coupled mode if 
executing an ILS. However, close monitoring of the 
autopilot is required because of abnormal elevator pitch 
forces expected in a windshear. Insure autopilot 
disconnect passing 200 feet AGL. 

7) Attain a stabilized airspeed approach before passing 
1000 feet AGL. Minimize power lever movement 
beyond this point. Maintain glideslope with pitch. The 
airspeed indicator will serve as a windshear indicator. 

8) Strong consideration should be given to executing the 
microburst escape procedure if one of the indications is 
observed: 

a) a rapid and sustained airspeed loss of 15 knots below 

Vre {', 

b) a descent rate 500 feet per minute greater than the 
predetermined value; 

c) greater than 1 dot low on the ILS glideslope associated 
with airspeed 10 knots below V re f; or 

d) a "well below glideslope" call on a PAR associated 
with airspeed 10 knots below V re f. 

Escape procedure: 

1) Apply maximum power (power levers to the stops). 

2) Set and maintain a pitch attitude of 10° on the AI. DO 
NOT attempt to recover airspeed. 

3) Once a positive rate of climb is established, select 
landing gear up. 

4) Do not raise the flaps until the airspeed has increased 
above 120 knots, indicating exit from the windshear. 
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2. Takeoff Microburst Encounter 



Precautions: 

1) Delay the takeoff If able. 

2) Select the longest suitable runway. 

3) Perform takeoff planning for adverse conditions as 
prescribed by the NATOPS Flight Manual. 

4) Increase rotate airspeed to 120 knots, or the airspeed 
that is attained with 2000 feet of runway remaining. 

5) Thoroughly brief the takeoff procedure, voice calls, and 
windshear indications among the flight crew. 
Penetration of a windshear during takeoff will be 
indicated by a loss of airspeed or no airspeed 
acceleration. 

Procedure: 

1) Abort the' takeoff if windshear is indicated before 
reaching 91 knots airspeed. 

2) Continue the takeoff if windshear is experienced after 91 
knots airspeed. Rotate the nose when: 

a) 120 knots; or 

b) 2000 feet of runway remain. 

3) Rotate to a pitch attitude of 15°. DO NOT delay rotating 
the nose because of low airspeed. 

4) Increase pitch attitude toward 20° if ground impact is 
imminent. 

5) Raise the landing gear when above 100 feet AGL with a 
positive rate of climb. 
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V. RECOMMENDATIONS 



Recommendations are broken into two sections. The first set are 
for analysis refinement and further research. The second set are 
directed to the P-3 operators. 



A. ANALYSIS REFINEMENT AND CONTINUED RESEARCH 

There were some weak areas of this analysis that were pointed 
out earlier. Although valid results were obtained, refinement of the 
windshear model and the aircraft equations of motion would further 
explore the edges of the operating envelopes. The following 
recommendations are provided for the microburst windshear model: 

1) Vortex wind field models seem to provide the closest 
representation of a microburst windshear, short of applying 
Navier-Stokes principles. Schultz's multiple vortex model 
combined with source flows should be combined and fitted to 
other previously recorded microbursts. 

2) A more sophisticated parameter estimation scheme should be 
applied to enhance and accelerate parameter fits of recorded 
data to the vortex/source flow model. 

3) Realistic boundary conditions to the windshear model should be 
defined. This would help explore the entry/exit aircraft 
response. 

The following recommendations are for the aircraft performance 
models: 
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1) The lift and drag equations should be improved from second to 
fourth order equations. This change will enhance the 
representation of the aircraft performance at low airspeed/high 
angle of attack. 

2) An available thrust algorithm should be integrated into the 
aircraft performance model. Comparison of turbojet, turbofan, 
turboprop, and reciprocating engine propulsion may lead to 
significant effects on performance through windshear. 

3) A state-space matrix of aircraft dynamics should be integrated 
into the performance algorithm. Some aircraft may exhibit 
unexpected dynamics associated with microburst windshear 
and escape maneuvers. 



B. P-3 OPERATIONS AND FLIGHT CREW TRAINING 

An engineering analysis provides the performance specifications 
for penetrating a microburst and escaping. However, certain issues 
must be addressed by the P-3 community so as to successfully 
employ the results: 

1) Flight crew coordination training is essential for successful 
employment of microburst escape procedures. All members of 
the flight deck must understand the teamwork required to 
execute these procedures. Therefore, comprehensive windshear 
training should be developed. The FAA Windshear Training Aid 
and this thesis provides an Initial foundation. 

2) A windshear algorithm should be developed for the P-3 Flight 
Simulators. There Is no effective inflight means to expose 
flight crews to microburst effects. Nor can escape procedure be 
effectively practiced. The only proven instructional means are 
with a flight simulator. 

3) A brief dissertation on microburst and gust front windshear 
should be provided in the flight station NATOPS manual. 
NATOPS change recommendations reflecting the conclusions of 
this study should be adopted. 
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APPENDIX A 



WINDSHEAR MODEL ALGORITHMS 



The following are MatLab® functions which model the winds 
experienced at DFW on August 2, 1985 by two separate aircraft. 
AA539 penetrated the microburst at 2500 feet. DAL191 penetrated 
the same microburst while attempting a landing. 

A. AA539 WINDSHEAR MODEL ALGORITHM 



function [ux,uz]*inicroburst(ze,x) 



X INITIAL UALUES AND INPUTS 
X Earth position 
Xze ■ altitude (ft) 

Xx horizontal position (ft) 

X Uortex Aing Dimensions 
AL-8503.3; 
r i L-2001 . 1 ; 
kL-131968.8; 

roiL-riL/0.371 ; 
yL=3350.4; 

LUoralt-3100.6; 

As-1701 .7; 
ris-323.9; 
ks-57201.9; 

roes«ris/0.371 j 

ys»830.9; 

sUoralt*2333.6; 



Xradius of ring filament (ft) 

Xcore radius (ft) 

Xvortex strength of the ring 
Xf i lament(ft A 2/sec) 

Xrelated core radius 
Xlarge ring y lateral position (ft) 
X large vortex altitude A6L (ft) 
Xradius of ring filament (ft) 

Xcore radius (ft) 

Xvortex strength of the ring 
Xfi lament (ft A 2/sec) 

Xrelated core radius 

Xsmall ring y lateral position (ft) 

Xsmall vortex altitude AGL (ft) 



no 



X PRItIRRV LRR6E RING CRLCULRTIONS 
z*ze-LUoralt; 
if x=*0, 
if yL=*0, 
uxL e O; 

uzL=-kL*RL A 2/(2*(RL A 2+z A 2) A 1.5); 

zerock s O; 

end 

else 

zerock=999; 
r»sqrt(x A 2+yL A 2); 
rl-sqrt (z A 2+(r-RL) A 2) ; 



r2=8qrt (z A 2+(r+RL) A 2) ; 



Xradial distance from the z axis 
Xclosest distance to the point of 
Xinterest from the ring filament. 
Xfarthest distance to the point of 
Xinterest from the ring filament. 
Xscal ing term 



Ia»bda*(r2-r1 )/(r2+r1 ); 
nl=(r+RL)/r2; 
n2=(r-RL)/r1; 
sig1=z/r2+z/r1 ; 
sig2=z/r2-z/r1 ; 

t oi*0 . 75*sqrt (1-1 ambda A 2)+0 . 25; 
delta1*(kL*lambda)/(pi*tov); 

de 1 1 a2» (0 . 2955*kl* I ambda A 3)/ (p i *t o« A 2*sqrt (1-1 ambda A 2) ) ; 
dS i dr=0 . 788*de I tal *(n1 -n2)-0 . 394* I ambda*de I t al *(n1 +n2) 

+de I ta2*( (nl -n2)- I a»bda*(n1 +n2) ) ; 

Xderiuatiue of the stream function vrt r 
dS i dz e de 1 1 al *(0 . 788*s i g2-0 . 394* I ambda*s i g 1 ) 

+delta2*(sig2-lambda*sig1 ); 

Xderiuatiue of the stream function vrt z 
dmpL=(1-exp(-r1 A 2/(0.1*rovL A 2))); 

Xuelocity damping factor 

uxL«(x/r)*(1/r)*dSidz; Xhorizontal uelocity in x direction 

uzL*-(1/r)*dSidr; Xuertical uelocity 

end Xif 
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X PRinRRV SI1RLL RING CfiLCULRT I OHS 
z=ze-sUoralt; 
i f x==0, 
if ys==0, 
oxs=0; 



uzs=-ks*Rs A 2/(2*(Rs / '2+z A 2) A 1 .5) ; 
end 
else 

r=sqrt((x-50) A 2+ys A 2) ; 
r1*sqrt(z A 2+(r-Rs) A 2) ; 



r2=sqrt(z A 2+(r+Rs) A 2) ; 



Xradial distance from the z axis 
Xclosest distance to the point of 
X interest from the ring filament. 
Xfarthest distance to the point of 
X interest from the ring filament. 
Xscal ing term 



Iambda*(r2-r1 )/(r2+r1 ) ; 
n1=(r+Rs)/r2; 
n2*(r-Rs)/r1 ; 
sig1*z/r2+z/r1 ; 
sig2 s z/r2-z/r1 ; 

t om»0 . 75*sqrt (1-1 ambda A 2) +0 . 25 ; 
delta1«(ks*lambda)/(pi*tom) j 

delta2«(0.2955*ks*lambda A 3)/(pi*toi A 2*sqrt(1-lambda A 2))j 
dS i dr=0 . 788*de 1 1 a 1 * (n 1 -n2)-0 . 394* I ambda*de 1 1 al 
*Cn1+n2)+delta2*((n1-n2)-lambda*(n1+n2)); 



Xderioative of the stream function »rt r 
dS i dz=de 1 1 al *(0 . 788*s i g2-0 . 394* I ambda*s i gl ) 

+delta2*(sig2-lambda*sig1 ); 

Xderiuatiue of the stream function «rt z 



if (r1 A 2/(0. 1*roms A 2))>444, 
dmps=1 ; 

else Xuelocity damping factor 

dmps*1 -exp(-r 1 A 2/(0 . 1 *roms A 2) ) ; 
end 

oxs«(x/r)*(1/r)*dSidz; Xhorizontal uelocity x direction 

uzs«-(1/r)*dSidr; Xoertical uelocity 

end Xif 
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X lOflGE LARGE RING CRLCULRTIONS 
x»x; 

z=-ze-LUoralt; 

I f x==0, 
if yL==0, 
vxLi*0; 

vzL i B -kL*RL A 2/(2*(RL A 2+z A 2) A 1 . 5) ; 



n1=(r+RL)/r2; 
n2=(r-RL)/r1 ; 
sig1«z/r2+z/r1 ; 
sig2*z/r2-z/r1 ; 

to»*0 . 75*sqrt (1-1 aebda A 2)+0 . 25; 
delta1 s (kL*laabda)/(pi*toi); 

de I ta2»(0 . 2955*kL* I aebda A 3)/(p i *t o« A 2*sqrt (1-1 aebda A 2 ) ) ; 
dS i dr B 0 . 788*de 1 1 a 1 * ( n 1 -n2 ) -0 . 391* I a»bda*de 1 1 a 1 
*(n1 +n2)+de I t a2*( (nl -n2)- I atbda*(n1 +n2) ) ; 



dS i dz-de 1 1 a 1 * (0 . 788*s i g2-0 . 391* I aibda*s i gl ) 

+de 1 1 a2* ( s i g2- 1 aebda*s i g 1 ) ; 

Xderivative of the streae function art z 
dtpL i -( 1 -exp(-r1 A 2/(0. 1 *ro«L A 2 ) ) ) ; 

Xvelocity daaping factor 

vxLi B (x/r)*(1/r)*dSidz; Xhorizontal velocity in x direction 

vzU B -(1/r)*dSldr ; Xvertical velocity 



end 

else 

r*sqrt(x A 2+yL A 2); 

r1*sqrt(z A 2+(r-RL) A 2); 



Xradial distance from the z axis 
Xclosest distance to the point of 
X interest fron the ring filaeent. 
Xfarthest distance to the point of 
X interest fro« the ring filaeent. 
Xscal ing tern 



r2-sqrt(z A 2+(r+RL) A 2); 



Ia»bda=(r2-r1 )/(r2+r1 ); 



Xderi native of the streai function «rt r 



end Xif 
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{ inRGE SHALL AIHG CALCULATIONS 
z=-ze-sUoralt; 
if x==0, 
if ys==0, 
vxsi s O; 

vzs i =-ks*As A 2/(2*(Rs A 2+z A 2) A 1 . 5) ; 
end 
else 

r*sqrt ( (x+50) A 2+ys A 2) ; 
r1=sqrt(z A 2+(r-Rs) A 2); 



r2=sqrt(z A 2+(r+Rs) A 2); 



{radial distance from the z axis 
{closest distance to the point of 
{interest from the ring filament, 
{farthest distance to the point of 
{interest from the ring filament, 
{seal ing term 



Iambda=(r2-r1 )/(r2+r1 ); 
n1=(r+Rs)/r2; 
n2*(r-Rs)/r1 ; 
sig1*z/r2+z/r1 ; 
sig2*z/r2-z/r1 ; 

to»=0 . 75*sqrt (1-1 ambda A 2)+0 . 25; 
delta1 e (ks*lambda)/(pi*to»); 

de I ta2=(0 . 2955*ks* I ambda A 3)/(p i *tom A 2*sqrt (1-1 ambda A 2) ) ; 
dS i dr«0 . 788*de 1 1 a 1 * ( n 1 -n2) -0 . 391* I ambda*de 1 1 a 1 
*(n1+n2)+delta2*((n1-n2)-lambda*(n1+n2)); 

{derivative of the stream function mrt r 
dS i dz=de I t al *(0 . 788*s i g2-0 . 391* I ambda*s i gl ) 

+delta2*(sig2-lambda*sig1 ); 

{derivative of the stream function mrt z 



if (r1 A 2/(0.1*roms A 2))>111, 
dmpsM ; 

else {velocity damping factor 

dmps i ■ 1 -exp(-r 1 A 2/(0 . 1 *roms A 2) ) ; 
end 



vxs i =(x/r)*( 1 /r)*dSidz; 
vzsi*-(1/r)*dSidr; 
end {if 



{horizontal velocity x direction 

{vertical velocity 
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x sumption 

if zerock==999, 
dmp=dsipL*draps*d»pL i *d»ps i ; 
ox=-(dmp*(uxL+uxs+uxL i +uxs i ) ) ; 
uz*dmp*(vzL+uzs-ozLi-uzsi ); 
else 

ox E -(uxL+uxs+uxL i +uxs i ) ; 
uz^uzL+uzs-uzLi-uzsi ; 
end 



B. DAL191 WINDSHEAR MODEL ALGORITHM 
function [ux > uz]*»icnoburst (ze,x) 



X INITIRL URLUES AND INPUTS 
X Earth position 
Xze s altitude (ft) 

Xx horizontal position (ft) 

X Uortex Ring Oi sens ions 
RL-7000; 
riL=2004. 1 ; 
kL«431968.8*l j 

romL-rlL/0.371; 

yL*-300; 

LUoralt*3400; 

Lxdisp=2500; 

Rs=1300; 

ris = 323.9; 

ks=57204.9*2.3; 

noms=nis/0.371 ; 
ys=i ; 

sUoralt=800; 

sxdisp E 300; 



Xradi'us of ring filanent (ft) 

Xcore radius (ft) 

Xuortex strength of the ring 
Xfi latent (ft A 2/sec) 

Xrelated core radius 
X large ring y lateral position (ft) 
Xlarge uortex altitude RGL (ft) 

X large ring displacement in x (ft) 

Xradius of ring filament (ft) 

Xcore radius (ft) 

Xuortex strength of the ring 
Xfi lament (ft A 2/sec) 

Xrelated core radius 
Xsmall ring y lateral position (ft) 
Xsmall uortex altitude RGL (ft) 
Xsmall ring displacement in x (ft) 
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X PRinflRV LARGE RING CRLCULRTIONS 
z=ze-LUonalt ; 
i f x bb O, 
if yL-O, 
uxL=0; 

uzL=-kL*RL A 2/(2*(RL A 2+z A 2) A 1 .5); 



zerock x O; 

end 

else 

zerock=999; 
p»sqpt(x A 2+yL A 2)j 
pi =sqpt (z A 2+(p-RL) A 2) ; 

p2*sqpt(z A 2+(p+RL) A 2); 



Xpadial distance fpom the z axis 
Xclosest distance to the point of 
Xintepest fpo» the ping filaaent. 
Xfapthest distance to the point of 
Xintepest fpoa the ping filaaent. 
Xscal ing tepa 



Iaabda*(p2-pl )/(p2+p1 ); 
n1*(p+RL)/p2; 
n2*(p-RL)/p1 ; 
sig1»z/p2+z/p1 ; 
sig2 x z/p2-z/p1 ; 

t oa*0 . 75*sqpt (1-1 aabda A 2)+0 . 25 ; 
delta1 B (kL*laabda)/(pi*to«); 

de I ta2 B (0 . 2955*kL* I aebda A 3)/(p i *toa A 2< t sqpt (1-1 aabda A 2) ) ; 
dSidp«0.788*delta1*(n1-n2)-0.391*laabda*delta1*(n1+n2) 
♦delta2*((n1-n2)-laabda*(n1+n2)); 

Xdepiuatiue of the stpeaa function apt p 
dS i dz*de 1 1 a 1 *(0 . 788*s i g2-0 . 391* I aabda*s i g 1 ) 

+delta2*(sig2-laabda*sig1 ) ; 

Xdepiuatiue of the stpeaa function apt z 
dapL-( 1 -exp (-p 1 A 2/(0 . 1 *poaL A 2) ) ) ; 

Xuelocity daaping fact op 

uxL*(x/p)*(1/p)*dSidzj Xhopizontal uelocity in x dipection 

uzL*-(1/p)*dSidp; Xueptical uelocity 

end Xif 
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X PRIMfiRV SMALL RING CRLCULRTIONS 
z a ze-sUoralt ; 
i f x**0, 
if y8"0, 
vxs=0; 

vzs*-ks*Rs A 2/(2*(Rs A 2+z A 2) A 1 .5) ; 
end 
else 

r=sqrt ( (x-50) A 2+ys A 2) ; 
r1*sqrt(z A 2+(r-Rs) A 2); 



r2*sqrt (z A 2+(r+Rs) A 2) j 



Xradial distance from the z axis 
Xclosest distance to the point of 
^interest from the ring filament, 
^farthest distance to the point of 
^interest from the ring filament. 
Xscal ing term 



Iambda*(r2-r1 )/(r2+r1 ); 
n1*(r+Rs)/r2; 
n2*(r-Rs)/r1 ; 
sig1=z/r2+z/r1 ; 
sig2*z/r2-z/r1 ; 

t oi*0 . 75*sqrt (1-1 ambda A 2 )+0 . 25; 
delta1=(ks*lambda)/(pi*toe); 

de 1 1 a2=(0 . 2955*ks* I ambda A 3)/(p i *to» A 2*sqrt (1-1 ambda A 2) ) ; 
dS I dr-0 . 788*de I t al *(n1 -n2)-0 . 394* I ambda*de I t al 
*(n1 +n2)+de I ta2*( (nl -n2)- I ambda*(n1 +n2) ) ; 

Xderi votive of the stream function mrt r 
dS i dz e de 1 1 a 1 *(0 . 788*s i g2-0 . 394* I ambda*s i gl ) 

+de 1 1 a2*(s I g2- 1 ambda*s i g 1 ) ; 

^derivative of the stream function mrt z 

if (r1 A 2/(0. 1*rois A 2))>444, 
dmps=1 ; 

else ^velocity damping factor 

dmps-1 -exp(-r 1 A 2/(0 . 1 *ro«s A 2) ) ; 
end 

vxs»(x/r)*(1/r)*dSidz; Xhorizontal velocity x direction 

vzs*-(1/r)*dSidr; ^vertical velocity 

end % i f 



117 



{ I URGE LARGE RING CRLCULr T I OHS 
x e x; 

z=-ze-LUoralt; 
if x==0, 
if yL==0, 
vxLi*0; 

vzL i »-kL*RL A 2/(2*(RL A 2+z A 2) A 1 . 5) ; 
end 
else 

r*sqrt(x A 2+yL A 2); 
rl B sqrt (z A 2+(r-RL) A 2) ; 



r2*sqrt(z A 2+(r+RL) A 2); 



{radial distance front the z axis 
{closest distance to the point of 
{interest front the ring filament, 
{farthest distance to the point of 
{interest front the ring filament, 
{seal ing term 



Iambda*(r2-r1 )/(r2+r1 ); 
nl*(r+RL)/r2; 
n2*(r-RL)/r1 ; 
sigl=z/r2+z/r1 ; 
sig2=z/r2-z/r1 ; 

tom=0 . 75*sqrt (1-1 ambda A 2)+0 . 25; 
delta1*(kL*lambda)/(pi*toi); 

del ta2«(0. 2955*kL* I ambda A 3)/(p i *to« A 2*sqrt ( 1 - I ambda A 2) ) ; 
dS i dr*0 . ?88*de I tal *(n1 -n2)-0 . 391* I ambda*de i t al 
*(n1 +n2)+de I ta2*( (nl -n2)- I ambda*(nl +n2) ) ; 

{derivative of the stream function frt r 
dS i dz*de 1 1 a 1 *( 0 . 788*s i g2-0 . 391* i ambda*s i g 1 ) 

+deita2*(sig2-iambda*sigl ); 

{derivative of the stream function mrt z 
dmpL i =( 1 -exp(-r 1 A 2/(0 . 1 *ro«L A 2) ) ) ; 

{velocity damping factor 

vxLi«(x/r)*(l/r)*dSidz; {horizontal velocity in x direction 

vzLi*-(l/r)*dSidr; {vertical velocity 

end {if 
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X IflfiGE SHALL RING CALCULATIONS 
z=-ze-sUoralt; 
if x==0, 
if ys='0, 
vxsi=0; 

vzs i *-k3*Rs*2/(2*(Rs*2+z'‘2) A 1 . 5) ; 
end 
else 

r*sqrt((x+50K2+ys A 2); Xnadial distance from the z axis 

rl s *sqrt(z^2+(r-Rs)'‘2); Xclosest distance to the point of 



n1*(r+Rs)/r2; 
n2*(r-Rs)/r1 ; 
sig1*z/r2+z/r1 ; 
sig2 s z/r2-z/r1 ; 

to«-0 . 75*sqrt (1-1 ambda^2)+0 . 25; 
delta1*(ks*lambda)/(pi*to»); 

de I t a2" (0 . 2955*ks* I ambda A 3)/(p i *t om A 2*sqrt (1-1 ambda A 2) ) ; 
dS i dr*0 . 788*de I tal *(n1 -n2)-0 . 394* I ambda*de I tal 
*(n1+n2)+delta2*((n1-n2)-lambda*(n1+n2)); 



dS i dz»de 1 1 a 1 * ( 0 . 788*s i g2-0 . 394* I ambda*s i g 1 ) 

+de I ta2*(sig2- lambda’s I gl ) ; 

Xderivative of the stream function «rt z 

if (rP2/(0. 1*roms A 2))>444, 
dmpsi=1 ; 

else Xvelocity damping factor 

dmps i « 1 -exp ( -r P 2/ ( 0 . 1 *roms A 2 ) ) j 



r2=sqrt(z A 2+(r+Rs)^2); 



I ambda»(r2-r1 )/(r2+r1 ) ; 



X interest from the ring filament. 
Xfarthest distance to the point of 
X interest from the ring filament. 
Xscal ing term 



Xderivative of the stream function mrt r 



end 



vxsl»(x/r)*(1/r)*dSidz; 
vzsi B -(1/r)*dSidr; 
end X i f 



Xhorizontal velocity x direction 

Xvertical velocity 
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X SOURCE PROFILE FOR UZ CORRECTION 

Sa 1 1 * 1 0OOO j Xsource alt itude 

gagaz=2772000; Xsource strength/2» 

gagax® 1355396; Xsource strength/2w far 

uzsource B 0.0; 

uxsource=0.0; 

i f ze>800, 

uzsource*-1 5*(LUora I t-ze)/( (LUora I t-ze) A 2+x A 2) A 0 .5; 
else 

Radius1=(x A 2+(Salt-ze) A 2) A 0.5; 

Xradius fro* primary source to the flight pt 
Radius2 a (x A 2+(Salt+ze) A 2) A 0.5; 

Iradius fro* iaage source to the flight pt 
uzsource=gagaz*( (Sa I t +ze)/Rad i us2 A 2-(Sa I t-ze)/Rad i usl A 2) ; 
end 

uxsource a gagax*x*( 1 /Rad i usl A 2+1 /Rad I us2 A 2) ; 

x sunnRTiON 
If zerock”999, 
d*p s d*pL*d*ps*dapL I *daps i ; 

If x<«2000, 

ux“(dap*(uxL+uxs+uxLi+uxsi )); 
else 

ux«uxsource; 

end 

uz a -(dap*(uzL+uzs-uzL i -uzs I )-uzsource) ; 
else 

ux«- ( uxL+uxs+uxL I +uxs I ) ; 
vz*uzL+uzs-uzL I -uzs i ; 
end 
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APPENDIX B 



PARAMETER SENSITIVITY ALGORITHM 



The following are MatLab® m-file and function algorithms used to 
determine the parameters for the windshear model that 
approximates the DAL191 recorded winds. 



A. M-FILE DRIVER 

^Parameter sensitivity program for the DRL191 mind model. 

clear 

clc 

^Recorded flight data from DRL191. 

Href number, a 1 1 i tude( ft RGL), dist frm HB center(ft), x mind 
(ft/s), y mind (ft/s)]. 



DRL-C1 


987 -8612 


2 10; 


2 


975 


-8361 


-1 


11; 


3 


962 


-8080 


-3 


11; 


1 


918 


-7799 


-1 


11; 


5 


935 


-7518 


-6 


11; 


6 


922 


-7237 


-9 


15; 


7 


909 


-6956 


-12 


12; 


8 


891 


-6676 


-12 


3; 


9 


875 


-6396 


-9 


6; 


10 


656 


-6117 


-11 


9; 


11 


812 


-5838 


-16 


1; 


12 


823 


-5561 


-25 


-18; 


13 


807 


-5287 


-36 


-11; 


11 


790 


-5016 


-39 


-13; 


15 


777 


-1719 


-29 


-20; 


16 


761 


-1186 


-26 


-17; 


17 


755 


-1229 


-22 


-17; 
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18 


747 


-3977 


-29 


-1?; 


19 


712 


-3730 


-20 


-11; 


20 


738 


-3487 


-17 


-23; 


21 


730 


-3249 


-12 


-30; 


22 


714 


-3013 


-8 


-36; 


23 


698 


-2778 


0 


-34; 


21 


660 


-2542 


3 


-32; 


25 


662 


-2306 


0 


-31; 


26 


647 


-2069 


-8 


-31; 


2? 


635 


-1831 


-18 


-24; 


28 


627 


-1593 


-20 


-26; 


29 


620 


-1355 


-10 


-36; 


30 


608 


-1117 


-8 


-40; 


31 


592 


-880 


18 


-4; 


32 


585 


-639 


32 


19; 


33 


582 


-391 


19 


16; 


34 


580 


-139 


18 


21; 


35 


573 
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19 


-3; 


36 


518 


385 


38 


-21; 


3? 


511 


656 


54 


-19; 


38 


462 


936 


65 


9; 


39 


411 


1228 


66 


9; 


40 


351 


1531 


72 


-13; 


41 


277 


1815 


71 


-12; 


42 


193 


2169 


72 


-33; 


43 


119 


2501 


76 


-2; 


44 


70 


2851 


78 


15; 


45 


12 


3211 


61 


17; 


46 


20 


3579 


81 


1; 


47 


0 


3954 


85 


2] 



XParaaeter matrix: 



[2300 


2500 


2800; 


XLarge ring x displacement 


1250 


1300 


1350; 


XSaall ring core radius 


57204.9*2 


57204.9*2.3 


57204.9*2.5; 


XSaall ring cortex strength 


-400 


-300 


-200; 


XLarge ring y posit 


300 


500 


700; 


XSaall ring x displacement 


-150 


1 


150; 


XSaall ring y displacement 
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1386000*1.8 1386000*2 1386000*2.2; XSource z strength 

193628*6.5 193628*7 193628*7.5]; XSource x strength 



uxerr B 0.0;uzerr B 0.0;an“0;bestras a 999; Unit i a I izat ion 

i-i; 

z=i; 

for k*1 :3, 
hole 

i 

z 

k 

bestns 

for 1*1:3, for ■■1:3, for n«1:3, 
for o*1 :3, for p*1 :3, 
uxerr a 0 . 0 ; uzerr a 0 . 0 ; 
for c B 1 :17, 

[Ux,Uh]"p«icroburst(DflL(c,2),DflL(c,3),Ps(1 , i ),Ps(2,z), 
Ps(3,k),Ps(1, l),Ps(5,«),Ps(6,n),Ps(7,o),Ps(8,p)); 
uxerr B uxerr+(0flL(c,1)-Ux) A 2; 
uzerr«uzer r+ ( DflL ( c , 5 ) -Uh ) A 2 ; 
end 

rasux B (uxerr/40) A 0.5; 
rasuz B (uzerr/10) A 0.5; 
r»st ot ■ ( (uxerr+uzerr ) /80) A 0 . 5 ; 

■n-an+l ; 

Pf1(«n, : )"[Ps(1 , i ),Ps(2,z),Ps(3,k),Ps(1, 1 ), 

Ps(5,i),Ps(6,n),Ps(7,o),Ps(8,p),r«sox,r«svz,r«stot]; 
if bestras>r»stot, 

BESTP-Pn(tn,:); 
bestr«s“r«stot; 
saue BEST 
end 

end, end, end, end, end, end, 
saue paraaetersll 

clear ex eh rasux rastot an PD BESTP 
uxerr a 0 . 0 ; uzerr B 0 .0; an»0 ; % i n i t i a I i zat i on 
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1 - 1 ; 

z=2; 

for k=1 :3, 
hone 

i 

z 

k 

bestras 

for 1 = 1:3, for m=1:3, for n=1:3, 
for o=1 :3, for p=1 :3, 
vxerr*0 . 0 ; vzerr=0 . 0 ; 
for c=1 :47, 

[Ux,Uh]*paicroburst(DRL(c,2),DflL(c,3),P8(1 , i),P«(2,z), 
Ps(3,k),Ps(4, l),P3(5,i),P3(6,n),p3(7,o),Ps(8,p)); 
vxerr=vxerr+(DRL(c,4)-IJx) A 2; 
vzerr*vzerr+(DRL(c,5)-Uh) A 2; 
end 

ra3vx=(vxerr/40) A 0 .5; 
ra3vz“(vzerr/10) A 0.5; 
ra3tot*((vxerr+vzerr)/80) A 0.5; 
an*an+1 ; 

Pn(m, : )=[Ps(1 , i ),P8(2,z),Ps(3,k),Ps(4, 1 ), 

P3(5,i),P3(6,n),P3(7,o),P3(8,p),r«3ux,r«3uz,r*3tot] 
if be8tr»8>r»8tot, 

BESTP=Pn(m,:)j 
bestras^rastot; 
save BEST 
end 

end, end, end, end, end, end, 
save paraaeters12 

clear «x «h r»svx nstot an PH BESTP 

vxerr=0.0;vzerr»0.0;an*0; Unit i a I izat ion 

M; 

z-3; 

for k=1 :3, 
hoae 



124 



z 

k 

best ms 

for 1*1:3, for #=1:3, for n=1:3, 
for o*1 : 3, for p*1 :3, 
uxerr*0 . 0 ; ozerr=0 . 0 ; 
for c=1 :17, 

[Ux,Uh]=paicroburst(DRL(c,2),DfiL(c,3),Ps(1 , i ),Ps(2,z), 
Ps(3,k),Ps(1, I ),Ps(5,#),Ps(6,n),Ps(7,o) ,Ps(8,p)); 
oxerr=yxerr+(DflL(c,1)-Ux) A 2; 
uzerr*uzerr+(DRL(c, 5)-Uh) > '2; 
end 

r#3ux=(vxerr/40) A 0 ,5; 
r#syz*(uzerr/40) / '0 .5; 
r#stot*((uxerr+uzerr)/80) / '0.5; 

#n*#n+1 ; 

PfKan, : )*[Ps(1 , i ),Ps(2,z),Ps(3,k),Ps(4, 1 ), 

Ps(5,i),Ps(6,n),Ps(7,o),Ps(8,p),r#sux,r#suz,r#stot]; 
if bestr#s>r#stot, 

BESTP*Pf1(an, : ); 
be3trie*r#stot; 
saue BEST 
end 

end, end, end, end, end, end, 
saue parameters 13 

clear «x «h rasux ristot an PH BESTP 

uxerr*0.0;uzerr=0.0;an*0; Unit i a I Izat Ion 

i*2; 

z*1; 

for k=1 :3, 
hoae 
i 

z 

k 

bestras 

for 1*1:3, for #=1:3, for n=1:3, 
for o*1 :3, for p*1 :3, 
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uxerr=0 . 0 ; yzerr c O . 0 ; 
for c=1 : 47, 

[Ux,Uh]*p»icpobupst(DflL(c,2),DflL(c,3),P3( 1 , i ),Ps(2,z), 
Ps(3,k),Ps(4, l),P3(5,a),Ps(6,n),P3(7,o),Ps(8,p)) ; 
vxerr=uxerr+(DflL(c,1)-Hx) A 2; 
vzepp=vzenp+(DflL(c,5)-llh) A 2; 
end 

pasux*(vxepp/40K0 .5; 
pi3vz=(vzepp/10) *0 . 5 ; 
past ot = ( ( uxepp+uzepp ) /80 ) A 0 . 5 ; 
an*an+1 ; 

PIKan, :)*[Ps(1, i),P3(2,z),Ps(3,k),Ps(4, 1), 

Ps(5,B),P3(6,n),P3(7,o),P3(8,p),pasvx,pni3uz,pastot]; 
if be3tpas>PB3tot > 

BESTP«Pn(an,:); 

be3tPi3»PB3tot ; 

3ave BEST 
end 

end, end, end, end, end, end, 
saue papaBeteps21 

clean «x «h pbsox pastot an PM BESTP 

uxepp*0.0;yzepp*0.0;an*0; Unit i a I izat ion 

i-2; 

z-2; 

fop k*1 :3, 
hoae 



k 

best pbs 

fop 1*1:3, fop b* 1:3, fop n*1:3, 
fop o*1 : 3, fop p*1 :3, 
yxepp*0 .0; uzepp*0 .0; 
fop c*1 :47, 

[Ux,Uh]*pBicpobupst(DflL(c,2),DflL(c,3),P3(1 , i ) , Ps(2, z) , 
Ps(3,k),P3(4,l),p3(5,B),P3(6,n),Ps(7,o),P3(8,p)); 
yxepp*uxepp+(DflL(c,4)-Ux) A 2; 
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ozerr a uzerr+(DRL(c,5)-Uh) A 2; 

end 

rasux* ( uxerr/40 ) A 0 . 5 ; 
rasvz*(uzerr/40) A 0.5j 
rastot*((uxerr+yzerr)/80) A 0.5; 

»n=an+l ; 

Pn(in,:MPs(1,i),P3(2,z),Ps(3,k),Ps(4,l), 

P3(5,»),P3(6,n),Ps(7,o),P3(8,p),ra3ux,rB3uz,rB3tot ] ; 
i f bestras>ra3tot, 

BESTP*PM(an, : ); 
bestrBs=rastot; 
saue BEST 
end 

end, end, end, end, end, end, 
saue paraoeters22 

clean ox oh rasux rastot on PD BESTP 

uxerr«0.0juzerr*0.0;an*0; X i n i t ial izat Ion 

i-2; 

z-3; 

for k*1 :3, 
hoae 
I 

z 

k 

bestras 

for 1*1:3, for a*1:3, for n*1:3, 
for o*1 :3, for p*1 :3, 
uxerr*0 . 0 ; uzerr*0 .0; 
for c*1 :47, 

[Ux , Uh]*pa i croburst (DRL(c, 2) , DRL(c, 3) , Ps( 1 , i ) , Ps(2, z) , 

p3(3,k),Ps(4, l),P3(5,a),Ps(6,n),p3(7,o),P3(8,p)); 
vxerr*vxerr+(DRL(c,4)-Ux) A 2; 
uzerr*vzerr+ (DRL ( c , 5 ) -llh ) A 2 ; 
end 

rasvx* ( uxerr/40 ) A 0 . 5 ; 
rasuz*(vzerr/40) A Q.5; 
ra3tot*((uxerr+vzerr)/80) A 0.5; 
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nn=nn+l ; 

PM(nn, : )*CPs(1 , i),Ps(2,z),Ps(3,k),P3(4, I), 

Ps(5,n),Ps(6,n),P3(7,o),P3(8,p),rnsvx,riJisuz,riii3tot ]; 
if bestrins>rBi8tot, 

BESTP=PH(nn, : ); 
be3tr«s=rn3tot; 
saue BEST 
end 

end, end, end, end, end, end, 
saue paraneter323 

clear t»x #h rnsux r»stot nn PD BESTP 

uxerr*0.0;uzerr*0.0;qin=0j X in i t ial izat ion 

i=3; 

z=1; 

for k*1 :3, 
hone 
I 

z 

k 

best ms 

for 1*1:3, for n*1:3, for n*1:3, 
for o*1 :3, for p*1 :3, 
vxerr*0 . 0 ; uzerr*0 . 0 ; 
for c*1 :17, 

[Ux,Uh]*pnicroburst(DRL(c,2),DflL(c,3),P3(1 , i ) , Ps(2, z ) , 
Ps(3,k),p3(1, l),P3(5,n),Ps(6,n),Ps(7,o),P3(8,p)); 
uxerr*vxerr+(DRL(c, 4)-Ux) A 2 ; 
ozerr*uzerr+(DRL(c,5)-Uh) A 2; 
end 

rn3vx*(vxerr/40) A 0.5; 
r«3vz*(vzerr/40) A 0 , 5 j 
ri3tot*((yxerr+vzerr)/80) A 0.5; 

■n*nn+1 ; 

Pn(nn,:)-[P3(1,i),Ps(2,z),P3(3,k),Ps(1,l), 

Ps(5,n),Pe(6,n),P3(7,o),P3(8,p), rnsux, rn3uz,rn3tot]j 
if be3trn3>rn3tot, 

BESTP-PfKnn, : ) ; 
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be3tr»3*ri3tot ; 
saue BEST 
end 

end, end, end, end, end, end, 
saue parameter331 

clear #x §h rasux rastot an Pfl BESTP 

uxerr a O.O;uzerr s O.O;an at O; Xinit i a I izat ion 

i=3; 

z“2j 

for k=1 :3, 
hoae 

I 

z 

k 

best r #3 

for 1*1:3, for a a 1:3, for n a 1:3, 
for o*l :3, for p*l :3, 
uxerr a O . 0; uzerr a 0 . 0 ; 
for c*l :47, 

[Ux,Uh]«p«icrobur3t(DRL(c,2),DflL(c,3),P3(1,i),P3(2,z), 

Ps(3,k),Ps(4, l),Ps(5,«),P3(6,n),p3(7,o),P3(8,p)); 
uxerr a uxerr+(DRL(c,4)-Ux) A 2; 
uzerr a uzerr+(DRL(c,5)-Uh) A 2j 
end 

ra3ux a (uxerr/40) A 0 .5; 
rasuz a (uzerr/40) A 0 . 5 ; 
ra3tot a ((uxerr+uzerr)/80) A 0.5; 
an a an+l ; 

PH(an, : )*[Ps(l , i ),Ps(2,z),P3(3,k),Ps(4, I ), 

p3(5,a),P3(6,n),P3(7,o),P3(8,p),r«3vx,ra3uz,r«3tot] 
if bestr«3>r«3tot , 

BESTP*Pn(»n, : ) ; 

bestr»3*r«3tot ; 
saue BEST 
end 

end, end, end, end, end, end, 
saue paraaeters32 
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clean »x <nh rasux rastot an PM BESTP 



uxerr=0.0;uzerr=0.0;an=0; % i n i t ializat ion 

i=3; 

z*3; 

for k=1 :3, 

hoae 

I 

z 

k 

bestras 

for 1=1:3, for a*1:3, for n=1:3, 
for o=1 :3, for p=1 :3, 
uxerr«0 . 0 ; ozerr=0 .0; 
for c«1 '. 47 , 

[Ux,Uh]=paicrobur3t(DRL(c,2),DRL(c,3),P3(1 , i ),Ps(2,z), 
P3(3,k),Ps(4, l),P3(5,a),P3(6,n),P3(7,o),P3(8,p)); 
uxerr *uxerr+ (DRL(c, 4 )-Ux) A 2; 
uzerr*vzerr+(DRL(c,5)-Uh) A 2; 
end 

ra3ux=(vxerr/40) A 0.5; 

ra3uz=(vzerr/40) A 0.5; 

ra3tot=((uxerr+uzerr)/80) A 0.5; 

•n=an+1 ; 

PH(an, :)«[Ps(1, i ),Ps(2,z),P3(3,k),P3(4, I), 

P3(5,a),Ps(6,n),P3(7,o),P3(8,p),ra3vx,ra3vz,ra3tot] 
if bestras>ra3tot, 

BESTP=Pfl(an, : ); 
bestra3“ra3tot; 
saue BEST 
end 

end, end, end, end, end, end, 

save paraaeters33 

clear 

load BEST 

uxerr«0 .0; vzerr»0 . 0 ; 
for i»1 :47, 
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hone 

i 

X S DRL( i ,3); 

[Ux,Uh]=p»icpobupst (DRL( i ,2),DRL( i , 3) , BESTP( 1 ),BESTP(2), 

BESTP(3),BESTP(-») ,BESTP(5),BESTP(6) J BESTP(?),BESTP(8)); 
node I ( i , : )=[ i ,lJx,Uh]; 
uxern=uxern+(DRL( i ,4)-tlx) A 2; 
uzerr*uzern+(DRL( i , 5)-Uh) A 2; 
end 

rnsox=(oxenr/40) A 0.5; 
rnsyz=(uzenp/40) / '0 .5; 
rnstot s ((uxerp+ozerr)/80) A 0.5; 
ansiep*0; 

•hile an3»ep"®999, 
dispCThe hopizontal RMS is: ’) 

disp(pnsvx) 

dispCThe veptical RMS is:') 
disp(pBsoz) 

dispCThe total systen RRS is:’) 
disp(pBstot) 

dispCThe foil on i ng plots ape available:') 
dispC 1»hopizontal ninds') 
dispC 2*veptical »inds’) 
dispC 999»end’) 
ansnep* input ('plot nunbep '); 
if ansiep"*1, 
p I ot (0RL( : , 3) , DRL( : , 4) ) 
hold 

plot(DRL( : ,3),»odel ( : ,2), '--’) 

pause 

hold 

el seif ansnep M 2, 
plot(DRL(:,3),DRL(:,5)) 
hold 

plot (DRL( : ,3),nodel ( : ,3), ' — ' ) 
pause 
hold 
end, end 
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B. WINDSHEAR FUNCTION 



funct i on [vx, vz]-pa i cpobupat (ze, x, Lxd i ap, Ra, ka, yL, axd i ap, RL, kL) 
X INITIAL URLUES RND INPUTS 



X Eapth poait Ion 
Xze * alt itude (ft) 

Xx hoplzontal poait ion (ft) 




X Uoptex Ring Diaenaiona 


XRL 


Xpadiua of ping fl laaent (ft) 


pIL«2004.1j 


Xcopo padiua (ft) 


XkL 


Xuoptex atpength of the ping 
Xfi laaent (ft A 2/aec) 


po*L-piL/0.3?1 ; 


Xpelated cope padiua 


*yL 


Xlapge ping y latepial poait ion (ft) 


LUopalt-3400; 


Xlapge voptex altitude RGL (ft) 


XLxdiap 


Xlapge ping diaplaceaent in x (ft) 


XRs 


Xpadiua of ping f i laaent (ft) 


p i a-323. 9; 


Xcope padiua (ft) 


Xka 


Xvoptex atpength of the ping 
Xfi laaent (ft A 2/aec) 


poaa-pia/0.371 ; 


Xpelated cope padiua 


ya*830.9/830.9; 


Xaaall ping y latepial poait ion (ft) 


aUopalt-800; 


Xaaall uoptex altitude RGL (ft) 


Xaxdiap 


Xaaall ping diaplaceaent in x (ft) 



X PRIURRY LARGE RING CRLCULRTIONS 

z«ze-LUopalt; 

zepock-999; 

p-aqpt ( (x-Lxd I ap) A 2+yL A 2) ; 
pi -aqpt (z A 2+(p-RL) A 2) ; 
p2-aqpt (z A 2+(r+RL) A 2) ; 

Iaabda«(p2-p1 )/(p2+p 1 )j 
n1«(p+RL)/p2j 
n2«(p-RL)/pl ; 
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3ig1*z/r2+z/rl ; 

8ig2=z/r2-z/r1 ; 

to#=0 . 75*sqrt (1-1 a»bdcT2)+0 . 25; 
delta1=(kL*la»bda)/(pi*to«») ; 

de I ta2=(0 . 2955*kL* I anibda A 3)/(p i *t o» A 2*sqrt (1-1 aibda A 2) ) ; 
dS i dr=0 . 788*de 1 1 a 1 *(n1 -n2) -0 . 394* I a»bda*de 1 1 al 
*(n1+n2)+de Ita2*((n1-n2)-l aibda*(n1 +n2) ) ; 
dS i dz=de I tal *(0 . 788*s i g2-0 . 394* I aibda*3 i gl ) 
+delta2*(sig2-laiabda*3ig1 ); 
d»pL=(1 -exp(-nP2/(0 . 1 *ro»L A 2) ) ) ; 
uxL s (x/r)*(1/r)*dSidz; 
uzL=-(1/r)*dSidr; 

* PRIIIfiRV SHALL RING CRLCULRTIONS 
z=ze-sUoralt; 

r=3qrt ( (x-sxd i sp) A 2+ys A 2) ; 
r 1 =3qrt (z A 2+(r-Rs) A 2) ; 
r2*3qrt (z A 2+(r+R3) A 2) ; 

Ia»bda=(r2-r1 )/(r2+r1 ) ; 
n1*(r+Rs)/r2; 
n2*(r-Rs)/r1 ; 
sig1*z/r2+z/r1 ; 
sig2*z/r2-z/r1 ; 

toi*0 . 75*sqrt (1-1 a»bda A 2)+0 . 25; 
delta1'(k3*laibda)/(pi*to«); 

de I ta2=(0 . 2955*ks* I aibda A 3)/(p i *toi A 2*sqrt (1-1 aibda A 2) ) ; 
dS i dr=0 . 788*de 1 1 a 1 *(n 1 -n2)-0 . 394* I a»bda*de 1 1 a 1 
*(n1+n2)+delta2*((n1-n2)- Ia»bda*(n1+n2)); 
dS i dz=de 1 1 a 1 * ( 0 . 788*3 i g2-0 . 394* I aibda*3 i g 1 ) 
+delta2*(3lg2-la»bda*3lg1 ); 
if (r1 A 2/(0. 1*rois A 2))>444, 
d»ps B 1 ; 
else 

dips* 1 -exp(-r 1 A 2/(0 . 1 5*rois A 2) ) ; 
end 

ux3»(x/r)*(1/r)*dSidz; 

uz3=-(1/r)*dSidr; 
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X IfIRGE LRRGE RING CRLCULRT I OHS 
z a -ze-LUopalt ; 



p=sqpt ( (x-Lxd i sp) A 2+yL A 2) ; 
p1=3qpt(z A 2+(p-RL) A 2); 
n2=sqrt (z a 2+(p+RL) a 2) ; 

Ia#bda=(r2-r1 )/(p2+p1 ); 
n1»(r+RL)/r2; 
n2 a (r-RL)/r1 ; 
sig1*z/p2+z/p1 ; 
sig2=z/p2-z/p1 ; 

t o»»0 . 75*sqpt (1-1 a»bda' s 2)+0 . 25; 
delta1«(kL*la«bda)/(pi*toi); 

de Ita2*(0. 2955*kL* I a*bda A 3)/(p i *to«T2*3qpt (1-1 a«bda A 2) ) ; 

dSidp=0.788*deltal*(n1-n2)-0.391*la«bda*delta1*(n1+n2)+delta2*((n1-n 

2)-laibda*(n1+n2)); 

dS i dz»de 1 1 a 1 *(0 . 788*s i g2-0 . 391* I aibda*s i g 1 ) 

+delta2*(sig2-lanbda*3ig1 ) ; 
dipL i -( 1 -exp(-p1 A 2/(0 . 1 *po»L a 2) ) ) ; 
uxLi*(x/p)*(1/p)*dSidz; 
vzLi“-(1/p)*dSidp; 

X IfIRGE SdRLL RING CRLCULRT IONS 

z--ze-3Uopalt; 

p*sqpt ( (x-sxd i sp) A 2+ys A 2) ; 

p1=sqpt(z A 2+(p-R3) A 2); 

p2=sqpt(z A 2+(p+R3) A 2); 

I aibda«(p2-p1 )/(p2+p1 ) ; 
n1“(p+Rs)/p2; 
n2*(p-Rs)/p1 ; 
sig1=z/p2+z/p1 ; 
sig2«z/p2-z/p1 ; 

ton-0 . 75*sqpt (1-1 aibda A 2)+0 . 25; 
delta1-(k3*la«bda)/(pi*toi) ; 

de I ta2=(0 . 2955*ks* I a«bdcT3)/(p i *to« A 2*sqpt (1-1 a*bda A 2) ) ; 
dS i dp-0 . 788*de I tal *(n1 -n2)-0 . 391* I aibda*de Ital 
*(n1+n2)+delta2*((n1-n2)- Iaibda*(n1+n2)); 
dS i dz-de I t al *(0 . 788*3 i g2-0 , 391* I aibda*3 I g 1 ) 
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+delta2*(sig2-la#bda*sig1 ); 
if (pP2/(0. 1*ro«is A 2))>m, 
d*psi s 1 ; 
else 

d®psi=1-exp(-r1 A 2/(0. 15*ro»s A 2)); 
end 

vxsi=(x/r)*(1/r)*dSidz; 

uzsi»-(1/r)*dSidr; 

* SOURCE PROFILE FOR UZ CORRECT I OH 

SalMOOOO; Xsource alt I tude 

Xgagaz Xsource strength/2i 

Xgagax Xsource strengthen far 

vzsource=0.0; 

uxsource s O.O; 

if ze>800, 

uzsource--15*(LUoralt-ze)/((LUoralt-ze) A 2+x A 2) A 0.5; 

else 

Radius1«(x A 2+(Salt-ze) A 2) A 0.5; 

Rad i us2«(x A 2+(Sa I t+ze) A 2) A 0. 5; 

uzsource B gagaz*( (Sa I t +ze)/Rad i us2 A 2-(Sa I t-ze)/Rad I usl A 2) ; 
end 

uxsource s gagax*x*( 1 /Rad i usl A 2+1 /Rad i us2 A 2) ; 

x sunnfiTioH 

d»p“d»pL*dips*d«pU*d«psi ; 
if x<«2000, 

ux»(d«p*(uxL+uxs+vxL I +uxs i ) ) ; 
else 

vx*vxsource; 

end 

vz“-(d«p*(uzL+uzs-uzL I -uzs I )-uzsource) ; 
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APPENDIX C 



LANDING APPROACH ALGORITHMS 



Contained in this appendix are the MatLab® programs utilized in 
calculating the aircraft response given a windshear and a particular 
escape maneuver. They are listed in three sections. The first section 
contains the driver routine and the functions required to solve the 
differential equations. The second and third sections contain the 
specific aircraft parameter functions and escape maneuver functions 
respectively. The windshear function is listed in Appendix A. 



A. DRIVER AND DIFFERENTIAL EQUATION SOLVER 

g***************************** 

X Descent Profile Hodel * 

g***************************** 

X This prograv el 1 1 calculate the inertial flight path angle for 
X a given aircraft, it utilizes the function 'vlcroburst' thlch 
X aodels the aind shear encountered during the DFU Delta Accident. 

X This prograv uses energy height theory to calculate inertial space 
position. 

clear 

cic 

foraat bank 
forvat compact 

X 

liNITIRLIZRTION AND ENU I RONflENTAL CONSTANTS 

% The foi loving initialize constants and variables for the vodei. 
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Xfl I RCRRFT CONSTANTS 

X The following constants are for the particular aircraft being 
studied. 

dispC LANDING APPAOACH ttODEL * ) 

dispC ') 

answer=0.0; 

disp( 'Type aircraft .* ’ ) 
dispC I. P3 at 89,5001b') 

disp( ' 2. P3 at 1 M, 5001b' ) 

dispC 3. P3 at 89,5001b 3 engine') 

dispC 4. L101 1 ' ) 

dispC 5. T41' ) 

answer* i nput ( ’ Enter desired aircraft by nuaber: '); 

If answer**1, 

[S,K,CDo,Taax,Ut,CLo,CLalpha,alphaaax,q,U1 ,T,theta]=P3LtApp; 

XP3 § 89,500lbs, app config. 

el seif answer*=2, 

[S,K,CDo,Taax,Ut,CLo,CLalpha,alphaaax,q,Ul ,T,theta]*P3HvApp; 

XP3 e 1M,500lbs, app config. 

el seif answer==3, 

[S,K,CDo,Twax,Ut,CLo,CLalpha,alpha«ax,q,U1 ,T,theta]=P3Eout ; 

XP3 with one engine out. 

el seif answer**^, 

[S,K,CDo,Taax,Ut,CLo,CLalpha,alphaaax,q,Ul ,T, theta ]*L1 01 1 App; XL1 01 1 
el self answer**5, 

[S,K,CDo,Taax,Ut,CLo,CLalpha,alphaaax,q,U1 ,T,theta]*T14App;XT44 
end 

XS*plane fora surface area In ft A 2 
XK*coef f icient for CD calculation 
XCDo*coef f icient of drag 
XTaax s aax thrust auailable in ft-lbs 
XUt*aircraft weight in lbs 
fl*Ut/32.1?4; X aass in slugs 
XCLo*zero A0A lift coefficient 
XCLalpha*l i ft curve slope 
Xalphaaax*stal I buffet alpha w/ approach flaps 
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lq s nominai pitch rate of 5deg/sec 
IU1=approach equivalent airspeed in ft/s 
JKT® initial thrust in ibf 



I 



l#*#p r0 g raniet j initial values**************************************** 

iters*30; Inumber of iterations 

deltat'1.0; Itime increment in sec 

nBcenter s 10300; Idistance fIB center is 

from runway 

xstart=-500; 1st art ing distance from 

MB center 

g******************************************************************* 



hstart B (f1Bcenter-xstart)*tan(0.052)j Istarts on a 3deg 

gl idesiope 

dispC Initial Input Section') 
dispC ') 

dispC The following values are preprogrammed: ' ) 
dispC a. The approach is based upon a 3deg gi idesiope. 1 ) 
dispC b. The starting distance (distance from HB center) is:') 
disp(xstart) 

dispC c. The fIB center is this from the end of the runway:') 
disp(riBcenter) 

dispC d. The starting altitude is:’) 
dlsp(hstart) 

dispC e. The time step in sec:') 
disp(deitat) 

dispC f. The number of iterations is:') 
dispC iters) 
dispC ') 



TU B Twax/Ut; It hrust /weight ratio 

UL*Ut/S; Iwing ioading 

dispCThe following values were suppi ied/calcuiated from aircraft 
constants' ) 

dispC a. Uref in ft/s (1 .35Uso+5kts for P3 and l.3Uso for other 

aircraft : ' ) 

dlsp(UI) 

dispC b. Thrust to weight is:’) 
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disp(TU) 

dispC c. Uing loading is:') 

disp(UL) 

dispC ') 

■aneuver=0.0; 

dispC Select the escape saneuver:') 



dispC 1. 
sindshear' ) 


Stay on 


glide 


slope, go 


■issed at 


200ft 


R6L; NO 


dispC 2. 
sindshear' ) 


Stay on 


glide 


slope, go 


sissed at 


200 ft 


RGLjUiTH 



dispC 3. Escape slth constant airspeed') 
dispC 4. Escape sith constant altitude') 
dispC 5. Escape sith constant theta') 
dispC 6. Escape sith constant alpha') 
dispC 7. Escape sith sax alpha') 

dispC 8. Esuiate DRL191 final fes seconds using recorded theta’) 
dispC 9. Esuiate DRL19I final fes seconds using recorded theta and 
sinds' ) 

saneuver® inputC Enter selection nusber: '); 

if saneuver BB 5, 

dlspCUhat value for escape theta?') 
thetaais B input( ' Enter in radians: ’); 
el seif saneuver BB 6, 
dlspCUhat value for escape alpha?') 
aiphaais B input( ’ Enter in radians: '); 
el seif saneuver BB 7, 
aiphaais B aiphasax; 
el self saneuver > B 8, 
hstart=hstart+50; 
end 

x1 B xstart ; 
h1 B hstart ; 

Ug-Ul ; 

gasal«-.052; 
rad (-3deg) 



Xinlt ial Izat ion 
Unit ial izat ion 

Xground speed initialization in ft/s 
X in it ial inertial flight path angle in 
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poifO. 002377; 
tr igger=0.0; 
check1 a 0.0; 
check2 a 0.0; 

X 

XflLGORITHn 



lair density in slugs/ft A 3 § S.L. 
Xutiiity check number 
lut i I ity check nutber 
Xutiiity check number 



X The actual algorithm is broken into seueral smaller dioisions and 
functions: 

X starting ualues - as it implies 
X loop > beginning of iterating process 
X nindshear inputs - uses the function ’microburst' to giue 
X point x and z axis mind speeds 

X predictor - first guess at the solution of the gouerning OEs 

X corrector - corrected solutions to the gouerning DEs and 

X final point ualues. 

X output file management - as it implies 
X aircraft control - uses multiple control functions for 
X control inputs to theta, thrust, and ROR 

X final output - as it implies 

X 



Xstart ing ualues: 
pt-1; 

Q«0.5*rom*Ur2; 

CL-Ut/(Q*S); 

CD«CDo+K*CL A 2; 

a I pha-(CL-CLo)/CLa I pha j 

gamaa-theta-alpha; 

Es«Ur2/(2*32.174)+h1; 

Esdot«-l2.32j 



Xcounter 

Xdynamic pressure 

Xcoeff icient of i i ft 

Xcoefficient of drag 

XROR required for 1g flight 

Xairmass flight path angle 

Xspecific energy in ft 

Xtime rate change in specific energy 



y(1)-x1; 

y(2)-h1; 

y(5)«U1; 

y(6)-gamaa; 



Xy(1 )*»x position in ft 
Xy(2)»aircraft altitude in ft 
Xy(5)*aircraft airspeed In ft/s 
Xy(6)«aircraft airmass fit. path (rad) 



Aircraft (1 , : Mxl ,h1 ,U1 ,-12.36, theta, a I pha, 1 , T] ; 
inert ial (1 , : )“[Ug,0,gamaa,gamai ,0,0, 1 ,Es,Esdot]; 

X start ing output files 

X 



MO 



XI oop control 
Xpresent pt in space 



X I oop : 
clc 

mhi le theta~*999, 
pt°pt+1 ; 

% 

Xmind shear inputs: 
x2=Ug*deltat+x1 ; Xapproximate horizontal displacement 

h2=tan(gamai )*(x2-x1 )+h1 ; Xapproximate uertical displacement 

if maneuver” 1 , 

Uxl 1 *0; Ux21 =0; Uxl 2=0 ; Ux22=0 ; Uhl 1 =0; Uh21 =0; Uh 1 2*0 ; Uh22=0 ; 
el seif maneuver==9, 

[Uxl 1 ,Uh1 I ]=recminds(pt); Xrecorded minds from DRL191 
[Ux22,Uh22]=recminds(pt+1 ); 

Ux12=(Ux1 1+Ux22)/2; 

11x21=11x12; 

Uh12=(Uh1 1+Uh22)/2; 

Uh21=Uh12; 

else 

[Ux11,Uh11]=microburst(h1,x1); XThese four calls to the function 
[Ux21,Uh21]*microburst(h2,x1); Xmicroburst finds the point mind 
[Ux12,Uh12]=microburst(h1 ,x2); Xshear for determination of the 
[Ux22,Uh22]*microburst(h2,x2); Xauerage block mind shear. 

XThe above is used in conjunction 
Xother calls mhen no mind shear is 
Xdesired. 
end 

UindX(pt)=Ux1 1; XUind shear matrices 

UindH(pt)=Uh1 1 ; 

dUxdx=( (Uxl 2-Uxl 1 )+(Ux22-Ux21 ) )/(2*(x2-x1 ) ) ; XThese formulas 
dUxdh=((Ux22-Ux12)+(Ux21-Ux1 1 ))/(2*(h2-h1 ));Xcalculate the 
dUhdx*((Uh12-Uh1 1 )+(Uh22-Uh21 ))/(2*(x2-x1 ));Xdi fferent ial change 
dUhdh=((Uh22-Uh12)+(Uh21-Uh1 1 ))/(2*(h2-h1 ));Xin mind shear. 

Xy(3)=mind shear in the x dir. ft/s 
Xy(4)=mind shear in the z dir. ft/s 



y(3)=Ux11; 

y(4)=Uh11; 

X 



141 



{predictor and intermediate values 
ydot=ypri»es(y,dUxdx,dUxdh,dUhdx,dUhdh > theta,T,S,n,Q,CL,CD, 
alpha, a I phamax) ; 

{cal I to the OE function 



for n=1 :6, 

yl (n)®y(n)+deltat*ydot (n); {Euler first step fed predictor 
yhalf(n)»y(n)+deltat/2*ydot(n); {Euler half step fed predictor 
end 



alpha®theta-y(6) ; 
CL=CLo+a I pha*CLa I pha ; 
CD«CDo+K*CL A 2; 
Q=0.5*roe*yhal f(5) A 2; 

X 



XflOfl 

Icoeff icient of lift 
Xcoefficient of drag 
Xdynaeic pressure 



{corrector and nee value formulation 
ydot*ypr imes(yhal f,dHxdx,dl)xdh,dHhdx,dHhdh,theta,T,S,f1,Q,CL,CD, 
a I pha, a I phamax); 



for n®1 ;6, 

y2(n)®yhal f(n)+deltat/2*ydot(n);XRichardson extrapolat ion 
y(n)»2*y2(n)-y1(n); {corrector scheme, 

end 



x1»y(1); 
hi «y(2) ; 

Ug-y(5)*cos(y(6))+y(3); 
R0C»y(5)*sin(y(6))+y(4); 
gama i *atan(ROC/Ug) ; 
Es«y(5) A 2/64.348+y(2); 



{nee x posit in ft 
{nee alt itude in ft RGL 

{ground speed in ft /sec 
{Rate of Climb in ft/sec 
{nee airmass FP angle in rad 
{nee specific energy 



Esdot*y(5)*ydot(5)/32, 174+ydot(2);{nee trc of specific energy 



alpha*theta-y(6) ; 
CL*CLo+a I pha*CLa I pha; 
CD«CDo+K*CL A 2; 
Q«0.5*roe*y(5) A 2; 

X 



{nee ROR 

{nee coefficient of lift 
{nee coefficient of drag 
{nee dynamic pressure 
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lout put file lanageient: 

Aircraft(pt, : )=[y(1 ) ,y(2),y(5), ydot (5), theta, alpha, pt , T 3 ; 

Inert ial (pt , : ) = [Ug,R0C,y(6),gaiai ,y(3),y(1) ,pt,Es,Esdot]; 

hoie 

dispC pt x h U 

airspeed accel ’) 

opl*[pt, Aircraft (pt , 1 ) ,Rircraft(pt,2), Aircraft (pt, 3), 

Aircraft(pt,1)]; 
disptopl ) 

dispC pt theta alpha gaiaa gaiai') 

op2=[pt,Rircraft(pt ,5), Aircraft (pt, 6), Inert ial (pt,3), 

Inert ial (pt, 1)]; 
dlsp(op2) 

dispC pt Ug ROC Thrust') 

op3*[pt, Inert ial(pt, 1 ), Inert ial (pt,2), Aircraft (pt,8)]; 
disp(op3) 

dispC pt Es Esdot') 

op5°[pt, Inert ial (pt, 8), Inert ial (pt, 9)]; 

disp(op5) 

dispC pt Ux Uh ’ ) 

op4*[pt, Inert ial (pt, 5), Inert Ial (pt, 6)]; 
disp(op4) 

% 



laircraft control: 

^functions called for theta, thrust 
land alpha control. 

if laneuuer <* 2, 
if y(2)<200, 
trigger *! ; 

[theta,T]*f1flcontrol(q,deltat,Tiax,theta,T,alpha,alphaiax,y); 
elsei f trigger ~* 1 , 
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[theta > Tl^GPcontro I (q^ de 1 1 at , ROC, Tatax^ theta, T, a I pha, a I phanax > U1 , 

hstart,xstart,y); 

elsei f trigger == 1 , 

[theta,T]-MRcontrol (q, del tat ,T«ax,thet a,T, alpha, a I pha»ax,y); 
end 

elsei f »aneuyer<8, 
if y(5) < (U1-33.33), 
trigger-1 ; 
end 
else 

if y(2)<-1004, 

[theta,T,checkl]-DRLtheta(T»ax,theta,T,U1 ,y, check 1 ); 
else 

[theta,T]-GPcontrol (q, del tat ,R0C,T«ax,theta,T,alpha,alpha»ax,U1 1 

hstart,xstart,y); 

end 

end 

if naneuoer >*3, 
i f trigger *= 1 , 
if taneuyer *■ 3, 

[ t het a,T, check 1 ]*CRSesc(q,deltat,Tiax,theta,T,alpha,alphaiiax, 

U1 ,y,checkl ) 

el seif «aneuver ■■ 1, 

[theta,T,check1 ]-CRLTesc(q,deltat ,T«ax,theta,T,alpha,alphaiax, 

y, check) ,ydot) 

el seif ianeuuer ■» 5, 

[theta,T]-CTHesc(q,deltat,Tiax,theta,T,alpha,alpha»ax,thetaai«) 
el seif «aneuver -= 6, 

[ t het a , T ] -CROflesc ( q , de 1 1 at , T tax , t het a , T , a I pha , a I phaeax , a I phaa i n ) 
el seif ■aneuver — ?, 

[theta,T]-f1RXesc(q,deltat ,T«ax,theta,T,alpha,alphaaii) 
end 
end 
end 

if pt— iters, Xnutber of iterations 

theta-999; 
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Xstop at ground itpact 



elsei f y(2)<0, 
theta°999; 
else 
end 

end Xwhile algoritha loop 

X - 

Xoutput files: 

answer 3 input ('Do you want a flight path plot? 1=VES ’); 
if ans»er==1 , 

plot(Rircraft( : , 1 ), Ri rcraf t( : ,2)) 
pause 

answer 3 0.O; 

end 

answer 3 input ('Do you want x, h, U, and airspeed accel? 1 3 VES '); 
if answer 33 1, 

dispC pt x h U airspeed accel ') 

op1*[Rircraft(:,7),Rircraft(: , 1 ), Aircraft (: ,2), Ri rcraft (: ,3), 
Rircraft(:,4)]; 
disp(opl) 
else 
end 

answer*lnput( 'Do you want theta, alpha, gaaaa, and gaaai? 1 *VES '); 
if answer 33 !, 

dispC pt theta alpha gaaaa gaaai') 

op2«[R I rcraf t ( : , 7) , R i rcraft ( : , 5) , R i rcraf t ( : , 6) , I nert i a I ( : , 3) , 

Inert ial ( : ,1)]; 
disp(op2) 
else 
end 

answer 3 Input ('Do you want Ug, ROC, and Thrust? 1*VES '); 
if answer 33 !, 

dispC pt Ug ROC Thrust’) 

op3 3 [Rircraft(: ,7), Inert ial (:,!), Inert ial ( : , 2) , fl i rcraft (:, 8)]; 
disp(op3) 
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else 

end 

answer* input( 'Do you want Es and Esdot? 1=VES ’); 
if answer**1, 

dispC pt Es Esdot') 

op5=[flincnaft( : ,7), Inert ial ( : ,8), Inert ial ( : ,9)]; 

disp(op5) 

else 

end 

answer* input( 'Do you want Ux and Uh? 1*YES '); 
if answer **1, 

dispC pt Ux Uh') 

op4*[flircraft( : ,7), Inert ial (:>5), Inert ial ( : ,6)]; 

disp(op4) 

else 

end 
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fund i on [ydot ]*ypr i ies(y, dUxdx, dUxdh, dllhdx, dUhdh, theta, 

T > S > n,Q > CL,CD 4 alpha,alphanax) 

X This function calculates the y’ solutions to the 

X coupled systei equations: 

X y(l ) '=xdot=U*cos(ga#aa)+Ux 

X y(2) ‘=hdot=U*sin(gai»aa)+Uh 

X y(3) ,= Uxdot=dllx/dx*xdot*dllx/dh*hdot 

X y(0 '=Uhdot=dUh/dx*xdot*dUh/dh*hdot 

X y(5) '=Udot=(T*cos(alpha))/H-D/ri-g*sin(gaiaa)-Uxdot*cos(ga*aa) 

X -Uhdot*sin(ga«aa) 

X y(6) ,e T*sin(alpha)/(H*U)+L/(f1*U)-g*cos(ga»aa)/U 

X +Hxdot*sin(gaB>aa)/U-Uhdot*cos(ga»aa)/U 



ydot(1)«y(5)*cos(y(6))+y(3)j 
ydot(2) e y(5)*sin(y(6))+y(4)j 
ydot (3)=dUxdx*ydot ( 1 )+dllxdh*ydot (2) ; 
ydot (1)«dUhdx*ydot ( 1 )+dUhdh*ydot (2) ; 
if alpha<alphatax, 

ydot (5)»T/n*cos(theta-y(6))-Q*S*CD/n-32 , 1 74*s I n(y(6) ) 
-ydot(3)*cos(y(6))-ydot(4)*sin(y(6))j 
ydot(6)“T/(n*y(5))*sin(theta-y(6))+Q*S*CL/(n*y(5)) 

-32. 174/y(5)*cos(y(6))+ydot(3)*sin(y(6))/y(5)- 

ydot(4)/y(5)*cos(y(6)); 

else 

ydot (5)-T/n*cos(a I pha»ax)-l}*S*CD/n-32 . 1 74*s i n(y (6) ) 
-ydot(3)*cos(y(6))-ydot(D* sin(y(6)); 
ydot(6)-T/(n*y(5))*sln(alpha«ax)+Q*S*CL/(n*y(5))-32.171/y(5) 

*cos(y(6))+ydot(3)* sln(y(6))/y(5)-ydot(1)/y(5)*cos(y(6)); 
end 



147 



B. AIRCRAFT PERFORMANCE PARAMETER FUNCTIONS 



fund lon[S,K,CDo,Tnax, 
=P3Ltfipp() 

XP3 at 89, 500 lbs, app 
S=1 300; 

K«0. 05041 ; 

CDo=0.0567; 

Teax®33400; 

Ut=69500; 

CLo a 0.800; 
CLalpha a 5.73; 
alpha»ax s 0. 244; 
q=0 . 0873; 

U1=236j 

T=5000; 

theta a -0.013; 



lit , CLo, CLa I pha, a I phaiax, q, U1 , T , theta ] 
flaps, gear doen 

Xplane for* surface area in ft A 2 
^coefficient for CO calculation 
Xcoefficient of drag 
Xeax thrust auailable in ft-lbs 
Xaircraft eeight in lbs 
Xzero flOfl lift coefficient 
XI i ft curve slope 

Xstall buffet alpha e/ approach flaps 
Xnoeinal pitch rate of 5deg/sec 
Xapproach equivalent airspeed in ft/s 
Xinitial thrust in Ibf 
Xinitial eater line deck angle in rad 



funct lon[S,K,C0o,Teax,Ut,CLo, 
-P3Hvflpp() 

XP3 at 114,000lbs, app flaps, 
S-1300; 

K-0. 05041; 

CDo-0.0567; 

Teax*33400; 

Ut-1 14500; 

CLo«0.800; 

CLalpha*5.73; 

alpha«ax*0.244; 

q*0.0873; 

U1-262; 

T-8000; 

theta»-0.013; 



CLa I pha, a I phaiax, q, U1 ,T, theta] 
gear doen 

Xplane for* surface area In ft A 2 
Xcoeff Sclent for CD calculation 
Xcoefflcient of drag 
Xeax thrust available in ft-lbs 
Xalrcraft eelght In lbs 
Xzero flOfl lift coefficient 
XI I ft curve slope 

Xstall buffet alpha «/ approach flaps 
Xnoeinal pitch rate of 5deg/sec 
Xapproach equivalent airspeed In ft/s 
Xinitial thrust in Ibf 
Xinitial eater line deck angle In rad 
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fund ion[S,K,CDo,Tiax,Ut ,CLo,CLalpha,alphaiax,q,U1 , T, theta] 
=L101 IflppO 

XL1011 at 362,000lbs, flpp flaps, gear down 



S=1578; 

K=0.059; 

CDo»0. 1 08 ; 

Tiax»126000; 

Ut-362000; 

CLo“0 . 532; 

CLalpha»1.96; 

aiphaiax»0.311; 


Xplane fori surface area in ft A 2 
Xcoefficient for CD calculation 
^coefficient of drag 
Xiax thrust aoailable in ft- 1 bs 
Xaircraft eeight in lbs 
Xzero flOR lift coefficient 
Xi i ft curue si ope 

Xstaii earning alpha e/ approach flaps 

X(16deg) 


q=0.0873; 

U1=227; 


Xnoninal pitch rate of 5deg/sec 
Xapproach equivalent airspeed in ft/s 

X(136kts) 


T*11500; 
theta«0. 118; 


Xinitial thrust in Ibf 
X Initial eater line deck angle in rad 



fund ion[S,K,C0o,Tiax,Ut,CLo,CLalpha,alphaiax,q,U1 ,T, theta] 
■TlIRppO 

XT11 at 6260 lbs, app flaps, gear do«n 



S=210; 

K=0.0503; 

CDo®0. 120; 
Tiax=3023; 

Ut **8280; 

CLo*0.587; 

CLaipha*6.2t; 

aiphaiax*0.211; 

q»0 . 0873; 

U1-203; 

T*1250; 

theta*-0.02; 


Xpiane fori surface area in ft *2 
Xcoefficient for CD calculation 
Xcoefficient of drag 
Xiax thrust available in ft-ibs 
Xaircraft eeight in lbs 
Xzero BOR lift coefficient 
XI i ft curve slope 

Xstall buffet alpha e/ approach flaps 
Xnoiinal pitch rate of 5deg/sec 
Xapproach equivalent airspeed in ft/s 
Xinitial thrust in Ibf 
Xinitial eater line deck angle in rad 
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funct ion[S,K,CDo,T»ax,lJt ,CLo,CLalpha,alpha»ax,q,U1 ,T, theta] 
=P3Eout() 

IP3 at 89,500 1 bs, app flaps, gear down, one engine out. 



S-1300; 

K-0. 05011; 

CDo-0.0630; 

Tmax=33400*(3/4); 

Ut -89500; 

CLo«0.800; 

CLalpha=5.73; 

alphawax=0.2‘H; 

q=0.0873; 

01=236; 

T=5000; 

theta*-0.013; 


Ip lane fore surface area In ft A 2 
Icoefflclent for CD calculation 
{coefficient of drag (inc for eng. out) 
Iwax thrust available in ft- lbs 
{aircraft weight in lbs 
{zero flOfl lift coefficient 
II i ft curve slope 

{stall buffet alpha w/ approach flaps 
Inowinal pitch rate of 5deg/sec 
{approach equivalent airspeed in ft/s 
{initial thrust in Ibf 
{Initial water line deck angle in rad 
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C. FLIGHT PATH CONTROL AND ESCAPE MANEUVER FUNCTIONS 



funct ionfthetoout ,Tout ]=GPcontro I (q, del tot , ROC, T#ox, theta, T, 

a I pha , a I phawax, U1 , hetart , xetart , y ) 
XThis control function tries to naintain a 3deg glide path. 
Tout'T; 

thetaout'theta; 

wantedh*hstart-(abs(xstart-y ( 1 ) )*0 . 0521) ; 



Xdeslred altitude for 3deg GP 
%a I t difference (neg when low) 
Xairspeed diff (neg when slow) 



hdi ff B y(2)-wantedh; 

Udi f f=y(5)-U1 ; 

X theta input: 
if alpha<0, 

thetaout'theta-a I pha ; 
el 8e If aipha<aipha«ax, 
if R0D0, 

thetaout B theta-q*de I tat ; 
el se if flOC<-25 . Oj, 

thetaout B theta+q*de I tat ; 
elseif hdi ff>100, 

thetaout B theta-q/2*de I tat ; 
eiseif hdi ff<-100, 

thetaout«theta+q/2*de I tat ; 
end 

If thetaout>0.524, 
thetaout *0.524; 
elseif thetaout<-0. 175, 
thetaout B -0. 175; 
end 
else 

thetaout»theta-(a I pho-a I phawax) ; 
end 

XThrust Input: 

If Udl f f >6.5, 

Tout*T-0 . 20*T«ax*de I tat ; 
if Tout<0, 

Tout»0; 

else 



X(1)check for negative ROR 

X(2)check for wax available ROR 
X(3)check for ROC > Ofpu 

Xif so pitch dwn q*t lie 
X(4)check for ROC < 1500fp» 

XI f so pitch up q*t ine 
X(5)check for alt >100 diff 

XI f so pitch dwn 1/2*q*tiwe 
X(6)check for alt <100 diff 

Xif so pitch up 1/2*q*ti»e 

Xnow check for theta linit of 
X -lOdeg to +30 deg 



Xcheck speed diff 
Xif diff > 5kts above app speed 
Xreduce thrust by 20X of Tiax 
Xbut not less than 0. 
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end 



elseif Udiff<-8.5, 

Tout“T+0 . 25*T«ax*de I tat ; 
if Tout>Tsax, 



XI f diff < 5kts beloi app speed 
Xinc thrust by 25X of Tnax 
Xbut not tore than Tnax. 



Tout=T»ax; 



else 

end 



end 



funct i on[ thetaout, Tout ]=f1flcontro I (q,deltat,Tiax,theta,T,alpha, 



XThis control function initiates a noma I sissed approach. 

X flax poser is added and theta increased to lOdeg. flax flOfl 
Xis applied shen required. 

Tout“T; 

thetaout “theta; 

X theta input: 

if alpha<alpha«ax, X(1)check for tax available flOfl 

if thetaout>“0. 175, X(2)check for lOdeg deck angle 

thetaout “0. 175; 

else 

thetaout a theta+q/2*deltat; X(3)lf theta not lOdeg, pitch to 

lOdeg 

end 

else 

t het aout “t het a- (a I pha-a I phanax ) ; 
end 

XThrust input: 
if Tout<Tiax, 

Tout-T+Tsax/2*de I tat ; 
if Tout>Tsax, 

TouWsax; 

end 

else 

end 



alphanax,y) 



152 



fund i ontthetaout , Tout , checkl >CflSesc(q, de I tat , Tmax, theta, 1 , 

a I pha, a I phamax, U1 ,y, check 1 ) 

Xflircraft escape maneuver consisting of constant airspeed. 
thetaout*thetaj 
if check1"0, 
if y(5XU1, 
thetaout s 0; 
else 

thetaout»q*deltat; 
checkM ; 
end 
else 

if y(5)>-(U1+8.4), 
thetaout B theta+q*deltat j 
elseif y(5X-(U1-B.-f) 
thetaout»theta-q*deltat ; 
end 
end 

if alpha>alphamax, 
thet aout B theta-(a I pha-a I phamax) ; 
end 

if T<Tmax, 

Tout“T+Tmax*0. 2*de I tat ; 
else 

Tout*Tmax; 

end 
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funct ion[thetaout , Tout , check 1 ]=CRLTesc(q,deltat ,T«ax, theta, T, 

alpha,alpha»ax,y,check1 ,ydot) 

XRircraft control using constant altitude escape naneuuer. 

thetaout*theta; 

if check 1»*0; 

check1°y(2); 

end 

if y(2)<check1-20, 
i f ydot(2)<«0, 
thetaout*theta+q*deltat ; 
end 

elseif y(2)>check1+20, 
if ydot(2)>*0, 
thetaout«theta-q*deltat ; 
end 
end 

if alpha>alphaaax, 
t hetaout *thet a-(a I pha-a I phaaax) ; 
end 

if T<T«ax, 

TouW+T*ax*0,2*deltat ; 
else 

Tout*T«ax; 

end 
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funct ion[thet aout, Tout ] a CTHesc(q, del tat, Tnax, theta, T, 

alpha,alpha«ax,thetaalft) 
{Constant theta escape naneuverv 
if theta<thetaaift, 
t het aout a t het a+q*de 1 1 at ; 
else 

thetaout a thetaai«; 

end 

if alpha>alpha»ax, 
t het aout =t het a-(a I pha-a I phanax) ; 
end 

if KTftax, 

T out *T +T «ax*0 . 2*de 1 1 at ; 
else 

TouWnax; 

end 



funct ion[thetaout,Tout] a CflOflesc(q,deltat,Tftax,theta,T, 

a I pha, a I phanax, a I phaa in) 

{Constant alpha escape naneuuer. 
if alpha<alphaai», 
t het aout a t het a+(a I phaa i n-a I pha ) ; 
if thetaout>(theta+q*deltat), 
t het aout »t het a+q*de 1 1 at ; 
end 

else 

t het aout a t het a- ( a I pha-a I phaa in); 
if thetaout<(theta-q*deltat), 
thetaout a theta-q*del tat ; 
end, end 

if alpha>alphanax, 
t het aout a t het a- ( a I pha-a I phanax ) ; 
end 

if KTftax, 

Tout=T+Tftax*0.2*deltat ; 
else, Tout a Tnax; end 
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funct ion[thetaout, Tout ]=nflXesc(q, del tat, T«ax, theta, T, alpha, a I phaai») 
Xriaxinui alpha escape ■aneuuer. 
if alpha<alphaai«, 
thetaout B theta+(alphaai»i-alpha) ; 
if thetaout>(theta+q*deltat), 
thetaout B theta+q*de I tat ; 
end 

else 

t het aout «t het a- ( a I pha-a I phaa i a ) ; 

If thetaout<(theta-q*deltat), 
thetaout B theta-q*deltat ; 
end 
end 

if T<T«ax, 

Tout B T+T«ax*0.2*deltat ; 
else 

Tout«T«axj 

end 
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funct ion[thetaout , Tout , check lout I'DRLthetaCFnax, theta, T,U1 ,y,check1 ) 
% This control function tries to enulate DHL 191 



Hast 18 sec of flight, 
if check1««0, ^starting point 

check 1=32; 
end 



X theta input : (theta, Xpoeer) 
DRL191-C0. 069813170071; 
0.0628318530 71; 0.051105206871; 

72; 0.017123889873; 0 
71; 0.022689280271; 0 
69; 0.051105206869; 0 
68; 0.123918376867; 0 
67; 0.185001900767; 0 
67; 0.200712863969; 0 
72; 0.253072711585; 0 
92; 0.271016692592; 0 
92; 0.237361778292; 0 
87; 0.253072711587; 0 
93; 0.062631853098; 0 
99; -0.005235987100; 
100; -0.130899693100; 
100; 0.0925021503100; 
100; 0.0511052068100; 
thetaout=DRL191 (check 1+1 , 1 ); 
Tout=DAL191(check1+1,2)/100*T«ax; 



0.0171238898 

0.0311159265 

0.0511052068 

0.0785398163 

0.1623156201 

0.2007128639 

0.2231021112 

0.2710166925 

0.2600510585 

0.2530727115 

0.1553313031 

0.0226892802 

-0.111862327 

0.0698131700 

-0.005235987 



.038397213573 
.038397213569 
.051105206870 
.132615023167 
.191986217767 
.200712863971 
.267035375589 
.267035375595 
.237361778292 
.223102111292 
.031115926596; 
-0.068067810100; 

-0.020913951100; 

0.0311159265100; 

0.0311159265100]; 



check1out**check1+1 ; 
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APPENDIX D 



TAKEOFF ALGORITHMS 



Contained in this appendix are the MatLab® programs utilized in 
calculating the aircraft response given a windshear and a particular 
escape maneuver. They are listed in three sections. The first section 
contains the driver routine. The second and third sections contain 
the specific aircraft parameter functions and escape maneuver 
functions respectively. The windshear function Is listed In Appendix 
A. The differential equation solving routine is listed in Appendix C. 

A. DRIVER 

g************************»**** 

X Takeoff Profile Mode I * 

g***************************** 

X This prograa ell I calculate the Inertial flight path angle for 
X a g i yen aircraft. It utilizes the function 'eicroburst' ehich 
X aodels the eind shear encountered during the DFU Delta Reel dent. 

X This prograa uses energy height theory to calculate inertial space 
position. 

clear 

clc 

foreat bank 
foreat coapact 

X 

XINITIALIZRTION AND EHU I ROHMEHTAL CONSTANTS 

X The fo I loving Initialize constants and variables for the aodei. 
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Xfl I RCRflFT CONSTANTS 

X The fol loving constants are for the particular aircraft being 
studied. 

dispC TAKEOFF MODEL') 

dispC ’) 

ansver-0.0; 

disp( 'Type aircraft : ' ) 
dispC 1. P3 at 90,0001b') 

dispC 2. P3 at 120,0001b') 

dispC 3. P3 at 135,0001b') 

dispC 4. LI 0 1 1 ' ) 

dispC 5. T44') 

dispC 6. P3 at 120,0001b, Ur-140') 

ansver- i nput ( ' Enter desired aircraft by nuvber: '); 

if ansver— 1, 

[S,K,CDo,Ut ,CLo,CLalpha,alphavax,q,U1 ,U2,T,thetacl ivb]-P3LtToff ; 
theta-0.087; Xinit ial izat ion 

el seif ansver— 2, 

[S,K,CDo,Ut,CLo,CLalpha,alpha«ax,q,U1 ,U2,T,thetacl l«b]=P3HuToff ; 
theta-0.130; Xinit ial izat ion 

el seif ansver— 3, 

[S,K,CDo,Ut,CLo,CLalpha,alphavax,q,U1,U2,T,thetacl lvb]-P3HuHuToff ; 
theta-0.133; Xlnlt Ial izat Ion 

el seif ansver— 4, 

[S,K,CDo,Ut,CLo,CLalpha,alphavax,q,U1 ,U2,T,thetacl I vb]-L1 01 1Toff;XL1 
011 

theta-0.2; Xinit ial izat ion 

el seif ansver— 5, 

[S,K,CDo,Ut,CLo,CLalpha,alphavax,q,U1 ,U2,T,thetacl ivb]-T44Toff; 

XT44 

theta-0.200; Xlnlt ial izat ion 

el seif ansver— 6, 

[S,K,CDo,Ut,CLo,CLalpha,alphavax,q,U1 ,U2,T,thetacl ivb]-aodP3HuTo; 
theta-0.130; Xinit Ial izat ion 

end 



159 



XS=piane fora surface area in ft ^2 
XK*coeff icient for CD calculation 
XCDo=coeff icient of drag 
XUt*alrcraft weight in ibs 
n*Ut/32.171; X Mass in slugs 
XCLo=zero flOfl i i ft coefficient 
XCLalpha=l i ft curve slope 
Xalphaaax=stai I buffet aipha w/ approach flaps 
Xq«noainai pitch rate of 5deg/sec 
XU Wot ate speed in ft/s 
XU2"takeoff safety airspeed In ft/s 
XT*thrust In Ibf 

X 



X***progra«ed Initial values*************************************** 
iters"65; Xnuaber of iterations 

del tat *0.5; Xtiae increaent in sec 

xstart*-1200; Xstarting distance froa 

riB center 

hstart=2.0; Xstarts at rotate 

j****************************************************************** 



disp(' initial input Section') 
dispC ') 

dispCThe foi lowing values are preprograaaed: 1 ) 

dispC a. The flight path Is based on wax power.') 

dlsp(' b. The starting distance (distance froa H6 center) is:') 

disp(xstart) 

dispC c. The flB center is this froa the end of the runway:’) 

disp(abs(xstart)) 

dispC d. The tiae step in sec: 1 ) 

disp(deltat) 

dispC e. The nuaber of iterations is:’) 
dispC Iters) 
dispC ') 



TU*T/Ut; Xthrust/weight ratio 

UL*Ut/S; Xwing loading 

dispCThe following values were supplied/calculated froa aircraft 
constants' ) 
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dispC a. Takeoff Safety Speed:’) 
disp(U2) 

dispC b. Thrust to weight is;') 
disp(TU) 

dispC c. U i ng loading is:') 

disp(UL) 

dispC ') 

aaneuoer=0.0; 

dispCSelect the escape aaneuoer:’) 
dispC 1. Takeoff flight path; HO sindshear') 
dispC 2. Takeoff flight path; UiTH windshear’) 
dispC 3. Escape with constant airspeed 1 ) 
dispC 4. Escape with constant altitude') 
dispC 5. Escape with constant theta') 
dispC 6. Escape with constant alpha') 
aaneuoer- inputC Enter selection nuaber: '); 

If aaneuuer--5, 

dispCUhat ualue for escape theta?') 
thetaala-InputC Enter in radians: '); 
else if aaneuoer— 6, 
dispCUhat ualue for escape alpha?') 
alphaaia-input( ' Enter In radians: '); 



end 




xl-xstart; 


Unit iai izat ion 


hl-hstart; 


{Initialization 


Ug»Ul ; 


{ground speed Initialization In ft/s 


gaaai *0.001 ; 
rad (Ideg) 


{initial inertial flight path angle 


row-0. 002377; 


{air density In slugs/ft A 3 t S.L. 


trigger-0.0; 


{utility check nuaber 


check! -0.0; 


{utility check nuaber 


check2-0.0; 

X 

XRLGORiTHn 


{utility check nuaber 
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X The actual algorithm is broken into several smaller divisions and 
functions: 

X starting values - as it implies 
X loop - beginning of iterating process 
X mindshear inputs - uses the function ’microburst' to give 
X point x and z axis mind speeds 

X predictor - first guess at the solution of the governing DEs 

X correcter - corrected solutions to the governing DEs and 

X final point values. 

X output file management - as it implies 
X aircraft control - uses multiple control functions for 
X control inputs to theta, thrust, and flOR 

X final output - as it implies 



X 

Xstart ing values: 
pt-1 ; 


Xcounter 


Q"0.5*rof*U1 A 2j 


Xdynamic pressure 


CL-Ut/(q*S)j 


Xcoeff icient of 1 1 ft 


CD«CDo+K*CL A 2; 


Xcoefflcient of drag 


a 1 pha*(CL-CLo)/CLa 1 pha; 


XROR required for 1g flight 


gamaa B theta-aiphaj 


Xairmass flight path angle 


Es-lir2/(2*32.1?4)+h1j 


Xspecific energy in ft 


Esdot B -12.32; 


Xtlme rate change In specific energy 


y(1) B x1; 


Xy(1 ) B x posit ion in ft 


y(2)-h1; 


Xy(2) B alrcraft altitude in ft 


y(5)-U1; 


Xy(5) a aircraft airspeed in ft/s 


y(6)-gamaa; 


Xy(6) B aircraft airmass flight path in 


rad 

Aircraft (1 , : )-[x1 ,h1 ,U1 , 


-12. 36, theta, a 1 pha, 1 ,T]j 


Inert laid, : MUg,0,gamaa,gamal ,0,0, 1 ,Es,Esdot]j 


% 

Xloop: 

clc 


X start ing output files 


•hile theta*’=999, 


Xloop control 


pt-pt+l ; 

X 


Xpresent pt in space 
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%• i nd shear inputs: 

x2-Ug*deltat+x1 ; lapproximate horizontal displacement 

h2-tan(gamai )*(x2-xl )+hl ; lapproximate vertical displacement 

if maneuver~-1, 

[Uxl 1 ,Uh1 I ]-microburst(h1 ,x1 ); IThese four calls to the function 
[11x21 ,Uh21 ]=microburst(h2,x1 ); Imicroburst finds the point mind 
[Ux12,Uh12]-microburst(h1 ,x2); Ishear for determination of the 
[Ux22,Uh22]=microburst(h2,x2); laverage block mind shear, 
else 

Uxl 1*0 ; Ux2 1 -0 ; Ux 1 2-0 ; Ux22«0 ; Uh 1 1 -0 ; Uh2 1 -0 ; Uh 12-0; Uh22«0 ; 

XThe above Is used in conjunction 
lother calls mhen no mind shear is 
Ides i red. 
end 

UlndX(pt)-Ux1 1 ; XU I nd shear matrices 

UindH(pt)«Uh11; 

dUxdx-((Ux12-Ux1 1 )+(Ux22-Ux21 ))/(2*(x2-x1 ));IThese formulas 
dUxdh«( (Ux22-Ux1 2)+(Ux21 -Uxl 1 ) )/(2*(h2-h1 ) ) ; lea I cu I at e the 
dUhdx*((Uh12-Uh1 1 )+(Uh22-Uh21 ))/(2*(x2-x1 ));ldl fferent ial change 
dUhdh*((Uh22-Uh12)+(Uh21-Uh1 1 ))/(2*(h2-h1 ));lin mind shear. 

y(3)-Ux11; ly(3)-mind shear in the x direction 

ft/s 

y(1)-Uh11; ly(4)-mlnd shear in the z direction 

ft/s 

X 

Ipredlctor and Intermediate values 

ydot-ypr I mes(y, dUxdx, dUxdh, dUhdx, dUhdh, theta, T, S, U, Q, CL, CD, a I pha, a I p 
hamax); 

leal I to the DE function 



for n-1 :6, 

yl (n)-y(n)+deltat*ydot(n); lEuler first step fmd predictor 
yhal f(n)-y(n)+deltat/2*ydot(n); lEuler half step fmd predictor 
end 
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aipha=theta-y(6); 
CL=CLo+a I pha*CLa I pha; 
CD«CDo+K*CL A 2; 

0*0 . 5*roe*yha I f (5) A 2; 

X 



XflOfi 

Xcoef f icient of lift 
Xcoefficient of drag 
Xdynaaic pressure 



Xcorrector and nee value fonuiation 
ydot=ypri*es(yhal f,dllxdx,dUxdh,dUhdx,dUhdh, theta, T,S,f1,Q,CL,CD,ai pha 
, alpha* ax); 



for n*1 :6, 

y2(n)=yha I f (n)+de I tat/2*ydot (n) ; XR i chardson extrapo I at i on 
y(n)*2*y2(n)-y1 (n); Xcorrecter scheie. 

end 



x1-y(1 ); 
h1-y(2); 

Ug-y(5)*cos(y(6))+y(3); 

R0C-y(5)*sin(y(6))+y(4); 

ga*ai-atan(ROC/Ug); 

Es*y(5) A 2/64.348+y(2); 



Xne« x posit in ft 
Xnee altitude in ft RGL 

Xground speed in ft /sec 
XRate of C i I *b in ft/sec 
Xne« ainass FP angle in rad 
Xne* specific energy 



Esdot*y(5)*ydot(5)/32. 174+ydot(2); Xnee trc of specific energy 



aipha*theta-y(6); 
CL*CLo+a I pha*CLa i pha ; 
CD*CDo+K*CL A 2; 
Q-0.5*roe*y(5) A 2; 

X 



Xnee ROR 

Xnee coefficient of lift 
Xnee coefficient of drag 
Xnee dynaeic pressure 



Xoutput file eanageaent: 

Rircraft(pt, : )*[y(1 ),y(2),y(5),ydot(5),theta,alpha,pt ,T] ; 

Inert iai(pt, : )-[Ug,R0C,y(6),ga*ai ,y(3),y(4),pt,Es,Esdot]; 

ho»e 

dlspC pt x h U 

airspeed accei ') 

opMptjfllrcraftfpt, 1 ), Aircraft (pt , 2), Aircraft (pt, 3), Aircraft (pt, 4)] 

« 

J 
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disp(opl ) 

dispC pt theta alpha ganaa 

ganai ' ) 

op2*[pt, Aircraft (pt, 5), Aircraft (pt ,6), Inert ial (pt ,3), Inert ia I (pt , 

1 ) 1 ; 

disp(op2) 

dispC pt Ug AOC Thrust') 

op3=[pt , Inert ial (pt , 1 ), Inert ia I (pt , 2 ), Aircraft (pt ,8)]; 
disp(op3) 

dispC pt Es Esdot') 

op5-[pt, Inert ia I (pt ,8), Inert ial (pt , 9)]; 

disp(op5) 

dispC pt Ux Uh' ) 

op4=[pt , Inert ial (pt ,5), Inert ial (pt ,6)]; 
disp(op4) 

% 



Jaircraft control: 

Jfunctions called for theta, thrust 
If y(2)>50, land alpha control, 

if ROC <-0, 
trigger-1; 
elseif y(5)<U2, 
trigger- 1 ; 
end 

elseif y(5)<U1, 
trigger-1 ; 
end 

if trigger—0, 
thet a-theta+q*0 . 0*de I tat ; 
i f theta>thetacl iab, 
theta-thetacl i»b; 
end 
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else 

If aaneuuer“=3, 

[ theta, check 1 ] iE ToffCRS(q,deltat,theta,alpha,alphaaax,U2,y, check! ) 
el seif aaneuuer*M, 

[theta, checkl ]=Tof fCRLT (q, de 1 1 at , theta, a I pha, a I phaaax, y, check 1 , ydot ) 
el seif aaneuver s!B 5, 

[ t het a] =To f f CTH ( q , de 1 1 at , t het a , a I pha , a I phaaax , t het aa i a ) 
el seif aaneuuer==6, 

[theta]*ToffCROfi(q,deltat,theta,alpha,alpha»ax,alphaaia) 
end 
end 

If pt“*lters, 
theta 1B 999j 
elseif pt >5, 
if y(2)<0.0, 
theta*999; 
end 
end 

end 

% 

Xoutput files: 

ansaer- Input (’Do you aant a flight path plot? 1«VES ’); 
if ansier”1, 

plot(Rircraft(:,1),Rlncraft(:,2)) 

pause 

ansver B 0.0; 

end 

ans«er B lnput( ’Do you aant x, h, U, and airspeed accel? 1*VES ’); 
if ansaer”1, 

dlspC pt x h U 

airspeed accel ') 

op1"[Rlrcraft( : ,?), Aircraft ( : , 1 ),Rlrcraft( : ,2),fllrcraft( : ,3),fllrcraf 
t ( : , 4 ) ] ; 



Xnuaber of Iterations 



Xahi le aigoritha loop 
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disp(opl ) 

else 

end 

anseer*input( ’Do you want theta, alpha, gaeaa, and ganai? 1 b VES '); 
if answer==l , 

dispC pt theta alpha ganaa 

gaia i ' ) 

op2»[flircraft( : ,7),flircraft( : ,5),flircraft( : ,6), Inert ial ( : ,3), Inertia 
!(:,«)]; 
disp(op2) 
else 
end 

anseer= input ('Do you eant Ug, ROC, and Thrust? I=VES '); 
if ans»er==l, 

disp(' pt Ug ROC Thrust') 

op3 s [fiircraft( : ,7), Inert ial ( : , 1 ), Inert ial ( : ,2),flircraft( : ,8)); 

disp(op3) 

else 

end 

anseer- input ('Do you eant Es and Esdot? 1 “VES '); 
if anseer«“1, 

dispC pt Es Esdot') 

op5*[flircraft( j , 7), inert iai ( : ,8), Inert ial ( : ,9)] ; 

disp(op5) 

else 

end 

anseer*input( 'Do you eant Ux and Uh? 1 “VES ’); 
if anseer»“1, 

diep(' pt Ux Uh') 

op4*[flircraft( : ,7), Inert ial ( : ,5), Inert ial ( : ,6)]; 

disp(op4) 

else 

end 
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B. AIRCRAFT PERFORMANCE PARAMETER FUNCTIONS 



funct ion[S,K,C0o,Ut,CLo,CLalpha,alphaiax,q,01 , 02, T, thetacl lab] 



«P3LtToff() 

XP3 at 90,000168, app flaps, 
S=1300; 

K=0. 05041 ; 

C0o*0.0551 ; 

Ut =90000; 

CLo-1 .000; 

CLalpha*5.73; 

alphaiax=0.244; 

q=0.0873; 

U1-204; 

02 = 220 ; 

T=33400; 

thetacl iib=0. 175; 



gear up 

{plane fori surface area in ft A 2 
{coefficient for CD calculation 
{coefficient of drag 
{aircraft lelght In lbs 
{zero ROfl lift coefficient 
{lift curue elope 

{stall buffet alpha •/ approach flape 
{noilnal pitch rate of 5deg/eec 
{Liftoff epeed ( 1 21 kte) 

{takeoff safety alrepeed In ft/s 
{thrust in Ibf 

{eater line deck angle In rad 



functlon[S,K,C0o,Ut 
«P3HvToff() 
{P3 at 120, 000 1 be, 
S-1300; 

K-0.05041 ; 
CDo-0.0551; 

Ut=1 20000; 

CLo= 1.000; 
CLalpha*5,73; 
alphaiax=0.244; 
q=0,0873; 

01=214; 

02=227; 

T*33400; 

thetacl lib=0. 175; 



, CLo , CLa I pha , a I phaiax , q , 0 1 , 02 , T , t het ac 1 1 lb ] 

app flaps, gear up 

{plane fori surface area in ft A 2 
{coefficient for CO calculation 
{coefficient of drag 
{aircraft eeight in lbs 
{zero R0R lift coefficient 
{lift curue elope 

{stall buffet alpha •/ approach flape 
{noilnal pitch rate of 5deg/sec 
{Liftoff epeed (12?kte) 

{takeoff safety alrepeed in ft/s 
{thrust In Ibf 

{•ater line deck angle In rad 
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fund I on[S, K, COo^ Ut ,Clo 
«P3HuHuToff() 
XP3 at 135,000lbs, app 
S-1300; 

K=0. 05041 j 
CDo=0.0551 ; 

Ut « 1 35000; 

CLo- 1,000; 

CLalpha=5.73; 

alphaiax=0.244; 

q*0.0873; 

01=229:, 

02=239; 

T=33100; 

thetacl iib=0. 1745; 



,CLalpha,alphaiax,q,01 , 02, T, thetacl itb] 
flaps, gear up 

Xplane fori surface area in ft ^2 
Xcoefficient for CD calculation 
Xcoefficient of drag 
Xaircraft weight In lbs 
Xzero flOfl lift coefficient 
XI i ft curue slope 

Xstall buffet alpha w/ approach flaps 
Xnoiinal pitch rate of 5deg/sec 
XL i ft of f speed (J36kts). 

Xtakeoff safety airspeed in ft/s 
Xthrust in Ibf 

Xiater line deck angle in rad 



funct ion[S,K,C0o,Ut,CLo 
=iodP3HuTo() 

XP3 at 120, 000 lbs, app 
S-1300; 

K-0. 05041; 

CDo-0.0551; 

Ut-120000; 

CLo-1.000; 

CLalpha=5.73; 

alphaiax=0.244; 

q*0.0673; 

01=242; 

02-254; 

T=33400; 

thetacl lib-0. 175; 



,CLalpha,alphaiax,q,01 ,02, T, thetacl iib] 

flaps, gear up. Increased rotate speed. 
Xplane fori surface area In ft A 2 
Xcoefficient for CO calculation 
Xcoefficient of drag 
Xaircraft weight in lbs 
Xzero flOfl lift coefficient 
XI I ft curue slope 

Xstall buffet alpha w/ approach flaps 
Xnoiinal pitch rate of 5deg/sec 
XLIftoff speed (145kts) 

Xtakeoff safety airspeed In ft/s 
Xthrust in Ibf 

Xwater line deck angle In rad 
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funct ion[S,K,CDo,Ut,CLo,CLaipha,aiphaaax,q,U1 ,U2,T, thetaci iab] 
=L101 ITof f() 

XL 1 0 1 1 at 462,000ibs, flpp flaps, gear up 



S-4578; 

K=0.059; 

CDo*0.098; 

Ut=462000; 

CLo*0.532; 

CLalpha=4.96; 

alphaaax=0.3H; 


Xpiane fora surface area In ft A 2 
Xcoefficient for CD calculation 
Xcoefficient of drag 
Xaircraft eeight in lbs 
Xzero ROfl lift coefficient 
XI I ft curve slope 

Xstali earning alpha «/ approach flaps 

X(18deg) 


q=0.0873; 

U 1=238; 

U2*255; 

T= 126000; 
thetaci iab=0.209; 


Xnoainai pitch rate of 5deg/sec 
Xrotate speed in ft/sec (Mlkts) 
Xtakeoff safety airspeed ft/s ( 1 51 kts) 
X thrust in Ibf 
Xinitiai deck angle rad(12deg) 



funct ion[S,K,CDo,Ut,CLo,CLalpha,alphaaax,q,U1 ,U2,T, thetaci lab] 
•TMToffO 

XT44 at ?81?lb3, app flaps, gear doin 

S=210; Xpiane fora surface area In ft ^2 



K=0.040; 

CDo=0.010; 

Ut*7817; 

CLo-0.0524; 

CLalpha=6.24; 

aiphaaax=0.227; 

q=0.0873; 

U1-152; 

U2-180; 

T-3023; 

thetaci iab=0. 262; 


Xcoefficient for CD calculation 

Xcoefficient of drag 

Xaircraft eeight in ibs 

Xzero R0R lift coefficient 

Xi i ft curve slope 

Xstail buffet alpha e/fiaps up 

Xnoainai pitch rate of 5deg/sec 

Xilftoff airspeed In ft/s(90kts) 

Xtakeoff safety speed ft/s (107kts) 

Xinitiai thrust in Ibf 

Xinitiai deck angle in rad (15 deg) 
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C. FLIGHT PATH CONTROL AND ESCAPE MANEUVER FUNCTIONS 



funct lon[thetaout,check1 ] 

■Tof fCflS(q,deltat,theta,alpha,alpha*ax,U2,y,check1 ) 
Xfllrcraft escape Maneuver consltlng of constant airspeed, 
thetaout “theta; 

If check1“*0, 

If y(5)<U2, 

thetaout*theta-q*deltat; 

If thetaout<0.0, 
thetaout 3 0.0; 
end 



else 

thetaout»q*deltat; 

checkl-1; 

end 

else 



If y(5)>-(U2*8.1), 



t het aout -t het a+q*de It at ; 
elself y(5)<-(U2-8.4) 
thetaout»theta-q*deltat; 
end 
end 

If alpha>alpha«ax, 
thetaout>theta-(alpha-alphaaax); 
end 
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funct ion[thetaout,checkl ] 

■Tof fCflLT (q, de I tat , theta, a I pha, a I phaaax, y, checkl , ydot ) 
Xflircraft control using constant altitude escape naneuuer. 
thetaout *theta; 
if check1*=0; 

check 1 *y (2)+20 j 
end 

if y(2)<check1-IO, 
if ydot(2)<=0, 
thetaout*t heta+q*de I tat ; 
end 

elseif y(2)>check1+10, 

If ydot(2)>-0, 
thetaout"theta-q*deltat; 
end 
end 

if thetaout<0, 
thetaout'O; 
end 

if alpha>alphaiax, 
thetaout"theta-(a Ipha-a Iphasax) ; 
end 



funct ion[thetaout]"ToffCTH(q,deltat,theta,alpha,alpha«ax,thetaai«) 
^Constant theta escape Baneuuer. 
if theta<thetaaiB, 
thetaout»theta+q*deltat; 
else 

thetaout»thetaai«; 

end 

if alpha>alphaeax, 
thetaout*thet a-(a I pha-a I phaaax) ; 
end 
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fund ion[thetaout ]*ToffCROR(q,deltat , theta,alpha,alphaaax,alphaaia) 
IConstant alpha escape eaneuuer. 
if alpha<alphaaia, 
thetaout*theta+(alphaaia-alpha); 
if thetaout>(theta+q*deltat ) , 
thetaout=theta+q*deltat ; 
end 

else 

thetaout=theta-(a I pha-a Iphaa in); 
if thetaout<(theta-q*deltat), 
thetaout s theta-q*deltat ; 
end 
end 

if alpha>alphaaax, 
thetaout *t heta-(a I pha-a I phaaax) ; 
end 
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APPENDIX E 



ON-STATION ALGORITHMS 

Contained in this appendix are the MatLab® programs utilized in 
calculating the response of a P-3 encountering a windshear while 
operating close to the surface. They are listed in three sections. The 
first section contains the driver routine. The second and third 
sections contain the aircraft parameter function and flight path 
control function respectively. The windshear function is listed in 
Appendix A. The differential equation solving routine is listed in 
Appendix C. 

A. DRIVER 

g*****t*tt*t****t**t«*t*t***** 

X P3 Loiter Profile flodel * 

g***************************** 

clear 

clc 

foraat bank 
foreat coepact 

X 

XINITIRLIZRTION AND EHU I ROHHENTAL CONSTRNTS 

X The folloeing initialize conetante and variables for the eodel. 
IRIRCRRFT CONSTANTS 

X The folloeing conetante are for the particular aircraft being 
studied. 

diepC P3 ON-STRTION' ) 
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dispC ') 
check 1 “0.0; 

U=0.0; 

ansser =1 .0; 

check 1 s i nput ( ’ Enter on-station loiter alt itude(ft) : ’); 

U=inputC Enter loiter airspeed in knots: '); 

U“U*I.687; Xconuert knots to ft/sec 

ansser* Input ( 'Enter 0 for no sindshear: '); 

[S,K,CDo,Ut J CLo,CLalpha J alpha«ax J q l T i theta]“P3onsta(U) ; 



XS“plane fori surface area in ft A 2 
XK*coeff icient for CD calculation 
XCDo=coef f icient of drag 
XUt “aircraft ieight in lbs 
f1*Ut/32.l71; X Bass in slugs 
XCLo»zero ROR lift coefficient 
XCLalpha"! I ft curue slope 
XalphaBax“stal I buffet alpha s/ approach flaps 
Xq a noBinal pitch rate of 5deg/sec 
XU1 -rotate speed in ft/s 
XU2«takeoff safety airspeed In ft/s 
XT“thrust in Ibf 



g***pp 0 g raae( j initial ua I ues************* ************************** 



I****************************************************************** 

dispC Initial Input Section') 
dispC ') 

disp('The fol losing ualues are preprograBBed: ' ) 

dispC a. The flight path is based on constant poser. ') 

dispC b. The starting distance (distance froa HB center) is:') 

disp(xstart) 

dispC c. The tlae step in sec:') 



Xinit la I izat ion 



X 



iters“M5; 
del tat -0.25; 
xstart“-5000; 



Xnusber of iterations 
Xtlse incresent In sec 
Xstartlng distance fros 
flB center 



175 



disp(deltat) 

dispC d. The nuiiber of iterations is:') 
disp( iters) 
dispC ') 



TU-T/Ut; Xthrust/weight ratio 

UL»Ut/S; Xw i ng loading 

dispCThe following values were supplied/calculated from aircraft 
constants' ) 

dispC a. Loiter fiir Speed:') 
disp(U) 

dispC b. Thrust to weight is:') 
disp(TU) 



dispC c. U i ng loading is: 

disp(UL) 

dispC ') 

xl-xstart; 

hl-checkl ; 

Ul-U; 

Ug-U; 

gawai ■0.001 ; 
rad ( 1 deg) 
row-0.002377; 
trigger-0.0; 
check2-0.0; 



') 



Xinit iai izat ion 
X i n i t i a i i zat i on 
Xinit iai izat ion 

Xground speed initialization in ft/s 
Xinit iai inertial flight path angle in 

Xair density in siugs/ft A 3 • S.L. 
Xutiiity check nuaber 
Xut 1 1 ity check nuaber 



XRLGORiTHn 



X The actual aigoritha is broken into several saaiier divisions and 
funct ions: 

X starting values - as it iapiies 
X loop - beginning of iterating process 
X windshear inputs - uses the function ‘aicroburst* to give 
X point x and z axis wind speeds 

X predictor - first guess at the solution of the governing OEs 

X corrector - corrected solutions to the governing DEs and 

X final point values. 

X output file aanageaent - as it iapiies 
X aircraft control - uses auitipie control functions for 
X control Inputs to theta, thrust, and ROR 
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X final output - as it implies 



X 

{starting values: 
pt “ 1 ; 


{counter 


Q-0.5*roi*Ur2; 


{dynamic pressure 


CL-Ut/(Q*S); 


{coefficient of lift 


CD b CDo+K*CL / '2; 


{coefficient of drag 


a 1 pha=(CL-CLo)/CLa 1 pha; 


XAOA required for 1g flight 


gamaa s t het a-a 1 pha ; 


Xairmass flight path angle 


Es b U 1 A 2/(2*32. 171)+h1 ; 


{specific energy in ft 


tsdot=U; 


'{time rate change In speciTic energy 


y(1)=x1; 


Xy(1 )=x posit ion in ft 


u 

CM 

ZJ> 


{y(2)=aircraft altitude in ft 


y(5)»U1 ; 


Xy(5)*alrcraft airspeed in ft/s 


y(6)»gamaa; 

rad 


Xy(6) B aircraft airmass flight path in 


Aircraft (1 , : )»[x1 , hi ,U1 ,0, theta, alpha, 1 ,T]; 

Inert ial (1 , : )*[Ug,0,gamaa,gamai ,0,0, 1 « Es ^ Esdot ] ; 


X 

Xloop: 


X starting output flies 


clc 




mhile theta~»999, 
pt»pt+1; 

X 


{present pt In space 



{mind shear inputs: 

x2*Ug*deltat+x1 ; {approximate horizontal displacement 

h2"tan(gamal )*(x2-x1 )+h1 ; {approx i sate vert i cal displacement 



i f ansier"“0, 

[Uxl 1 ,Uh11 ]»mlcroburst(h1 ,x1); {These four calls to the function 
[Ux21 ,Uh21 ]*Bicroburst(h2 i x1 ); Xmicroburst finds the point Bind 
[Ux12,Uh12] B aicroburst(h1 ,x2); {shear for deteriinat ion of the 
[Ux22,Uh22]“microburst(h2 J x2) j {average block Bind shear, 
if pt<*8, {These next If logic commands 

Uxl 1*11x1 1*pt/8jUh1 1 “Uhl 1*pt/8; {control the entry and exit 
Ux12*Ux12*pt/8;Uh12*Uh12*pt/8; Xmindshear effects. 
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elseif y(1 )> s abs(xstant) < 
if check2==0, 
check2-pt; 
end 

if (8-pt+check2)>»0, 

Ux21=Ux2lM8-(pt-check2))/8;Uh2l»Uh21*(8-(pt-check2))/8; 

Ux22=Ux22*(8-(pt-check2))/8jUh22-Uh22*(8-(pt-check2))/8j 

else 

Ux 1 1 *0; Ux21 *0 ; Uxl 2*0; Ux22=0 ; Uhl 1 =0; Uh21 =0; Uhl 2=0; Uh22=0 ; 
end 
end 
else 

Ux 1 1 -0 ; Ux2 1 -0 ; Ux 1 2-0 ; Ux22-0 ; Uh 11 =0 ; Uh2 1 -0 ; Uh 1 2-0 ; Uh22»0 ; 

XThe above is used in conjunction 
Xother calls ehen no elnd shear Is 
Ides i red. 

end Xansier 

UindX(pt)-Uxl 1 ; XUInd shear latrices 

U I ndH(pt )— Uhl 1 ; 

dUxdx-((Ux!2-Ux1 1 )+(Ux22-Ux21 ))/(2*(x2-x1 ));XThese fonulas 
dUxdh-((Ux22-Ux12)+(Ux21-Uxl 1 ))/(2*(h2-h1 ));Xcalculate the 
dUhdx-( (Uhl 2-Uhl I )+(Uh22-Uh21 ) )/(2*(x2-x1 ) ) ; Xd I f f erent I a I change 
dUhdh-((Uh22-Uh12)+(Uh21-Uh1 1 ))/(2*(h2-h1 )) ;Xln elnd shear. 

y(3)-Ux11; Xy(3)-ilnd shear In the x direction 

ft/s 

y(1)-Uh11; Xy(1)»ilnd shear in the z direction 

ft/s 



Xpredictor and intermediate values 
ydot-ypriies(y,dUxdx,dUxdh,dUhdx,dUhdh,theta,T,S,n,QjCL,CD, 
alpha,alpha«ax); Xcal I to the DE function 

for n-1 :6j 

yl (n)-y(n)+deltat*ydot(n)j XEuler first step fid predictor 
yhal f(n)-y(n)+deltat/2*ydot(n); XEuler half step fid predictor 
end 
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alpha«theta-y(6); 
CL=CLo+a I pha*CLa I pha; 
CD«CDo+K*CL A 2; 
Q»0.5*roe*yhal f(5) A 2; 

{ 



{ROR 

{coefficient of lift 
{coefficient of drag 
{dynamic pressure 



{corrector and nee ualue fonulation 
ydot»ypriies(yhalf ,dUxdx,dUxdh,dHhdx,dUhdh, theta, T,S,M,Q> 

CL , CD, a I pha , a I phaaax ) ; 

for n=l :6, 

y2(n)*yhal f(n)+deltat/2*ydot(n);{Richardson extrapolation 
y(n)»2*y2(n)-y1 (n); {corrector scheae. 

end 



xl-y(l); 
hi »y(2) ; 

Ug»y(5)*cos(y(6))+y(3); 
R0C»y(5)*sin(y(6))+y(4); 
gaaa i »atan(ROC/Ug) ; 
Es«y(5) A 2/64.348+y(2); 



{nee x posit in ft 
{nee altitude in ft RCL 
{ground speed in ft/sec 
{Rate of Cl ieb in ft/sec 
{nee airaass FP angle in rad 
{nee specific energy 



Esdot=y(5)*ydot(5)/32. 1?4+ydot(2); {nee trc of specific energy 



alpha*theta-y(6); 
CL»CLo+a I pha*CLa i pha; 
CD*CDo+K*CL A 2; 
Q*0.5*roe*y(5) A 2; 

% 



{nee ROR 

{nee coefficient of lift 
{nee coefficient of drag 
{nee dynaaic pressure 



{output file aanageaent: 

RircraftCpt, : )»[y(1 ),y(2),y(5),ydot(5),theta,aipha,pt,T]; 

Inert ial(pt, ; )»[Ug,R0C,y(6),gaaal ,y(3),y(4),pt,Es,Esdot]; 
hoae 

disp(‘ pt x h 

airspeed accel ') 

op1»[pt,fllrcraft(pt , I ), Aircraft (pt,2), Aircraft (pt, 3), 

Aircraft (pt, 4)3 ; 
disptopl ) 

dispC pt theta alpha gaaaa 

gaaai ' ) 

op2»[pt,flircraft(pt,5),Rircraft(pt,6), inert ial (pt,3), 



U 
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Thrust ' ) 



InertiaKpt, 4)]; 
disp(op2) 

d!sp( 1 pt Ug ROC 

op3»[pt, Inert iai (pt, 1 ), Inert ial (pt, 2), Aircraft (pt,8)]j 
disp(op3) 

dispC pt Es Esdot') 

op5=[pt, Inert ial(pt,8), Inert iaKpt, 9)]; 
disp(op5) 

dispC pt Ux Uh ’ ) 

op4*[pt, Inert ial (pt,5), Inert iaKpt, 6)]; 
dlsp(op4) 

* 

Xaircraft control: 

[theta, T]»CRLT (q, de I tat , theta, a Ipha, a I phaeax, y , check 1 , 
ydot,T,U) 

if pt«*iters, Xnunber of iterations 

theta»999; 
el seif y C 1 ) >5000, 
theta*999; 
end 

end lehlle a Igor I the loop 

X 

Xoutput files: 

ansier"lnput( 'Do you eant a flight path plot? I»VES '); 
if ans«er»«l, 

p I ot ( R I rcraf t ( : , 1 ) , fl I rcra f t ( : , 2) ) 
pause 

ans»er=0.0; 
end 

ansfer*input( 'Do you eant x, h, U, and airspeed accel? 1«VES '); 
if ans*er s *1, 

dlepC pt x h U 

airspeed accel ') 

opl*[flircraft( : ,7), Aircraft ( : , 1 ),Aircraft( : ,2),Rircraft( : ,3), 
Rlrcraft( : ,4)]; 
dlsp(op1 ) 
else 
end 
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answer* input ('Do you want theta, alpha, gawaa, and gamai? 1*YES '); 
if answer**!, 

diepC pt theta alpha gawaa 

gawai ') 

op2*[Rircraft ( : ,7) , Aircraft ( : ,5),Rircraft( : ,6), Inert ial ( : ,3), 

Inert ial ( : ,4)]; 
disp(op2) 
else 
end 

answer* input ( 'Do you want Ug, ROC, and Thrust? 1*YES '); 
if answer**!, 

dispC pt Ug ROC Thrust') 

op3*[Rlrcraft( : ,?), Inert ial ( : , 1), Inert ial (.• ,2),flircraft(: ,8)]; 

disp(op3) 

else 

end 

answer* input ( 'Do you want Es and Esdot? I*YES '); 
if answer** 1, 

dispC pt Es Esdot') 

op5*[fllrcraft( : ,7), Inert Ial ( : ,8), Inert ial ( : ,9)]; 

disp(opS) 

else 

end 

answer* input ('Do you want Ux and Uh? 1*YES '); 
if answer**1, 

dispC pt Ux Uh') 

op4*[R I rcraft ( : , 7) , I nert i a I ( : , 5) , Inert i a I ( : , 6) ] ; 

disp(opl) 

else 

end 
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B. AIRCRAFT PERFORMANCE PARAMETER FUNCTION 



funct ion[S,K,CDo,Ut,CLo,CLalpha,alphaiax,q,T,theta] s P3onsta(U) 
XP3 at 120, 000 lbs, flaps and gear up. 

S=1300; Xplane for® surface area In ft ^2 

<■0.05041; ^coefficient for CD calculation 

CDo=0.0213; *coeff Iclent of drag 

Ut B 120000; laircraft weight in lbs 

CLo»0.350; Xzero flOfl lift coefficient 

CLalpha B 6.25; JMift curve slope 

alpha*ax B 0.209; Xstall buffet alpha •/ flaps up 

q B 0.0873; Xnominal pitch rate of 5deg/sec 

T- . 5*0 . 002377*lT2*S*(CDo+K*(Ut/( . 5* . 002377*U A 2*S) ) A 2) ; 

Ithrust in Ibf 



theta B (Ut/( . 5* . 002377*r2*S)-CLo)/CLa I pha ; 



Xvater 

angle 



I Ine deck 
in rad 
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C. FLIGHT PATH CONTROL FUNCTION 



fund lon[thetaout ,Tout ] B CflLT(q,deltat jtheta.alphQjalphainaXjy, 

check1,ydot,T,U) 

Xflircraft control using constant altitude escape naneuuer. 
thetaout'theta; 

Tout°T; 

if y(2)<check1-5, 
if ydot(2)<=0, 

thetaout=theta+q/2 . 5*de I tat j 
end 

elself y(2)>check1+5, 
if ydot(2)>*0, 

thetaout B theta-q/2.5*deltat ; 
end 
end 

I f abs(ydot(2)*deltat)>abs(y(2)-check1 ), 
thetaout»thetaout-sign(ydot(2))*q/1*deltat j 
end 

if alpha>alphaaax, 
thetaout"theta-(alpha-alpha*ax); 
end 

If y(5)<=(U-67.5), 

Tout c T+deltat/2*33400; 

end 

if Tout >33400, 

Tout "33400; 
end 
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APPENDIX F 



APPROACH TO LANDING ENCOUNTER GRAPHIC DATA 

This appendix contains the calculated performance of different 
aircraft and weight combinations upon encountering a microburst 
windshear during an approach to landing. Each figure is for a 
particular aircraft performing a specified escape maneuver. Each 
figure contains four graphs depicting altitude, theta, alpha, airspeed, 
specific energy, and time. The abscissa in all graphs is x (distance 
from microburst center). 
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Figure FI. P-3 at 89,5001bs gross weight performing a 
constant airspeed escape. 
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Figure F2. P-3 at 89,5001bs gross weight performing a 
constant altitude escape. 
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Figure F3. P-3 at 89,5001bs gross weight performing a 
constant 5° theta escape. 
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Figure F4. P-3 at 89,5001bs gross weight performing a 
constant 8« theta escape. 
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Figure F5. P-3 at 89,5001bs gross weight performing a 
constant 10° theta escape. 
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Figure F6. P-3 at 89,5001bs gross weight performing a 
constant 15° theta escape. 



190 



time(s) 








x(ft) 



Figure F7. P-3 at 89,5001bs gross weight performing a 
constant 12 unit AO A escape. 
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Figure F8. P-3 at 89,5001bs gross weight performing a 
constant 15 unit AOA escape. 
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Figure F9. P-3 at 89,5001bs gross weight performing a 
constant 20 unit AO A escape. 
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Figure F10. P-3 at 114,0001bs gross weight performing a 
constant airspeed escape. 
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Figure Fll. P-3 at 114,0001bs gross weight performing a 
constant altitude escape. 



195 



tim*(s) 







x(ft) 



Figure F12. P~3 at 114,0001bs gross weight performing a 
constant 5° theta escape. 
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Figure F13. P-3 at 114,0001bs gross weight performing a 
constant 8° theta escape. 
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Figure F14. P-3 at 114,0001bs gross weight performing a 
constant 10° theta escape. 
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Figure F15. P-3 at 114,0001bs gross weight performing a 
constant 15a theta escape. 
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Figure F16. P-3 at 114,0001bs gross weight performing a 
constant 12 unit AO A escape. 
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Figure F17. P-3 at 114,0001bs gross weight performing a 
constant 15 unit AOA escape. 



201 



tlme(s) 







x(ft) 



Figure F18. P-3 at 114,0001bs gross weight performing a 
constant 20 unit AOA escape. 



202 



tfme(s) 





20 

15 

* 10 
■$ 5 

0 



-500 0 500 1 000 1 500 2000 2500 3000 3500 4000 4500 5000 5500 



i j 1 j alpj 


k 1. L i | 


i j.'* fs'**. : j ..••j j 


[-—-j.. ! 


! T " \ ^ 




\ : yT j .7“ 


* 


I K ! 




y 1 i ; ; 




t — 1 — i — » — i — » — i — * — r 


— ■ — i — ■ — i 


i — ' — i — ■ — i — ■ — i — ■ — i — > — 1 





x(ft) 



Figure FI 9. Three-engine heavy airline transport 
performing a constant airspeed escape. 
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Figure F20. Three-engine heavy airline transport 
performing a constant altitude escape. 
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Figure F21. Three-engine heavy airline transport 
performing a constant 5° theta escape. 
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Figure F22. Three-engine heavy airline transport 
performing a constant 10° theta escape. 
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Figure F23. Three-engine heavy airline transport 
performing a constant 15° theta escape. 



207 






Es<«) Es(ft) h(ft) 







x(ft) 



Figure F24. T-44 performing a constant airspeed escape. 
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Figure F25. T-44 performing a constant altitude escape. 
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Figure F26. T-44 performing a constant 5° theta escape. 
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Figure F27. T-44 performing a constant 10° theta escape. 
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Figure F28. T-44 performing a constant 15° theta escape. 
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Figure F29. T-44 performing a constant 10 unit AO A 
escape. 
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Figure F30. T-44 performing a constant 16 unit AOA 
escape. 
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Figure F31. T-44 performing a constant 25 unit AO A 
escape. 
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Figure F32. P-3 at 89,5001bs gross weight with three 
engines operating performing a constant 5° theta escape. 
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Figure F33. P-3 at 89,5001bs gross weight with three 
engines operating performing a constant 10° theta escape. 
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Figure F34. P-3 at 89,5001bs gross weight with three 
engines operating performing a constant 15° theta escape. 
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APPENDIX G 



TAKEOFF ENCOUNTER GRAPHIC DATA 

This appendix contains the calculated performance of different 
aircraft and weight combinations upon encountering a microburst 
windshear at or immediately after lift off. Each figure is for a 
particular aircraft performing a specified escape maneuver. Each 
figure contains four graphs depicting altitude, theta, alpha, airspeed, 
specific energy, and time. The abscissa in all graphs is x (distance 
from microburst center). 
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Figure Gl. P-3 at 90,0001bs gross weight performing a 
constant airspeed maneuver during takeoff. 
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Figure G2. P-3 at 90,0001bs gross weight performing a 
constant 10o theta maneuver during takeoff. 
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Figure G3. P-3 at 90,0001bs gross weight performing a 
constant 12 unit AO A maneuver during takeoff. 
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Figure G4. P-3 at 90,0001bs gross weight performing a 
constant 15 unit AOA maneuver during takeoff. 
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Figure G5. P-3 at 120,0001bs gross weight performing a 
constant 10° theta maneuver during takeoff. 
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Figure G6. P-3 at 120,0001bs gross weight performing a 
constant 15° theta maneuver during takeoff. 
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Figure G7. P-3 at 120,0001bs gross weight with an increased 
rotate speed followed by a constant 10° theta climb. 
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APPENDIX H 



P-3 FLIGHT SIMULATOR DATA 
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FUEL 


FUEL 


ENGINE 


1 


2 


3 


4 


TANK 


LOAD 
















SHP 


483 


484 


499 


597 


1 


3937 
















TIT 


549 


549 


551 


564 


2 


3937 
















FUEL FLOU 


849 


849 


855 


893 


3 


3937 












4 


3935 












5 


996 


RUNUAY- 100 


MOFFETT 


NAS RUY 


32R 





BRAKES 

PILOT 

COPILOT 



FLIGHT 

IAS 115 
ALT 490 
HDG 265.9 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 



GROSS UEIGHT 
C.G. 

FLAP POSITION 
GEAR POSITION 



00: 00: 00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 



00: 03: 58 



89960 
20. 00 
19. 1 
1 . 0 



PRESSURE ALTITUDE 
CALIBRATED AIRSPD 
EQUIUALENT AIRSPD 
TRUE AIRSPD ( F/S ) 
MACH NUMBER 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH UELOCITY (D/S) 
ROLL UELOCITY (D/S) 
YAU UELOCITY (D/S) 
NORTH-SOUTH UELOCITY 
EAST-UEST UELOCITY 
UERTICAL UELOCITY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
UERTICAL ACCEL (G’S) 



FLIGHT/AERO 



524.0 
1 14. 7 
1 14.72 
195. 41 
0. 18 



ELEUATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 



7. 3 
6.4 
283.2 
-0.156 
-0. 164 
-0.273 
111.71 
-160.32 
-2.99 
0. 1216 
0.0117 
-0.9739 

CONTROL 

-5.43 

16.06 

-4.34 

-11.06 

-0.30 

- 2.20 



BANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM ) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
UERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 



-1.4 
-1 . 4 
181 

0. 0004 
-0. 0024 
-0. 0023 
-0. 373 
-0. 138 
0.579 
336 
-1664 
-3575 



LOADING 



ELEUATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TAB 
PEDAL POSITION 
AILERON TRIM TAB 
UHEEL POSITION 



12.69 
3.24 
3.79 
-0. 47 
0.42 
-1 . 16 



TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 



ENGINES 



I XX 
122 



INERTIA 

INERTIA 



(/ 

(/ 



1024) 

1024) 



5024 

46.4 

460 

UEIGHT 

725 

1549 



AND 



THRUST COEFFICIENT 
LATERAL T.C. . 

ENGINE T. I . T. 

BALANCE 

IYY INERTIA (/ 1024) 
CROSS PRODUCT/ INERT I A 



0. 09 
0. 08 
549 



862 

43240 



NOTE: UALUES INUALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRl 
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FUEL FUEL 
TANK LOAD 

1 3929 

2 3929 

3 3929 

4 3927 

5 996 



ENG I NE 


1 


2 


3 


4 


BRAKES 




SHP 


485 


486 


496 


600 


PILOT 




TIT 


550 


550 


551 


564 


COPILOT 


FUEL FLOU 


848 


848 


853 


893 


FLIGHT 






IAS 


108 












ALT 


490 


RUNUAY- 100 


MOFFETT 


NAS RUY 


32R 




HDG 


256. 8 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 

GROSS UEIGHT 
C G 

FLAP POSITION 
GEAR POSITION 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH UELOCITY (D/S) 
ROLL UELOCITY (D/S) 
YAU UELOCITY (D/S) 
NORTH-SOUTH UELOCITY 
EAST-UEST UELOCITY 
UERTICAL UELOCITY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
UERTICAL ACCEL (G’S) 

ELEUATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 

TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 

IXX INERTIA (/ 1024) 
122 INERTIA (/ 1024) 



00:00:00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 



89928 
20. 00 
19. 1 
1 . 0 



FLIGHT/ AERO 



11.9 
7.8 
273. 8 
0. 195 
0. 086 
-0. 023 
80. 19 
-166. 20 
-13.38 
0. 1999 
0. 0046 
-0. 9636 

CONTROL 

-6.40 
16.87 
-4.98 
-14.28 
0. 44 



LOADING 



0.03 

5232 
46. 4 
461 

UEIGHT 

725 

1549 



PRESSURE ALTITUDE 
CALIBRATED AIRSPD 
EQUIUALENT AIRSPD 
TRUE AIRSPD (F/S) 

MACH NUMBER 

BANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM ) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
UERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 

ELEUATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TAB 
PEDAL POSITION 
AILERON TRIM TAB 
UHEEL POSITION 



ENGINES 



THRUST COEFFICIENT 
LATERAL T.C. . 

ENGINE T. I . T. 

AND BALANCE 

IYY INERTIA (/ 1024) 
CROSS PRODUCT/ INERT I A 



00: 04:27 

524.0 
108. 6 
108. 59 
185. 00 
0. 17 

-0.3 
-1 . 0 
804 

- 0. 0002 
0 . 0003 
-0. 0007 
-0. 056 
0. 108 
0.455 
-176 
256 
-1092 

12.69 

2.76 

3.79 

-0.55 

0.42 

1.46 

0 . 10 
0 . 10 
550 

862 

43230 



NOTE: UALUES INUALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRUE 



FUEL 


FUEL 


ENGINE 


1 


2 


3 


4 


BRAKES 




TANK 


LOAD 


SHP 


486 


4B7 


505 


604 


PILOT 




1 


3923 


TIT 


549 


549 


551 


564 


COPILOT 




3923 
















2 


FUEL FLOU 


850 


850 


856 


894 


FLIGHT 




3 


3923 












IAS 


1 17 


4 


3921 












ALT 


490 
















5 


996 


RUNUAY- 100 


MOFFETT 


NAS RUY 


32R 




HDG 


256. 2 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 



GROSS UEIGHT 
C G 

FLAP POSITION 
GEAR POSITION 



00:00:00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 



00:04:49 



89904 
20. 00 
19. 1 
1.0 



PRESSURE ALTITUDE 
CALIBRATED AIRSPD 
EQUIUALENT AIRSPD 
TRUE AIRSPD (F/S) 
MACH NUMBER 



524.0 
1 17.4 
1 17.37 
200. 03 
0. 18 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH UELOCITY (D/S) 
ROLL UELOCITY (D/S) 
YAU UELOCITY (D/S) 
NORTH-SOUTH UELOCITY 
EAST-UEST UELOCITY 
UERTICAL UELOCITY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
UERTICAL ACCEL (G’S) 



ELEUATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 



FLIGHT/AERO 

2.4 
6.4 
273. 3 
-0. 016 
0. 000 
0. 023 
86 . 1 1 
-1B0. 05 
13. 63 
0.0376 
0.0046 
-1.0149 

CONTROL LOADING 



BANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM ) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
UERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 



- 0.2 
- 0.6 
-821 
— 0 . 0000 
-0 . 0008 
- 0. 0000 
-0.025 
0.210 
-0.4B2 
-176 
-640 
-16 



-5.43 
15.94 
-3.51 
-8.97 
0. 19 
-1.52 



TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 



ENGINES 



ELEUATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TAB 
PEDAL POSITION 
AILERON TRIM TAB 
UHEEL POSITION 



12.69 

3.25 

3.79 

-0.3B 

0.42 

0.52 



4B4B 

46.4 

463 



IXX INERTIA (/ 1024) 
IZZ INERTIA (/ 1024) 



THRUST COEFFICIENT 
LATERAL T.C. . 

ENGINE T. I . T. 

UEIGHT AND BALANCE 

725 IYY INERTIA (/ 1024) 

1549 CROSS PRODUCT/ INERT I A 



0.08 

0.08 

549 



861 

43222 



NOTE: UALUES INUALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRU 
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FUEL 


FUEL 


ENGINE 


1 


2 


3 


4 


BRAKES 




TANK 


LOAD 


SHP 


483 


484 


499 


597 


PILOT 




1 


3983 


TIT 


550 


550 


552 


564 


COPILOT 




3983 
















2 


FUEL FLOU 


849 


849 


855 


893 


FLIGHT 




3 


3983 












IAS 


1 14 


4 


3983 












ALT 


490 
















5 


997 


RUNUAY- 100 


MOFFETT 


NAS RUY 


32R 




HDG 


321.6 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 

GROSS UEIGHT 
C. G. 

FLAP POSITION 
GEAR POSITION 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH UELOCITY (D/S) 
ROLL UELOCITY (D/S) 
YAU UELOCITY (D/S) 
NORTH-SOUTH UELOCITY 
EAST-UEST UELOCITY 
UERTICAL UELOCITY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
UERTICAL ACCEL (G’S) 

ELEUATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 

TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 

IXX INERTIA (/ 1024) 
IZ2 INERTIA (/ 1024) 



00:00:00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 
90144 PRESSURE ALTITUDE 

24.65 CALIBRATED AIRSPD 

19. 1 EQUIUALENT AIRSPD 

1.0 TRUE AIRSPD ( F/S ) 

MACH NUMBER 

FLIGHT/AERO 



7. 7 
6.6 
338.8 
0. 250 
-0.141 
-0. 359 
195. 01 
5. 56 
-4.21 
0. 1284 
- 0.0010 
-0.9914 

CONTROL LOADING 



BANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM ) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
UERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 



-3.96 
8.94 
-3. 16 
-5. 06 
0.68 
-0. 72 



ELEUATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TAB 
PEDAL POSITION 
AILERON TRIM TAB 
UHEEL POSITION 



ENGINES 

5056 THRUST COEFFICIENT 

46.4 LATERAL T.C. 

460 ENGINE T. I . T.* 

UEIGHT AND BALANCE 

728 IYY INERTIA (/ 1024) 

1551 CROSS PRODUCT/ INERT I A 



00:01:13 

524. 0 
114.5 
114. 53 
195. 12 
0. 18 



- 1.8 

- 0.2 

252 

0. 0025 
-0. 0007 
0.001 1 
-0. 238 
-1.073 
-0.013 
2128 
-576 
1837 

12.69 
4.05 
3. 79 
-0.34 
0.42 
2. 15 

0.09 
0. 08 
550 

862 

43270 



NOTE: UALUES INUALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRUE 
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FUEL 


FUEL 


ENGINE 


1 


2 


3 


4 


BRAKES 




TANK 


LOAD 


SHP 


484 


483 


500 


597 


PILOT 




1 


3917 


TIT 


549 


549 


551 


564 


COPILOT 




3917 


















FUEL FLOU 


849 


849 


855 


893 


FLIGHT 




3 


3917 












IAS 


1 14 


4 


3915 












ALT 


490 
















5 


996 


RUNUAY- 100 


MOFFETT 


NAS RUY 


32R 




HDG 


255.8 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 

GROSS UEIGHT 
C G 

FLAP POSITION 
GEAR POSITION 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH VELOCITY (D/S) 
ROLL VELOCITY (D/S) 
YAU VELOCITY (D/S) 
NORTH-SOUTH VELOCITY 
EAST-UEST VELOCITY 
VERTICAL VELOCITY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
VERTICAL ACCEL (G’S) 

ELEVATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 

TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 

I XX INERTIA (/ 1024) 
IZ2 INERTIA (/ 1024) 



00:00:00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 



89880 
30.00 
19. 1 
1.0 

FLIGHT/AERO 

8 . 1 
6.2 
272.9 
-0.039 
-0. 086 
0.062 
84. 02 
-176. 15 
-6.78 
0. 1567 
0. 0002 
-0.9844 

CONTROL LOADING 

-1.54 
1 . 19 
-3. 10 
-5.31 
0. 19 
-1 . 14 



ENGINES 



5056 

46.4 

460 



PRESSURE ALTITUDE 
CALIBRATED AIRSPD 
EQUIVALENT AIRSPD 
TRUE AIRSPD ( F/S ) 

MACH NUM8ER 

BANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM ) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
VERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 

ELEVATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TAB 
PEDAL POSITION 
AILERON TRIM TAB 
UHEEL POSITION 

THRUST COEFFICIENT 
LATERAL T.C. 

ENGINE T. I.T.* 



UEIGHT AND BALANCE 



725 

1548 



IYY INERTIA (/ 1024) 
CROSS PRODUCT/ INERT I A 



00: 05: 10 

524.0 
114.7 
1 14. 72 
195.28 
0. 18 

0.4 
-0. 3 
407 

0.0028 
- 0. 001 1 
-0.0014 
0.438 
0.376 
0. 123 
2384 
-768 
-2103 

12.55 

5.33 

3.79 

-0.34 

0.42 

0.53 

0. 09 
0.08 
549 

861 

43214 



NOTE: VALUES INVALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRI 
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FUEL 


FUEL 


ENGINE 


1 


2 


3 


4 


BRAKES 




TANK 


LOAD 


SHP 


486 


485 


496 


601 


PILOT 




1 


3911 


TIT 


550 


550 


551 


564 


COP I LOT 




391 1 
















2 


FUEL FLOU 


848 


848 


853 


893 


FLIGHT 




3 


3910 












IAS 


108 


4 


3908 












ALT 


490 
















5 


996 


RUNUAY- 100 


MOFFETT 


NAS RUY 


32R 




HDG 


252. 4 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 

GROSS UEIGHT 
C.G. 

FLAP POSITION 
GEAR POSITION 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH UELOC I TV (D/S) 
ROLL UELOC I TY (D/S) 
YAU UELOC I TY (D/S) 
NORTH-SOUTH UELOC I TY 
EAST-UEST UELOC I TY 
UERTICAL UELOC I TY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
UERTICAL ACCEL (G’S) 

ELEUATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 

TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 

IXX INERTIA (/ 1024) 
122 INERTIA (/ 1024) 



00:00:00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 



89856 
30.00 
19. 1 
1.0 

FLIGHT/AERO 

12.8 
7.3 
269. 5 
-0.109 
0.070 
0. 062 
71.96 
-169.59 
-18.07 
0 . 2229 
- 0.0112 
-0. 9587 

CONTROL LOADING 

-1.72 
1.41 
-2.91 
-4.37 
1.38 
-1. 17 



ENGINES 



5168 

46.4 

463 



PRESSURE ALTITUDE 
CALIBRATED AIRSPD 
EQUIUALENT AIRSPD 
TRUE AIRSPD (F/S) 

MACH NUMBER 

BANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
UERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 

ELEUATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TAB 
PEDAL POSITION 
AILERON TRIM TAB 
UHEEL POSITION 

THRUST COEFFICIENT 
LATERAL T.C. 

ENGINE T. I.T.* 



UEIGHT AND BALANCE 



724 

1548 



IYY INERTIA (/ 1024) 
CROSS PRODUCT/ INERT I A 



00: 05: 32 



524.0 
108. 6 
108. 59 
185.09 
0. 17 

1 . 1 
0.7 
1085 
0. 0002 
-0. 0003 
- 0. 0002 
0. 255 
-0. 033 
0.483 
80 

-192 

-276 

12.42 
5.23 
3. 79 
-0.32 
0.42 
4.28 

0. 10 
0. 10 
550 

861 

43206 



NOTE: UALUES INUALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRUE 



FUEL 


FUEL 


ENGINE 


1 


2 


3 


4 


TANK 


LOAD 


SHP 


4744 


4744 


4744 


4744 


1 


1 0000 


TIT 


1075 


1075 


1075 


1075 


2 


10000 


FUEL FLOU 


2467 


2466 


2466 


2467 


3 


10000 












4 


10000 












5 


995 


RUNUAY- 100 


MOFFETT 1 


NAS RUY 


32R 





BRAKES 

PILOT 

COPILOT 



FLIGHT 

IAS 125 
ALT 490 
HDG 22B.6 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 



00:00:00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 



00: 07: 52 



GROSS UEIGHT 
C.G. 

FLAP POSITION 
GEAR POSITION 



1 14200 
26.01 
19. 1 
1.0 



PRESSURE ALTITUDE 
CALIBRATED AIRSPD 
EQUIUALENT AIRSPD 
TRUE AIRSPD (F/S) 
MACH NUMBER 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH UELOCITY (D/S) 
ROLL UELOCITY (D/S) 
YAU UELOCITY (D/S) 
NORTH-SOUTH UELOCITY 
EAST-UEST UELOCITY 
UERTICAL UELOCITY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
UERTICAL ACCEL (G'S) 



ELEUATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 



FLIGHT/AERO 

8.0 
4.6 
246. 1 
0.273 
-0.430 
-0.500 
-5.B2 
-213.54 
-12.72 
0. 1443 
-0.0097 
-1 . 0163 

CONTROL LOADING 



524.0 
125.6 
125.59 
213.97 
0. 19 



BANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM ) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
UERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 



-2.7 

- 0.8 

759 

0. 0000 
0.0010 
-0.0048 
-1.923 
0. 015 
-0.B01 
-176 
1344 
-9720 



-2.7B 

6.28 

-10.06 

-46.25 

1.29 

-1.56 



ELEUATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TAB 
PEDAL POSITION 
AILERON TRIM TAB 
UHEEL POSITION 



10.69 
4.66 
5. 0B 
- 1 . 12 
0.42 
4.04 



TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 



ENGINES 



292B0 

90.0 

4699 



THRUST COEFFICIENT 
LATERAL T.C. . 
ENGINE T. I . T. 



0.42 
0. 42 
1075 



IXX INERTIA (/ 1024) 
I2Z INERTIA (/ 1024) 



UEIGHT AND BALANCE 



1 154 
1984 



IYY INERTIA (/ 1024) 
CROSS PRODUCT/ INERT I A 



872 

47766 



NOTE: UALUES INUALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRU 
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FUEL 


FUEL 


ENG I NE 


1 


2 


3 


4 


BRAKES 




TANK 


LOAD 


SHP 


4720 


4722 


4722 


4720 


PILOT 




1 


1 0000 


TIT 


1075 


1075 


1075 


1075 


COPILOT 




1 0000 
















2 


FUEL FLOU 


2459 


2459 


2459 


2459 


FLIGHT 




3 


1 0000 












IAS 


113 


4 


1 0000 












ALT 


490 
















S 


995 


RUNUAY- 100 


MOFFETT 1 


NAS RUY 


32R 




HDG 


203. 1 



FLIGHT CONDITIONS PAGE 



FLIGHT TIMER 

GROSS UEIGHT 
C G 

FLAP POSITION 
GEAR POSITION 



PITCH ANGLE 
ANGLE OF ATTACK 
HEADING ANGLE 
PITCH VELOCITY (D/S) 
ROLL VELOCITY (D/S) 
YAU VELOCITY (D/S) 
NORTH-SOUTH VELOCITY 
EAST-UEST VELOCITY 
VERTICAL VELOCITY 
LONGITUDINAL ACCEL 
LATERAL ACCEL 
VERTICAL ACCEL (G’S) 

ELEVATOR POSITION 
COLUMN FORCE 
RUDDER POSITION 
PEDAL FORCE 
AILERON POSITION 
UHEEL FORCE 

TOTAL THRUST 
THROTTLE ANGLE 
ENGINE S.H.P. 

IXX INERTIA (/ 1024) 
I2Z INERTIA (/ 1024) 



00:00: 00 MET TIMER 

CONF I GURAT I ON/COND I T I ONS 



1 14208 
26. 01 
19. 1 
1 . 0 



FLIGHT/AERO 



12. 9 
5. 9 
220. 1 
-0.29? 

0. 102 
0. 039 
-92. 79 
-166.59 
-23.70 
0.2163 
0.0105 
-0.9460 

CONTROL 

-3. 15 
8.00 
-12.70 
-48.78 
0. 75 
-0.83 



LOADING 



PRESSURE ALTITUDE 
CALIBRATED AIRSPD 
EQUIVALENT AIRSPD 
TRUE AIRSPD (F/S) 

MACH NUM8ER 

8ANK ANGLE 
SIDESLIP 

RATE OF CLIMB ( FPM ) 
PITCH ACCELERATION 
ROLL ACCELERATION 
YAU ACCELERATION 
NORTH-SOUTH ACCEL 
EAST-UEST ACCELERATION 
VERTICAL ACCELERATION 
TOTAL PITCHING MOMENT 
TOTAL ROLLING MOMENT 
TOTAL YAUING MOMENT 

ELEVATOR TRIM TAB 
COLUMN POSITION 
RUDDER TRIM TA8 
PEDAL POSITION 
AILERON TRIM TA8 
UHEEL POSITION 



ENGINES 



30176 THRUST COEFFICIENT 

90.0 LATERAL T.C. 

4671 ENGINE T. I. T. 

UEIGHT AND BALANCE 

1154 IYY INERTIA (/ 1024) 

1984 CROSS PRODUCT/ INERT I A 



00: 08: 32 

524.0 
1 12. 7 
112.72 
192. 12 
0 . 1 ? 

-0.5 
-2.3 
1425 
0.0022 
-0.0014 
-0.0004 
-0.024 
-0. 234 
0.929 
1872 
-1728 
-710 

10. 69 
4.45 
5.08 
-1.41 
0.42 
2.34 

0.54 

0.54 

1075 

872 

47766 



NOTE: VALUES INVALID DURING ATG - TO USE COL MARKER SU FOR SNAPS SET COLSNP TRUE 
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FUEL FUEL 




ENGINE 


I 


2 3 4 


BRAKES 


TANK LOAD 


















SHP 


4743 


4742 4742 4743 


PILOT 


1 10000 


















TIT 


1075 


1075 1075 1075 


COPILOT 


2 10000 


















FUEL FLOU 


2465 


2465 2465 2465 


FLIGHT 


3 10000 
























IAS 


124 


4 10000 
























ALT 


490 


5 995 


RUNUAY- 100 MOFFETT NAS RUY 32R 
















HDG 


193. 7 






FLIGHT 


CONDITIONS PAGE 






FLIGHT TIMER 




00 ! 00 ! 00 




MET TIMER 


00: 09 : 06 






CONF I GURAT I ON/COND I T I ONS 






GROSS UEIGHT 




1 14208 




PRESSURE ALTITUDE 


524.0 


C.G. 




21 . 00 




CALIBRATED AIRSPD 


124.2 


FLAP POSITION 




19. 1 




EQUIUALENT AIRSPD 


124. 16 


GEAR POSITION 




1.0 




TRUE AIRSPD ( F/S ) 


211.66 










MACH NUMBER 




0. 19 






FLIGHT/AERO 






PITCH ANGLE 




8.5 




BANK ANGLE 




-1.2 


ANGLE OF ATTACK 




5. 1 




SIDESLIP 




-4. 1 


HEADING ANGLE 




210. 7 




RATE OF CLIMB 


( FPM ) 


757 


PITCH UELOCITY (D/S) 


0. 086 




PITCH ACCELERATION 


-0. 0016 


ROLL UELOCITY (D/S) 


0.250 




ROLL ACCELERATION 


-0.0037 


YAU UELOCITY (D/S) 


-0. 180 




YAU ACCELERATION 


-0. 0048 


NORTH-SOUTH UELOCITY 


-136.83 




NORTH-SOUTH ACCEL 


0.518 


EAST-UEST UELOCITY 


-160.99 




EAST-UEST ACCELERATION 


-0.250 


UERTICAL UELOCITY 


-12. 68 




UERTICAL ACCELERATION 


-0.932 


LONGITUDINAL ACCEL 


0. 1521 




TOTAL PITCHING 


MOMENT 


-1456 


LATERAL ACCEL 




0.0393 




TOTAL ROLLING 


MOMENT 


-4160 


UERTICAL ACCEL (G’S) 


-1.0176 




TOTAL YAUING MOMENT 


-9577 






CONTROL LOADING 






ELEUATOR POSITION 


-4.99 




ELEUATOR TRIM 


TAB 


10.69 


COLUMN FORCE 




16.50 




COLUMN POSITION 


3.46 


RUDDER POSITION 




-14.61 




RUDDER TRIM TAB 


8. 98 


PEDAL FORCE 




-52.44 




PEDAL POSITION 


-1.62 


AILERON POSITION 




0. 17 




AILERON TRIM TAB 


0.42 


UHEEL FORCE 




-1.41 




UHEEL POSITION 


0.45 








ENGINES 






TOTAL THRUST 




29376 




THRUST COEFFICIENT 


0. 43 


THROTTLE ANGLE 




90.0 




LATERAL T.C. . 




0.43 


ENGINE S.H.P. 




4693 




ENGINE T. I.T. 




1075 






UEIGHT AND BALANCE 






IXX INERTIA (/ 1024) 


1 154 




IYY INERTIA (/ 


1024) 


872 


IZZ INERTIA (/ 1024) 


1984 




CROSS PRODUCT/ INERTIA 


47766 


NOTE: UALUES INUALID 


DURING ATG 


- TO USE 


COL MARKER SU FOR 


SNAPS SET 


COLSNP TRUE 
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FUEL FUEL 




ENGINE 


12 3 4 


8RAKES 


TANK LOAD 
















SHP 


4722 4721 4721 4722 


PILOT 


1 10000 
















TIT 


1075 1075 1075 1075 


COPILOT 


2 10000 
















FUEL FLOU 


2459 2459 2459 2459 


FLIGHT 


3 10000 




















IAS 


1 13 


>4 10000 




















ALT 


490 


5 995 


RUNUAY- 100 MOFFETT NAS RUY 32R 














HDG 


171 . 9 






FLIGHT 


CONDITIONS PAGE 






FLIGHT TIMER 




00 • 00 : 00 


MET TIMER 


00: 10: 07 






CONF I GURAT I ON/COND I T I ONS 






GROSS UEIGHT 




1 14208 


PRESSURE ALTITUDE 


524.0 


C.G. 




21 . 00 


CALIBRATED AIRSPD 


1 12. 6 


FLAP POSITION 




19. 1 


EQUIUALENT AIRSPD 


1 12.56 


GEAR POSITION 




1.0 


TRUE AIRSPD ( F/S ) 


191.91 








MACH NUM8ER 




0. 17 






FLIGHT/AERO 






PITCH ANGLE 




12.6 


BANK ANGLE 




-0. 7 


ANGLE OF ATTACK 




6.8 


SIDESLIP 




-2. 5 


HEADING ANGLE 




188. 7 


RATE OF CLIM8 


( FPM ) 


1 169 


PITCH UELOCITY (D/S) 


-0.070 


PITCH ACCELERATION 


-0. 0010 


ROLL UELOCITY (D/S) 


0.273 


ROLL ACCELERATION 


-0.0012 


YAU UELOCITY (D/S) 


0.055 


YAU ACCELERATION 


0. 0032 


NORTH-SOUTH UELOCITY 


-166.58 


NORTH-SOUTH ACCEL 


-0.232 


EAST-UEST UELOCITY 


-93.36 


EAST-UEST ACCELERATION 


-0. 130 


UERTICAL UELOCITY 


-19.51 


UERTICAL ACCELERATION 


-0.237 


LONGITUDINAL ACCEL 


0. 2287 


TOTAL PITCHING 


MOMENT 


-944 


LATERAL ACCEL 




0.0112 


TOTAL ROLLING 


MOMENT 


-1600 


UERTICAL ACCEL (G’S) 


-0.9812 


TOTAL YAUING MOMENT 


6515 






CONTROL LOADING 






ELEUATOR POSITION 


-6.02 


ELEUATOR TRIM 


TA8 


10.69 


COLUMN FORCE 




15.87 


COLUMN POSITION 


2.96 


RUDDER POSITION 




-14. 19 


RUDDER TRIM TAB 


8.98 


PEDAL FORCE 




-41.87 


PEDAL POSITION 




-1.56 


AILERON POSITION 




0.90 


AILERON TRIM TA8 


0.42 


UHEEL FORCE 




1.61 


UHEEL POSITION 


3. 13 








ENGINES 






TOTAL THRUST 




30208 


THRUST COEFFICIENT 


0.54 


THROTTLE ANGLE 




90.0 


LATERAL T.C. . 




0.54 


ENGINE S.H.P. 




4674 


ENGINE T. I . T. 




1075 






UEIGHT AND BALANCE 






IXX INERTIA (/ 1024) 


1 154 


IYY INERTIA (/ 


1024) 


872 


IZZ INERTIA (/ 1024) 


1984 


CROSS PRODUCT/ INERT I A 


47766 


NOTE: UALUES INUALID 


DURING ATG • 


- TO USE COL MARKER SU FOR 


SNAPS SET 


COLSNP TRUE 
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